Search

Utilizing localized orbitals, local correlation theory can reduce the unphysically high system-size scaling of post-Hartree–Fock (post-HF) methods to linear scaling in insulating molecules. The sparsity of the four-index electron repulsion integral (ERI) tensor is central to achieving this reduction. For second-order Møller–Plesset theory (MP2), one of the simplest post-HF methods, only the (ia|jb) ERIs are needed, coupling occupied orbitals i, j and virtuals a, b. In this paper, we compare the numerical sparsity (called the "ragged list") and two other approaches revealing the low-rank sparsity of the ERI. The ragged list requires only one set of (localized) virtual orbitals, and we find that the orthogonal valence virtual-hard virtual set of virtuals originally proposed by Subotnik et al. gives the sparsest ERI tensor. To further compress the ERI tensor, the pair natural orbital (PNO) type representation uses different sets of virtual orbitals for different occupied orbital pairs, while the occupied-specific virtual (OSV) approach uses different virtuals for each occupied orbital. Our results indicate that while the low-rank PNO representation achieves significant rank reduction, it also requires more memory than the ragged list. The OSV approach requires similar memory to that of the ragged list, but it involves greater algorithmic complexity. An approximation (called the "fixed sparsity pattern") for solving the local MP2 equations using the numerically sparse ERI tensor is proposed and tested to be sufficiently accurate and to have highly controllable error. A low-scaling local MP2 algorithm based on the ragged list and the fixed sparsity pattern is therefore promising. [ABSTRACT FROM AUTHOR]

The excitation energy transfer (EET) process for photosynthetic antenna complexes consisting of subunits, each comprised of multiple chromophores, remains challenging to describe. The multichromophoric Förster resonance energy transfer theory is a popular method to describe the EET process in such systems. This paper presents a new time-domain method for calculating energy transfer based on the combination of multichromophoric Förster resonance energy transfer theory and the Numerical Integration of the Schrödinger Equation method. After validating the method on simple model systems, we apply it to the Light-Harvesting antenna 2 (LH2) complex, a light harvesting antenna found in purple bacteria. We use a simple model combining the overdamped Brownian oscillators to describe the dynamic disorder originating from the environmental fluctuations and the transition charge from the electrostatic potential coupling model to determine the interactions between chromophores. We demonstrate that with this model, both the calculated spectra and the EET rates between the two rings within the LH2 complex agree well with experimental results. We further find that the transfer between the strongly coupled rings of neighboring LH2 complexes can also be well described with our method. We conclude that our new method accurately describes the EET rate for biologically relevant multichromophoric systems, which are similar to the LH2 complex. Computationally, the new method is very tractable, especially for slow processes. We foresee that the method can be applied to efficiently calculate transfer in artificial systems as well and may pave the way for calculating multidimensional spectra of extensive multichromophoric systems in the future. [ABSTRACT FROM AUTHOR]

Over the last century, quantum theories have revolutionized our understanding of material properties. One of the most striking quantum phenomena occurring in heterogeneous media is the quantum tunneling effect, where carriers can tunnel through potential barriers even if the barrier height exceeds the carrier energy. Interestingly, the tunneling process can be accompanied by the absorption or emission of light. In most tunneling junctions made of noble metal electrodes, these optical phenomena are governed by plasmonic modes, i.e., light-driven collective oscillations of surface electrons. In the emission process, plasmon excitation via inelastic tunneling electrons can improve the efficiency of photon generation, resulting in bright nanoscale optical sources. On the other hand, the incident light can affect the tunneling behavior of plasmonic junctions as well, leading to phenomena such as optical rectification and induced photocurrent. Thus, plasmonic tunneling junctions provide a rich platform for investigating light–matter interactions, paving the way for various applications, including nanoscale light sources, sensors, and chemical reactors. In this paper, we will introduce recent research progress and promising applications based on plasmonic tunneling junctions. [ABSTRACT FROM AUTHOR]

Introducing nanorods into a polymeric matrix can enhance the physical and mechanical properties of the resulting material. In this paper, we focus on understanding the dispersion and orientation patterns of nanorods in an unentangled polymer melt, particularly as a function of nanorod concentration, using molecular dynamics simulations. The system is comprised of flexible polymer chains and multi-thread nanorods that are equilibrated in the NPT ensemble. All interactions are purely repulsive except for those between polymers and rods. Results with attractive vs repulsive polymer–rod interactions are compared and contrasted. The concentration of rods has a direct impact on the phase behavior of the system. At lower concentrations, rods phase separate into nematic clusters, whereas at higher concentrations more isotropic and less structured rod configurations are observed. A detailed examination of the conformation of the polymer chains near the rod surface shows extension of the chains along the director of the rods (especially within clusters). The dispersion and orientation of the nanorods are a result of the competition between depletion entropic forces responsible for the formation of rod clusters, the enthalpic effects that improve mixing of rods and polymer, and entropic losses of polymers interpenetrating rod clusters. [ABSTRACT FROM AUTHOR]

Many phase change materials (PCMs) are found to crystallize without exhibiting a glass transition endotherm upon reheating. In this paper, we review experimental evidence revealing that these PCMs and likely other hyperquenched molecular and metallic systems can crystallize from the glassy state when reheated at a standard rate. Among these evidences, PCMs annealed below the glass transition temperature Tg exhibit slower crystallization kinetics despite an increase in the number of sub-critical nuclei that should promote the crystallization speed. Flash calorimetry uncovers the glass transition endotherm hidden by crystallization and reveals a distinct change in kinetics when crystallization switches from the glassy to the supercooled liquid state. The resulting Tg value also rationalizes the presence of the pre-Tg relaxation exotherm ubiquitous of hyperquenched systems. Finally, the shift in crystallization temperature during annealing exhibits a non-exponential decay that is characteristic of structural relaxation in the glass. Modeling using a modified Turnbull equation for nucleation rate supports the existence of sub-Tg fast crystallization and emphasizes the benefit of a fragile-to-strong transition for PCM applications due to a reduction in crystallization at low temperature (improved data retention) and increasing its speed at high temperature (faster computing). [ABSTRACT FROM AUTHOR]

Electrochemical impedance spectroscopy (EIS) is a powerful characterization technique for the in-depth investigation of kinetic/transport parameters detection, reaction mechanism understanding, and degradation effects exploration in lithium-ion battery (LIB) systems. However, due to the lack of standardized criterion/paradigm, severe misinterpretations occur frequently during an EIS measurement. In this paper, the significance of instrumental accuracy is described and the character/principle of selection on the simulation model is illuminated/proposed, showing that an adequate precision device and an appropriate fitting model are a prerequisite for a correct EIS analysis. Moreover, the drawbacks of conventional two-electrode EIS experiments for typical coin-type cells are rigorously pointed out by comparison with the ideal three-electrode configuration, where the real impedance information of the cathode would be masked by the sum of both the anode film resistance response and the unavoidable inductive loop signal. The three-electrode case enables efficient accurate observations on individual electrodes, thus facilitating abundant and useful information acquisition. Consequently, devices with a sufficient accuracy, rational simulation models, and advanced three-electrode cells are distinctly illustrated as standardized criterion/paradigm for EIS characterizations, which are essentially important for electrode and interface modifications in LIBs. [ABSTRACT FROM AUTHOR]

Fewest-switches surface hopping (FSSH) has emerged as one of the leading methods for modeling the quantum dynamics of molecular systems. While its original formulation was limited to adiabatic populations, the growing interest in the application of FSSH to coherent phenomena prompts the question of how one should construct a complete density matrix based on FSSH trajectories. A straightforward solution is to define adiabatic coherences based on wavefunction coefficients. In this paper, we demonstrate that inconsistencies introduced in the density matrix through such treatment may lead to a violation of positivity. We furthermore show that a recently proposed coherent generalization of FSSH results in density matrices that satisfy positivity while yielding improved accuracy throughout much (but not all) of parameter space. [ABSTRACT FROM AUTHOR]

The transport of particles through channels is of paramount importance in physics, chemistry, and surface science due to its broad real world applications. Much insight can be gained by observing the transition paths of a particle through a channel and collecting statistics on the lifetimes in the channel or the escape probabilities from the channel. In this paper, we consider the diffusive transport through a narrow conical channel of a Brownian particle subject to intermittent dynamics, namely, stochastic resetting. As such, resetting brings the particle back to a desired location from where it resumes its diffusive phase. To this end, we extend the Fick–Jacobs theory of channel-facilitated diffusive transport to resetting-induced transport. Exact expressions for the conditional mean first passage times, escape probabilities, and the total average lifetime in the channel are obtained, and their behavior as a function of the resetting rate is highlighted. It is shown that resetting can expedite the transport through the channel—rigorous constraints for such conditions are then illustrated. Furthermore, we observe that a carefully chosen resetting rate can render the average lifetime of the particle inside the channel minimal. Interestingly, the optimal rate undergoes continuous and discontinuous transitions as some relevant system parameters are varied. The validity of our one-dimensional analysis and the corresponding theoretical predictions is supported by three-dimensional Brownian dynamics simulations. We thus believe that resetting can be useful to facilitate particle transport across biological membranes—a phenomenon that can spearhead further theoretical and experimental studies. [ABSTRACT FROM AUTHOR]

Using an approach based on the diffusion analog of the Cattaneo–Vernotte differential model, we find the exact analytical solution to the corresponding time-dependent linear hyperbolic initial boundary value problem, describing irreversible diffusion-controlled reactions under Smoluchowski's boundary condition on a spherical sink. By means of this solution, we extend exact analytical calculations for the time-dependent classical Smoluchowski rate coefficient to the case that includes the so-called inertial effects, occurring in the host media with finite relaxation times. We also present a brief survey of Smoluchowski's theory and its various subsequent refinements, including works devoted to the description of the short-time behavior of Brownian particles. In this paper, we managed to show that a known Rice's formula, commonly recognized earlier as an exact reaction rate coefficient for the case of hyperbolic diffusion, turned out to be only its approximation being a uniform upper bound of the exact value. Here, the obtained formula seems to be of great significance for bridging a known gap between an analytically estimated rate coefficient on the one hand and molecular dynamics simulations together with experimentally observed results for the short times regime on the other hand. A particular emphasis has been placed on the rigorous mathematical treatment and important properties of the relevant initial boundary value problems in parabolic and hyperbolic diffusion theories. [ABSTRACT FROM AUTHOR]

In this paper, we investigate prenucleation and heterogeneous nucleation in the liquid Pb/solid Al system as an example of systems with large lattice misfit using molecular dynamics simulation. Solid Pb and Al have a large positive lattice misfit (f) of 18.2% along the densely packed [110] direction. This study reveals that prenucleation occurs at 600 K (an undercooling of 15 K), and a 2-dimensional (2D) ordered structure forms at the interface with a coincidence site lattice (CSL) between the first Pb and first Al layers. The CSL accommodates the major part of the f, and only a small residual lattice misfit (fr) of 1.9% remains. The formation of the CSL transforms the original substrate into a considerably potent nucleant, where the first Pb layer becomes the new surface layer of the substrate. At an undercooling of about 22 K, nucleation proceeds by merging 2D ordered structure through structural templating: the second Pb layer is epitaxial to the CSL Pb layer, the third Pb layer largely accommodates the fr, and the fourth Pb layer is a nearly perfect crystalline plane. Further analysis indicates that the interface with the CSL has a lower interfacial energy than with a cube-to-cube orientation relationship. For the first time, we established that the CSL was an effective mechanism to accommodate the f for systems with a large positive misfits. Heterogeneous nucleation is governed not by a single mechanism (misfit dislocations in Turnbull's model), but instead by various mechanisms depending on f. This study sheds new light on the atomistic mechanism of heterogeneous nucleation. [ABSTRACT FROM AUTHOR]

For many laser-based methods of material characterization and analysis, a tensor formulation of theory is necessary, especially in techniques that exploit nonlinear or chiral optics. The fundamental interactions that underpin such methods offer various levels of approach to theory, but the most rigorous often lead to equations of considerable complexity. To compute the values for individual material parameters frequently demands making assumptions of extreme simplicity, overly dependent on calculational method, yet still providing unsatisfactory results. A pragmatic and entirely rigorous symmetry-based approach to the irreducible tensorial structures circumvents many of these problems, securing reliable results and guiding the pathway to applications. Instead of focusing on individual tensor components, such an approach can rapidly determine the number of linearly independent quantities—and hence the number of operationally different setups necessary for full characterization. By such means, one can directly ascertain how variations of optical polarization and beam geometry can reliably capture the response of any material system. The use of an irreducible tensor method operates independently of any means that might be chosen to calculate material properties. It removes the need for common simplifying assumptions, such as the approximation of tensorial structure by a scalar representation, adoption of a two-state model, or disregarding near-resonance damping. It also obviates any dependence on a choice of simulation package or quantum-calculational software. In this paper, the principles are set down and illustrated by application to experiments of varying degrees of complexity, including interactions of growing significance in the realm of chiral nonlinear optics. Limitations of this approach are also critically assessed. [ABSTRACT FROM AUTHOR]

Coarse-grained (CG) models facilitate an efficient exploration of complex systems by reducing the unnecessary degrees of freedom of the fine-grained (FG) system while recapitulating major structural correlations. Unlike structural properties, assessing dynamic properties in CG modeling is often unfeasible due to the accelerated dynamics of the CG models, which allows for more efficient structural sampling. Therefore, the ultimate goal of the present series of articles is to establish a better correspondence between the FG and CG dynamics. To assess and compare dynamical properties in the FG and the corresponding CG models, we utilize the excess entropy scaling relationship. For Paper I of this series, we provide evidence that the FG and the corresponding CG counterpart follow the same universal scaling relationship. By carefully reviewing and examining the literature, we develop a new theory to calculate excess entropies for the FG and CG systems while accounting for entropy representability. We demonstrate that the excess entropy scaling idea can be readily applied to liquid water and methanol systems at both the FG and CG resolutions. For both liquids, we reveal that the scaling exponents remain unchanged from the coarse-graining process, indicating that the scaling behavior is universal for the same underlying molecular systems. Combining this finding with the concept of mapping entropy in CG models, we show that the missing entropy plays an important role in accelerating the CG dynamics. [ABSTRACT FROM AUTHOR]

Simulating the unitary dynamics of a quantum system is a fundamental problem of quantum mechanics, in which quantum computers are believed to have significant advantage over their classical counterparts. One prominent such instance is the simulation of electronic dynamics, which plays an essential role in chemical reactions, non-equilibrium dynamics, and material design. These systems are time-dependent, which requires that the corresponding simulation algorithm can be successfully concatenated with itself over different time intervals to reproduce the overall coherent quantum dynamics of the system. In this paper, we quantify such simulation algorithms by the property of being fully-coherent: the algorithm succeeds with arbitrarily high success probability 1 − δ while only requiring a single copy of the initial state. We subsequently develop fully-coherent simulation algorithms based on quantum signal processing (QSP), including a novel algorithm that circumvents the use of amplitude amplification while also achieving a query complexity additive in time t, ln(1/δ), and ln(1/ϵ) for error tolerance ϵ: Θ ‖ H ‖ | t | + ln (1 / ϵ) + ln (1 / δ) . Furthermore, we numerically analyze these algorithms by applying them to the simulation of the spin dynamics of the Heisenberg model and the correlated electronic dynamics of an H2 molecule. Since any electronic Hamiltonian can be mapped to a spin Hamiltonian, our algorithm can efficiently simulate time-dependent ab initio electronic dynamics in the circuit model of quantum computation. Accordingly, it is also our hope that the present work serves as a bridge between QSP-based quantum algorithms and chemical dynamics, stimulating a cross-fertilization between these exciting fields. [ABSTRACT FROM AUTHOR]

Efficient energy transport in photosynthetic antenna is a long-standing source of inspiration for artificial light harvesting materials. However, characterizing the excited states of the constituent chromophores poses a considerable challenge to mainstream quantum chemical and semiempirical excited state methods due to their size and complexity and the accuracy required to describe small but functionally important changes in their properties. In this paper, we explore an alternative approach to calculating the excited states of large biochromophores, exemplified by a specific method for calculating the Qy transition of bacteriochlorophyll a, which we name Chl-xTB. Using a diagonally dominant approximation to the Casida equation and a bespoke parameterization scheme, Chl-xTB can match time-dependent density functional theory's accuracy and semiempirical speed for calculating the potential energy surfaces and absorption spectra of chlorophylls. We demonstrate that Chl-xTB (and other prospective realizations of our protocol) can be integrated into multiscale models, including concurrent excitonic and point-charge embedding frameworks, enabling the analysis of biochromophore networks in a native environment. We exploit this capability to probe the low-frequency spectral densities of excitonic energies and interchromophore interactions in the light harvesting antenna protein LH2 (light harvesting complex 2). The impact of low-frequency protein motion on interchromophore coupling and exciton transport has routinely been ignored due to the prohibitive costs of including it in simulations. Our results provide a more rigorous basis for continued use of this approximation by demonstrating that exciton transition energies are unaffected by low-frequency vibrational coupling to exciton interaction energies. [ABSTRACT FROM AUTHOR]

This paper presents a state-averaged complete active space self-consistent field (SA-CASSCF) in the atomic multipole optimized energetics for biomolecular application (AMOEBA) polarizable water model, which enables rigorous simulation of non-adiabatic molecular dynamics with nonequilibrium solvation effects. The molecular orbital and configuration interaction coefficients of the solute wavefunction, and the induced dipoles on solvent atoms, are solved by minimizing the state averaged energy variationally. In particular, by formulating AMOEBA water models and the polarizable continuum model (PCM) in a unified way, the algorithms developed for computing SA-CASSCF/PCM energies, analytical gradients, and non-adiabatic couplings in our previous work can be generalized to SA-CASSCF/AMOEBA by properly substituting a specific list of variables. Implementation of this method will be discussed with the emphasis on how the calculations of different terms are partitioned between the quantum chemistry and molecular mechanics codes. We will present and discuss results that demonstrate the accuracy and performance of the implementation. Next, we will discuss results that compare three solvent models that work with SA-CASSCF, i.e., PCM, fixed-charge force fields, and the newly implemented AMOEBA. Finally, the new SA-CASSCF/AMOEBA method has been interfaced with the ab initio multiple spawning method to carry out non-adiabatic molecular dynamics simulations. This method is demonstrated by simulating the photodynamics of the model retinal protonated Schiff base molecule in water. [ABSTRACT FROM AUTHOR]

Lyotropic liquid crystal phases (LCPs) are widely studied for diverse applications, including protein crystallization and drug delivery. The structure and properties of LCPs vary widely depending on the composition, concentration, temperature, pH, and pressure. High-throughput structural characterization approaches, such as small-angle x-ray scattering (SAXS), are important to cover meaningfully large compositional spaces. However, high-throughput LCP phase analysis for SAXS data is currently lacking, particularly for patterns of multiphase mixtures. In this paper, we develop semi-automated software for high throughput LCP phase identification from SAXS data. We validate the accuracy and time-savings of this software on a total of 668 SAXS patterns for the LCPs of the amphiphile hexadecyltrimethylammonium bromide (CTAB) in 53 acidic or basic ionic liquid derived solvents, within a temperature range of 25–75 °C. The solvents were derived from stoichiometric ethylammonium nitrate (EAN) or ethanolammonium nitrate (EtAN) by adding water to vary the ionicity, and adding precursor ions of ethylamine, ethanolamine, and nitric acid to vary the pH. The thermal stability ranges and lattice parameters for CTAB-based LCPs obtained from the semi-automated analysis showed equivalent accuracy to manual analysis, the results of which were previously published. A time comparison of 40 CTAB systems demonstrated that the automated phase identification procedure was more than 20 times faster than manual analysis. Moreover, the high throughput identification procedure was also applied to 300 unpublished scattering patterns of sodium dodecyl-sulfate in the same EAN and EtAN based solvents in this study, to construct phase diagrams that exhibit phase transitions from micellar, to hexagonal, cubic, and lamellar LCPs. The accuracy and significantly low analysis time of the high throughput identification procedure validates a new, rapid, unrestricted analytical method for the determination of LCPs. [ABSTRACT FROM AUTHOR]

This paper presents a theoretical investigation of the design of a new actuator type made of anisotropic colloidal particles grafted with stimuli-responsive polymer chains. These artificial muscles combine the osmotic actuation principle of stimuli-responsive hydrogels with the structural alignment of colloidal liquid crystals to achieve directional motion. The solubility of the stimuli-responsive polymer in the neutral state, its degree of polymerization, the salt concentration, and the grafting density of the polymer chains on the surface of the colloidal particles are investigated and identified as important for actuator performance and tunability. The computational results suggest that the proposed mechanically active material matches or exceeds the performances of natural muscles and provide the guidelines for the realization of artificial muscles with predetermined actuation properties. [ABSTRACT FROM AUTHOR]

Light-induced orientation of gas phase molecules is a long-pursued goal in physics and chemistry. Here, we experimentally demonstrate a six-fold increase in the terahertz-induced orientation of iodomethane (CH3I) molecules at room temperature, provided by rotational pre-excitation with a moderately intense near-IR pulse. The paper highlights the underlying interference of multiple coherent transition pathways within the rotational coherence manifold and is analyzed accordingly. Our experimental and theoretical results provide desirable and practical means for all-optical experiments on oriented molecular ensembles. [ABSTRACT FROM AUTHOR]

Ambient pressure x-ray photoelectron spectroscopy (APXPS) can provide a compelling platform for studying an analyte's oxidation and reduction reactions in solutions. This paper presents proof-of-principle operando measurements of a model organometallic complex, iron hexacyanide, in an aqueous solution using the dip-and-pull technique. The data demonstrates that the electrochemically active liquid meniscuses on the working electrodes can undergo controlled redox reactions which were observed using APXPS. A detailed discussion of several critical experimental considerations is included as guidance for anyone undertaking comparable experiments. [ABSTRACT FROM AUTHOR]

Quantum mechanics has a deterministic Schrödinger equation for the wave function. The Göttingen–Copenhagen statistical interpretation is based on the Born Rule that interprets the wave function as a "probability amplitude." A precept of this interpretation is the lack of determinism in quantum mechanics. The Bohm interpretation is that the wave function is a source of a field experienced by the electrons, thereby attributing determinism to quantum theory. In this paper, we present a new perspective on such determinism. The ideas are based on the equations of motion or "Quantal Newtonian" Laws obeyed by each electron. These Laws, derived from the temporal and stationary-state Schrödinger equation, are interpreted in terms of "classical" fields whose sources are quantal expectations of Hermitian operators taken with respect to the wave function. According to the Second Law, each electron experiences an external field—the quantal Coulomb-Lorentz law. It also experiences an internal field representative of properties of the system: correlations due to Coulomb repulsion and Pauli principle; the density; kinetic effects; and an internal magnetic field component. There is a response field. The First Law states that the sum of the external and internal fields experienced by each electron vanishes. These fields are akin to those of classical physics: They pervade all space; their structure is descriptive of the quantum system; the energy of the system is stored in these fields. It is in the classical behavior of these fields, which arise from quantal sources that one may then speak of determinism in quantum mechanics. [ABSTRACT FROM AUTHOR]

Spiky coatings (also known as nanoforests or Fakir-like surfaces) have found many applications in chemical physics, material sciences, and biotechnology, such as superhydrophobic materials, filtration and sensing systems, and selective protein separation, to name but a few. In this paper, we provide a systematic study of steady-state diffusion toward a periodic array of absorbing cylindrical pillars protruding from a flat base. We approximate a periodic cell of this system by a circular tube containing a single pillar, derive an exact solution of the underlying Laplace equation, and deduce a simple yet exact representation for the total flux of particles onto the pillar. The dependence of this flux on the geometric parameters of the model is thoroughly analyzed. In particular, we investigate several asymptotic regimes, such as a thin pillar limit, a disk-like pillar, and an infinitely long pillar. Our study sheds light onto the trapping efficiency of spiky coatings and reveals the roles of pillar anisotropy and diffusional screening. [ABSTRACT FROM AUTHOR]

In this study, based on the theory developed in Paper I, we explore the combined effects of molecular fluorescence and excitation energy transfer in a minimal model—a pair of single-vibration-mode chromophores coupled to surface plasmon polaritons. For the chromophores with zero Huang–Rhys factors and strong couplings to surface plasmon polaritons, we find that the frequencies of Rabi oscillations (the strengths of strong light–matter couplings) are associated with the initial excitation conditions. On the other hand, for the chromophores weakly coupled to surface plasmon polaritons, our numerical calculations together with analytical analysis elaborate on the conditions for the superradiant and subradiant decay behaviors. Moreover, we show that the modified decay rate constants can be explicitly expressed in terms of generalized spectral densities (or dyadic Green's functions), revealing a relationship between photonic environments and the collective effects such as superradiance and subradiance. For the chromophores with nonzero Huang–Rhys factors and strong coupling to surface plasmon polaritons, the effects of molecular vibrations emerge. We demonstrate that the low-frequency vibrational modes do not affect the excited state population dynamics, while the high-frequency vibrational modes can modify either the period of Rabi oscillation (Franck–Condon Rabi oscillation) or the amplitude of excited state population. Our study shows that the collective effects, including superradiance and subradiance, can be controlled via dielectric environments and initial excitation conditions, providing new insights into polariton chemistry and the design of quantum optical devices. [ABSTRACT FROM AUTHOR]

Some diarylethene molecular switches have a low quantum yield for cycloreversion when excited by a single photon, but react more efficiently following sequential two-photon excitation. The increase in reaction efficiency depends on both the relative time delay and the wavelength of the second photon. This paper examines the wavelength-dependent mechanism for sequential excitation using excited-state resonance Raman spectroscopy to probe the ultrafast (sub-30 fs) dynamics on the upper electronic state following secondary excitation. The approach uses femtosecond stimulated Raman scattering (FSRS) to measure the time-gated, excited-state resonance Raman spectrum in resonance with two different excited-state absorption bands. The relative intensities of the Raman bands reveal the initial dynamics in the higher-lying states, Sn, by providing information on the relative gradients of the potential energy surfaces that are accessed via secondary excitation. The excited-state resonance Raman spectra reveal specific modes that become enhanced depending on the Raman excitation wavelength, 750 or 400 nm. Many of the modes that become enhanced in the 750 nm FSRS spectrum are assigned as vibrational motions localized on the central cyclohexadiene ring. Many of the modes that become enhanced in the 400 nm FSRS spectrum are assigned as motions along the conjugated backbone and peripheral phenyl rings. These observations are consistent with earlier measurements that showed higher efficiency following secondary excitation into the lower excited-state absorption band and illustrate a powerful new way to probe the ultrafast dynamics of higher-lying excited states immediately following sequential two-photon excitation. [ABSTRACT FROM AUTHOR]

The formalism of the generalized quantum master equation (GQME) is an effective tool to simultaneously increase the accuracy and the efficiency of quasiclassical trajectory methods in the simulation of nonadiabatic quantum dynamics. The GQME expresses correlation functions in terms of a non-Markovian equation of motion, involving memory kernels that are typically fast-decaying and can therefore be computed by short-time quasiclassical trajectories. In this paper, we study the approximate solution of the GQME, obtained by calculating the kernels with two methods: Ehrenfest mean-field theory and spin-mapping. We test the approaches on a range of spin–boson models with increasing energy bias between the two electronic levels and place a particular focus on the long-time limits of the populations. We find that the accuracy of the predictions of the GQME depends strongly on the specific technique used to calculate the kernels. In particular, spin-mapping outperforms Ehrenfest for all the systems studied. The problem of unphysical negative electronic populations affecting spin-mapping is resolved by coupling the method with the master equation. Conversely, Ehrenfest in conjunction with the GQME can predict negative populations, despite the fact that the populations calculated from direct dynamics are positive definite. [ABSTRACT FROM AUTHOR]

In a recent paper [A. D. Becke, J. Chem. Phys. 156, 214101 (2022)], we compared two Kohn–Sham density functionals based on physical modeling and theory with the best density-functional power series fits in the literature. With only a handful of physically motivated pre-factors, our functionals matched, and even slightly exceeded, the performance of the best power-series functionals on the general main group thermochemistry, kinetics, and noncovalent interactions (GMTKN55) chemical database of Goerigk et al. [Phys. Chem. Chem. Phys. 19, 32184 (2017)]. This begs the question: how much can their performance be improved by adding power-series terms of our own? We address this question in the present work. First, we describe a series expansion variable that we believe contains more local physics than any other variable considered to date. Then we undertake modest, one-dimensional fits to the GMTKN55 data with our theory-based functional corrected by power-series exchange and dynamical correlation terms. We settle on 12 power-series terms (plus six parent terms) and achieve the lowest GMTKN55 "WTMAD2" error yet reported, by a substantial margin, for a hybrid Kohn–Sham density functional. The new functional is called "B22plus." [ABSTRACT FROM AUTHOR]

Recent experiments of chemical reactions in optical cavities have shown great promise to alter and steer chemical reactions, but still remain poorly understood theoretically. In particular, the origin of resonant effects between the cavity and certain vibrational modes in the collective limit is still subject to active research. In this paper, we study the unimolecular dissociation reactions of many molecules, collectively interacting with an infrared cavity mode, through their vibrational dipole moment. We find that the reaction rate can slow down by increasing the number of aligned molecules, if the cavity mode is resonant with a vibrational mode of the molecules. We also discover a simple scaling relation that scales with the collective Rabi splitting, to estimate the onset of reaction rate modification by collective vibrational strong coupling and numerically demonstrate these effects for up to 104 molecules. [ABSTRACT FROM AUTHOR]

In this paper, we report an open system consisting of three self-replicating peptides, in which peptide 1 inhibits the duplex template of peptide 2, peptide 2 inhibits duplex 3, and peptide 3 inhibits duplex 1 to complete the negative feedback loop. This interacting chemical network yields oscillations in the concentrations of all species over time and establishes a possible mechanism for pre-biotic chemical systems organization. The first focus of our analysis is the effect of altering rates of duplex formation and inhibition on oscillations. We then examine the autocatalytic rate constant in the symmetric and asymmetric cases. [ABSTRACT FROM AUTHOR]

This paper provides results of dielectric relaxation (DR) spectroscopy of aqueous solutions of tris(hydroxymethyl)aminomethane (Tris) covering frequencies of 0.05 ≤ ν/GHz ≤ 89. The DR spectra can be well fit by a sum of Cole–Cole relaxation, assigned to the solute, and 2 Debye modes already observed for neat water. Analysis of the amplitudes reveals that Tris is hydrated by 7 H2Os up to its solubility limit. However, the rather high effective solute dipole moment of ≈12 D suggests that H2O dipoles in contact with Tris should reorient independently from it. Accordingly, an alternative description of the DR spectra with a superposition of 4 Debyerelaxations was attempted. In this model, the slowest mode at ∼4 GHz arises from solute reorientation and that at ∼8 GHz was assigned to dynamically retarded hydration water, whereas relaxations at ∼18 and ∼500 GHz are again those of (rather unperturbed) bulk water. Analysis of the solvent-related modes shows that Tris indeed slows down 7–8 H2O molecules. However, the solute–solvent interaction strength is rather weak, excluding the rotation of an alleged Tris-(7–8) H2O cluster as an entity. The now derived effective dipole moment of (6.3 ± 0.5) D for the bare Tris molecule allows speculations on its conformation. With the help of computational methods, we suggest that Tris dissolved in water most likely possesses an intramolecular H-bond between the nitrogen and hydrogen atoms of amino and hydroxyl groups, respectively. In addition, computational results indicate that the seven hydration H2Os found by DR bind directly to the Tris OH groups. [ABSTRACT FROM AUTHOR]

Dilute magnetic semiconductors (DMSs) with both charge and spin degrees of freedom have emerged as promising candidates in the spintronic industry. However, the Curie temperature below room temperature and uncertainty about the origin of ferromagnetism hinder the application of DMSs. To address these issues, we explored a better SnO2-based co-doped method (Fe, V + VSn) using ab initio calculations. The calculation results show that the Sn13FeVO32 (Fe, V + VSn) has a high Curie temperature (716 K), good ferromagnetic properties, stronger covalency of bonds, and better optical transparency in the visible light range. In addition, the holes or electrons generated by the complexes in the (Fe, V + VO/VSn) co-doped system cause a spin-polarized double exchange effect in the Fe-3d, V-3d, and O-2p orbitals, which leads to magnetism of the co-doped systems. The static dielectric constant ɛ1(0) of the system increases after doping. Among them, Sn14FeVO31 (Fe, V + VO) has the largest ɛ1(0), indicating that Sn14FeVO31 has the strongest polarization ability and better photocatalytic properties. In Sn14FeVO31, the imaginary part of the dielectric function and the absorption spectrum all have new peaks in the low-energy region, which are caused by the jump of electrons from the guide band of the spin-polarized impurity energy level. This paper proposes a new method for preparing dilute magnetic semiconductors in spin electronic devices with high room temperature ferromagnetic properties and excellent optical properties through the (Fe, V + VO/VSn) co-doped SnO2. [ABSTRACT FROM AUTHOR]

Seniority-zero geminal wavefunctions are known to capture bond-breaking correlation. Among this class of wavefunctions, Richardson–Gaudin states stand out as they are eigenvectors of a model Hamiltonian. This provides a clear physical picture, clean expressions for reduced density matrix (RDM) elements, and systematic improvement (with a complete set of eigenvectors). Known expressions for the RDM elements require the computation of rapidities, which are obtained by first solving for the so-called eigenvalue based variables (EBV) and then root-finding a Lagrange interpolation polynomial. In this paper, we obtain expressions for the RDM elements directly in terms of the EBV. The final expressions can be computed at the same cost as the rapidity expressions. Therefore, except, in particular, circumstances, it is entirely unnecessary to compute rapidities at all. The RDM elements require numerically inverting a matrix, and while this is usually undesirable, we demonstrate that it is stable, except when there is degeneracy in the single-particle energies. In such cases, a different construction would be required. [ABSTRACT FROM AUTHOR]

This paper presents further development of the new semi-classical trajectory-based formalism described in Paper I [Chistikov et al., J. Chem. Phys. 151, 194106 (2019)]. We report the results of simulation and analysis of the low-frequency collision-induced absorption (CIA) in CO2–Ar, including its true dimer component. Our consideration relies on the use of ab initio intermolecular potential energy and induced dipole surfaces for CO2–Ar calculated in an assumption of a rigid CO2 structure using the CCSD(T) method. The theory, the details of which are reported in Paper I [Chistikov et al., J. Chem. Phys. 151, 194106 (2019)], permits taking into account the effect of unbound and quasi-bound classical trajectories on the CIA in the range of a rototranslational band. This theory is largely extended by trajectory-based simulation of the true bound dimer absorption in the present paper. The spectra are obtained from a statistical average over a vast ensemble of classical trajectories restricted by properly chosen domains in the phase space. Rigorous classical theory is developed for two low-order spectral moments interpreted as the Boltzmann-weighted average of the respective dipole functions. These spectral moments were then used to check the accuracy of our trajectory-based spectra, for which both spectral moments can be evaluated independently in terms of specific integrals over the trajectory-based calculated spectral profiles. Good agreement between the spectral moments calculated as integrals over the frequency domain or the phase space largely supports the reliability of our simulated CIA spectra, which conform with the available microwave and far-infrared observations. [ABSTRACT FROM AUTHOR]

This Editorial reports how the depletion force theory was originally developed by Sho Asakura and Fumio Oosawa and how their one-page paper was "rediscovered" about 20 years after the paper was published. The first part of this Editorial is mostly based on the lecture by Oosawa and his autobiographies, and the second part is written by one of two scientists who found the paper. The aim of this Editorial is to record the background of the discovery of the depletion force. We believe that this Editorial presents an interesting story showing how science develops. The story reminds us of the importance of basic education and continuous interests in unknown phenomena and interactions between people of different disciplines, although they are sometimes considered as separate elements of research. [ABSTRACT FROM AUTHOR]

One may obviously think that the best way to control magnetic properties relies on using a magnetic field. However, it is not convenient to focus a magnetic field on a small object, whereas it is much easier to do so with an electric field. Magnetoelectric coupling allows one to control the magnetization with the electric field and the polarization with the magnetic field and could therefore provide a solution to this problem. This paper aims at quantifying the impact of the electric field on both the isotropic magnetic exchange and the Dzyaloshinskii–Moriya interaction in the case of a binuclear system of S = 1/2 spins. This study follows previous studies that showed that very high Dzyaloshinskii–Moriya interaction, i.e., the antisymmetric exchange, can be generated when close to first order spin orbit coupling. We will, therefore, explore this regime in a model Cu(II) complex that exhibits a quasi-degeneracy of the d x 2 − y 2 and dxy orbitals. This situation is indeed the one that allows us to obtain the largest spin orbit couplings in transition metal complexes. We will show that both the magnetic exchange and the Dzyaloshinskii–Moriya interaction are very sensitive to the electric field and that it would therefore be possible to modulate and control magnetic properties by the electric field. Finally, rationalizations of the obtained results will be proposed. [ABSTRACT FROM AUTHOR]

Nitroaromatic compounds can photorelease nitric oxide after UV absorption. The efficiency of the photoreaction depends on the molecular structure, and two features have been pointed out as particularly important for the yield of the process: the presence of methyl groups at the ortho position with respect to the nitro group and the degree of conjugation of the molecule. In this paper, we provide a theoretical characterization at the CASPT2//CASSCF (complete active space second-order perturbation theory//complete active space self-consistent field) level of theory of the photorelease of NO for four molecules derived from nitrobenzene through the addition of ortho methyl groups and/or the elongation of the conjugation. Our previously described mechanism obtained for the photorelease of NO in nitrobenzene has been adopted as a model for the process. According to this model, the process proceeds through a reactive singlet–triplet crossing (STC) region that the system can reach from the triplet 3(πOπ*) minimum. The energy barrier that must be surmounted in order to populate the reactive STC can be associated with the efficiency of the photoreaction. Here, the obtained results display clear differences in the efficiency of the photoreaction in the studied systems and can be correlated with experimental results. Thus, the model proves its ability to highlight the differences in the photoreaction efficiency for the nitroaromatic compounds studied here. [ABSTRACT FROM AUTHOR]

We have studied the persistence probability p(t) of an active Brownian particle with shape asymmetry in two dimensions. The persistence probability is defined as the probability of a stochastic variable that has not changed its sign in a given fixed time interval. We have investigated two cases: (1) diffusion of a free active particle and (2) that of a harmonically trapped particle. In our earlier work, by Ghosh et al. [J. Chem. Phys. 152, 174901 (2020)], we had shown that p(t) can be used to determine the translational and rotational diffusion constant of an asymmetrically shaped particle. The method has the advantage that the measurement of the rotational motion of the anisotropic particle is not required. In this paper, we extend the study to an active anisotropic particle and show how the persistence probability of an anisotropic particle is modified in the presence of a propulsion velocity. Furthermore, we validate our analytical expression against the measured persistence probability from the numerical simulations of single particle Langevin dynamics and test whether the method proposed in our earlier work can help distinguish between active and passive anisotropic particles. [ABSTRACT FROM AUTHOR]

Adding carbon nanoparticles into organic phase change materials (PCMs) such as paraffin is a common way to enhance their thermal conductivity and to improve the efficiency of heat storage devices. However, the sedimentation stability of such blends can be low due to aggregation of aromatic carbon nanoparticles in the aliphatic paraffin environment. In this paper, we explore whether this important issue can be resolved by the introduction of a polymer agent such as poly(3-hexylthiophene) (P3HT) into the paraffin–nanoparticle blends: P3HT could ensure the compatibility of aromatic carbon nanoparticles with aliphatic paraffin chains. We employed a combination of experimental and computational approaches to determine the impact of P3HT addition on the properties of organic PCMs composed of paraffin and carbon nanoparticles (asphaltenes). Our findings clearly show an increase in the sedimentation stability of paraffin–asphaltene blends, when P3HT is added, through a decrease in average size of asphaltene aggregates as well as in an increase of the blends' viscosity. We also witness the appearance of the yield strength and gel-like behavior of the mixtures. At the same time, the presence of P3HT in the blends has almost no effect on their thermophysical properties. This implies that all properties of the blends, which are critical for heat storage applications, are well preserved. Thus, we demonstrated that adding polyalkylthiophenes to paraffin–asphaltene mixtures led to significant improvement in the performance characteristics of these systems. Therefore, the polymer additives can serve as promising compatibilizers for organic PCMs composed of paraffins and asphaltenes and other types of carbon nanoparticles. [ABSTRACT FROM AUTHOR]

Black-phase formamidinium lead iodide (FAPbI3), with a narrow bandgap and high thermal stability, has emerged as an in-demand material for highly efficient perovskite solar cells (PSCs). In a two-step sequential deposition, the PbI2 film plays an important role in the formation of a perovskite film with desirable qualities. This paper explores using N-methyl-2-pyrrolidone (NMP), a strong Lewis base, and N,N-dimethylformamide (DMF) as a mixed precursor solvent (DMF/NMP) of PbI2 and reports on preparing PbI2 films with a porous morphology by thermal treatment. Porous PbI2 films ensure the diffusion and sufficient reaction of the formamidinium iodide solution to form a smooth perovskite film. In addition, a dynamic spin coating method is also introduced to improve the uniformity of the perovskite film. Both methods yield a pure α-phase FAPbI3 film immediately in the unannealed state, which is necessary for the perovskite film to maintain phase stability. Finally, PSCs with a power conversion efficiency of 21.20% (0.13 cm2) are fabricated and optimized. The unencapsulated PSCs retain 90% of the initial efficiency for 1000 hours in dry air and exhibit a good thermal stability when heated to 85 °C. [ABSTRACT FROM AUTHOR]

More efficient excitation energy transfer is simulated for distances smaller than 5 nm, while at distances less than 1.5 nm, the opposite occurs, indicating an interplay between the excitation of the first nanoparticle and the back transfer from the second nanoparticle. In another article in this Special Issue, Herrera and Litinskaya[119] develop a macroscopic model of exciton-plasmon metamaterials composed of nanoparticle dimers with organic molecules (chromophores are considered in the paper) situated inside the gap. 155, 134117 (2021).10.1063/5.0060171 97 S. Wang, G. D. Scholes, and L.-Y. Hsu, "Quantum dynamics of a molecular emitter strongly coupled with surface plasmon polaritons: A macroscopic quantum electrodynamics approach", J. Chem. Phys. The field of plasmonics has been in the spotlight in physics,[1] materials science,[2] chemistry,[3] and biology[4] for more than a decade, with relevant applications in photovoltaics,[5] biosensing,[6] and nanoimaging.[7] Thanks to fascinating progress in nanofabrication, we now have tools to create artificial metal/dielectric interfaces with 1 nm precision,[[8]] while ultrafast optical detection techniques are now offering riveting time resolution as short as 2 fs.[10] Combined together, these tools are used to further advance our understanding of light-matter interaction at the nanoscale when the classical view collides and merges with quantum properties of matter and light.[11] The ability of plasmonic materials to sustain surface plasmon-polariton (SPP)[12] resonances results in a very small SPP mode volume. [Extracted from the article]

It is uncontested that a machine learning scheme cannot correctly reproduce physical properties that vary on a manifold in configuration space if the fingerprint, used as an input for the machine learning scheme, is constant on this manifold. In our original paper, we found several manifolds whose fingerprint variation is sufficiently small to prevent machine learning based on a standard training scheme even if some 40 configurations on the manifold are included in the training set. Response to "Comment on 'Manifolds of quasi-constant SOAP and ACSF fingerprints and the resulting failure to machine learn four-body interactions'" [J. Chem. Phys. [Extracted from the article]

In this paper, we have developed a unitary variant of a double exponential coupled cluster theory, which is capable of handling molecular strong correlation with arbitrary electronic complexity. With the Hartree–Fock determinant taken as the reference, we introduce a sequential product of parameterized unitary Ansätze. While the first unitary, containing the excitation operators, acts directly on the reference determinant, the second unitary, containing a set of rank-two, vacuum-annihilating scattering operators, has nontrivial action only on certain entangled states. We demonstrate the theoretical bottleneck of such an implementation in a classical computer, whereas the same is implemented in the hybrid quantum–classical variational quantum eigensolver framework with a reasonably shallow quantum circuit without any additional approximation. We have further introduced a number of variants of the proposed Ansatz with different degrees of sophistication by judiciously approximating the scattering operators. With a number of applications on strongly correlated molecules, we have shown that all our schemes can perform uniformly well throughout the molecular potential energy surface without significant additional implementation cost over the conventional unitary coupled cluster approach with single and double excitations. [ABSTRACT FROM AUTHOR]

Though jamming transitions are long studied in condensed matter physics and granular systems, much less is known about active jamming (or unjamming), which commonly takes place in living materials. In this paper, we explore, by molecular dynamic simulations, the jamming–unjamming transition in a dense system of active semiflexible filaments. In particular, we characterize the distinct impact of polar vs nematic driving for different filament rigidities and at varying densities. Our results show that high densities of dynamic active filaments can be achieved by only changing the nature of the active force, nematic or polar. Interestingly, while polar driving is more effective at unjamming the system at high densities below confluency, we find that at even higher densities, nematic driving enhances unjamming compared to its polar counterpart. The effect of varying the rigidity of filaments is also significantly different in the two cases: While for nematic driving, lowering the bending rigidity unjams the system, we find an intriguing reentrant jamming–unjamming–jamming transition for polar driving as the filament rigidity is lowered. While the first transition (unjamming) is driven by softening due to reduced rigidity, the second transition (jamming) is a cooperative effect of ordering and coincides with the emergence of nematic order in the system. Together, through a generic model of self-propelled flexible filaments, our results demonstrate how tuning the nature of self-propulsion and flexibility can be employed by active materials to achieve high densities without getting jammed. [ABSTRACT FROM AUTHOR]

This paper presents the experimental and theoretical analyses of ArICl(IP,vIP,nIP) states' population and decay at energies lower than the ArICl(E,vE = 0,nE) dissociation limit (IP = E0+, D′2, β1), vIP = 0, 1, and nIP are the quantum numbers of the van der Waals (vdW) modes. We have measured the excitation spectra of the ArICl(E,vE = 0,1,nE → X,vX,nX) and ArICl(β,0,nβ → A and/or D ′ , v D ′ , n D ′ → A ′ luminescence as well as luminescence spectra themselves. To construct potential energy surfaces (PESs) for valence (A1, A′2) and ion-pair (E, β, and D′) electronic states of the complex, we utilized the intermolecular diatomic-in-molecule perturbation theory first order method. The experimental and calculated spectroscopic characteristics of the T-shaped ArICl valence and E, β states agree well. The ArICl(D′) state PES has no vdW levels in the T-shaped configuration, and collinear ArICl(D′) binding energy is larger than that of the T-shaped ArICl(β) state. We calculated vibrational state energies and the ArICl(IP → valence states) luminescence excitation spectra, as well as luminescence spectra themselves, by using the Heidelberg MCTDH code. The comparison of the experimental and calculated excitation spectra shows that the latter describe their principal features. The bound–bound ArICl(E,0,nE → X and β,0,nβ → A) parts of experimental luminescence spectra are described adequately by the calculated spectra, whereas bound-free ArICl(E,0,nE → X, D′, 0, nD′ → A′) parts are not described since the bound-free transitions occur in repulsive parts of the ArICl(X, A ′ PESs, which we cannot describe accurately. [ABSTRACT FROM AUTHOR]

In this paper, we present the theory and implementation of nuclear magnetic resonance shielding constants with gauge-including atomic orbitals for the hybrid multiconfigurational short-range density functional theory model. As a special case, this implementation also includes Hartree–Fock srDFT (HF-srDFT). Choosing a complete-active space (CAS) wave function as the multiconfigurational parameterization of the wave function, we investigate how well CAS-srDFT reproduces experimental trends of nuclear shielding constants compared to DFT and complete active space self-consistent field (CASSCF). Calculations on the nucleobases adenine and thymine show that CAS-srDFT performs on average the best of the tested methods, much better than CASSCF but only marginally better than HF-srDFT. The performance, compared to regular DFT, is similar when functionals containing exact exchange are used. We generally find that the inclusion of exact exchange is important for an accurate description of the shielding constants. In cases where no exact exchange is included, we observe that the HF- and CAS-srDFT often outperform regular DFT. For calculations on transition metal nuclei in organometallic compounds with significant static correlation, the CAS-srDFT method again outperforms CASSCF compared to experimental shielding constants, and the change from HF-srDFT is substantial. In conclusion, the static correlation posed by the metal complexes seems to be captured by CAS-srDFT, which is promising since this type of correlation is not well described by regular DFT. [ABSTRACT FROM AUTHOR]

Polymer crystallization is a long-standing interesting problem both in fundamental polymer physics and in polymer manufacturing. Fundamentally, the connectivity of the macromolecules provides a unique feature for the study of nucleation and growth of crystals in contrast to the crystallization of smaller molecules. In addition, understanding the crystallization in polymers is industrially important due to the necessity of its control to achieve mechanically durable plastic materials. Molecular dynamics simulations offer a suitable way of studying this phenomenon due to their capability to probe the small time and length scales that are characteristic of nucleation. In this paper, we use a long alkane chain model to study nucleation and the growth of polyethylene crystals both within bulk and nano-droplets whose diameters range from about 6 to 11 nm. It is found that the droplets approach being more spherical during the nucleation whereas they deviate from this shape during the growth regime. Strikingly, a mean first passage time analysis indicates that the nucleation rate per unit volume decreases as the droplet size is increased. Finally, visual inspection of the simulation snapshots suggests that nuclei majorly emerged from the surface of the droplets. [ABSTRACT FROM AUTHOR]

Ytterbium monohydroxide is a promising molecule in the search for new physics. It is well known that levels of opposite parity, separated by energy split, the so-called l-doublets, define the experimental electric field strength required for molecule polarization. In addition, in our previous paper [Phys. Rev. A 105, L050801 (2022)], we have shown that the value of l-doubling directly influences the sensitivity of linear triatomic molecules toward the P , T -odd effects. In our work [J. Chem. Phys. 155, 164301 (2021)], we have calculated the value of l-doubling for YbOH molecules with the approximation of fixed O–H bond length. Taking the importance of this property into account, in the present study, we consider the additional degree of freedom corresponding to ligand (OH) deformation. [ABSTRACT FROM AUTHOR]

Experimental and theoretical studies have highlighted protonated water clusters (PWCs) as important models of the excess proton in aqueous systems. A significant focus has been characterizing the spectral signatures associated with different excess proton solvation motifs. Accurate vibrational frequency calculations are crucial for connecting the measured spectra to the structure of PWCs. In this paper, we extend and characterize a coupled local mode (CLM) approach for calculating the infrared spectra of PWCs using the H + ( H 2 O) 4 cluster as a benchmark system. The CLM method is relatively low cost and incorporates the anharmonicity and coupling of OH vibrations. Here, we demonstrate the accuracy of the technique compared to experiments. We also illustrate the dependence of calculated spectral features on the underlying electronic structure theory and basis sets used in the local mode frequency and coupling calculations. [ABSTRACT FROM AUTHOR]

We adopt a broken-symmetry strategy for evaluating effective magnetic constants J within the fully self-consistent GW method. To understand the degree of spin contamination present in broken-symmetry periodic solutions, we propose several extensive quantities demonstrating that the unrestricted self-consistent GW preserves the broken-symmetry character of the unrestricted Hartree–Fock solutions. The extracted J are close to the ones obtained from multireference wave-function calculations. In this paper, we establish a robust computational procedure for finding magnetic coupling constants from self-consistent GW calculations and apply it to solid antiferromagnetic nickel and manganese oxides. [ABSTRACT FROM AUTHOR]

Under certain conditions, the dynamics of coarse-grained models of solvated proteins can be described using a Markov state model, which tracks the evolution of populations of configurations. The transition rates among states that appear in the Markov model can be determined by computing the relative entropy of states and their mean first passage times. In this paper, we present an adaptive method to evaluate the configurational entropy and the mean first passage times for linear chain models with discontinuous potentials. The approach is based on event-driven dynamical sampling in a massively parallel architecture. Using the fact that the transition rate matrix can be calculated for any choice of interaction energies at any temperature, it is demonstrated how each state's energy can be chosen such that the average time to transition between any two states is minimized. The methods are used to analyze the optimization of the folding process of two protein systems: the crambin protein and a model with frustration and misfolding. It is shown that the folding pathways for both systems are comprised of two regimes: first, the rapid establishment of local bonds, followed by the subsequent formation of more distant contacts. The state energies that lead to the most rapid folding encourage multiple pathways, and they either penalize folding pathways through kinetic traps by raising the energies of trapping states or establish an escape route from the trapping states by lowering free energy barriers to other states that rapidly reach the native state. [ABSTRACT FROM AUTHOR]

Accurate theoretical prediction of the band offsets at interfaces of semiconductor heterostructures can often be quite challenging. Although density functional theory has been reasonably successful to carry out such calculations, efficient, accurate semilocal functionals are desirable to reduce the computational cost. In general, the semilocal functionals based on the generalized gradient approximation (GGA) significantly underestimate the bulk bandgaps. This, in turn, results in inaccurate estimates of the band offsets at the heterointerfaces. In this paper, we investigate the performance of several advanced meta-GGA functionals in the computational prediction of band offsets at semiconductor heterojunctions. In particular, we investigate the performance of r2SCAN (two times revised strongly constrained and appropriately normed functional), rMGGAC (revised semilocal functional based on cuspless hydrogen model and Pauli kinetic energy density functional), mTASK (modified Aschebrock and Kümmel meta-GGA functional), and local modified Becke–Johnson exchange-correlation functionals. Our results strongly suggest that these meta-GGA functionals for supercell calculations perform quite well, especially, when compared to computationally more demanding GW calculations. We also present band offsets calculated using ionization potentials and electron affinities, as well as band alignment via the branch point energies. Overall, our study shows that the aforementioned meta-GGA functionals can be used within the density functional theory framework to estimate the band offsets in semiconductor heterostructures with predictive accuracy. [ABSTRACT FROM AUTHOR]

Computational analysis of difference spectra between two analogous systems is a challenging issue, as reliable estimation of a tiny difference spectrum requires an extraordinary precision of the two original spectra. We have developed an alternative new method to calculate the difference spectra under background-free conditions, which greatly improved the efficiency of computation. In this paper, we report further improvement by using efficient parallel implementation and the time correlation formula based on time derivative quantities. As a consequence, the present work achieved further remarkable acceleration in the calculations of difference infrared and Raman spectra in the order of magnitude and thereby allowed us by analyzing these difference spectra at a practical cost of computation. [ABSTRACT FROM AUTHOR]