Search

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]

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]

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]

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]

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]

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]

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]

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]

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]

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]

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 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]

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]

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]

An error in the theoretical treatment by Hadinger et al., using the experimental data on HCl of Coxon and Ogilvie, is pointed out. [ABSTRACT FROM AUTHOR]

Presents a response to the commentary on the article 'Further Contributions to the Energy Levels of a Perturbed Anharmonic Oscillator: Application to Adiabatic Corrections.'

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]

The ability of biomolecules to exert forces on their surroundings or resist compression from the environment is essential in a variety of biologically relevant contexts. For filaments in the low-temperature limit and under a constant compressive force, Euler buckling theory predicts a sudden transition from a compressed state to a bent state in these slender rods. In this paper, we use a mean-field theory to show that if a semiflexible chain is compressed at a finite temperature with a fixed end-to-end distance (permitting fluctuations in the compressive forces), it exhibits a continuous phase transition to a buckled state at a critical level of compression. We determine a quantitatively accurate prediction of the transverse position distribution function of the midpoint of the chain that indicates this transition. We find that the mean compressive forces are non-monotonic as the extension of the filament varies, consistent with the observation that strongly buckled filaments are less able to bear an external load. We also find that for the fixed extension (isometric) ensemble, the buckling transition does not coincide with the local minimum of the mean force (in contrast to Euler buckling). We also show that the theory is highly sensitive to fluctuations in length in two dimensions and the buckling transition can still be accurately recovered by accounting for those fluctuations. These predictions may be useful in understanding the behavior of filamentous biomolecules compressed by fluctuating forces, relevant in a variety of biological contexts. [ABSTRACT FROM AUTHOR]

Quantum master equations (QMEs) provide a general framework for describing electronic dynamics within a complex molecular system. Off-diagonal QMEs (OD-QMEs) correspond to a family of QMEs that describe the electronic dynamics in the interaction picture based on treating the off-diagonal coupling terms between electronic states as a small perturbation within the framework of second-order perturbation theory. The fact that OD-QMEs are given in terms of the interaction picture makes it non-trivial to obtain Schrödinger picture electronic coherences from them. A key experimental quantity that relies on the ability to obtain accurate Schrödinger picture electronic coherences is the absorption spectrum. In this paper, we propose using a recently introduced procedure for extracting Schrödinger picture electronic coherences from interaction picture inputs to calculate electronic absorption spectra from the electronic dynamics generated by OD-QMEs. The accuracy of the absorption spectra obtained this way is studied in the context of a biexciton benchmark model, by comparing spectra calculated based on time-local and time-nonlocal OD-QMEs to spectra calculated based on a Redfield-type QME and the non-perturbative and quantum-mechanically exact hierarchical equations of motion method. [ABSTRACT FROM AUTHOR]

In this paper, we study the nuclear gradients of heat bath configuration interaction self-consistent field (HCISCF) wave functions and use them to optimize molecular geometries for various molecules. We show that HCISCF nuclear gradients are fairly insensitive to the size of the "selected" variational space, which allows us to reduce the computational cost without introducing significant errors. The ability of the HCISCF to treat larger active spaces combined with the flexibility for users to control the computational cost makes the method very attractive for studying strongly correlated systems, which require a larger active space than possible with a complete active space self-consistent field. Finally, we study the realistic catalyst, Fe(PDI), and highlight some of the challenges this system poses for density functional theory (DFT). We demonstrate how HCISCF can clarify the energetic stability of geometries obtained from DFT when the results are strongly dependent on the functional. We also use the HCISCF gradients to optimize geometries for this species and study the adiabatic singlet–triplet gap. During geometry optimization, we find that multiple near-degenerate local minima exist on the triplet potential energy surface. [ABSTRACT FROM AUTHOR]

In Paper I [J. R. Elliott, A. J. Schultz, and D. A. Kofke, J. Chem. Phys. 143, 114110 (2015)] of this series, a methodology was presented for computing the coefficients of a power series of the Helmholtz energy in reciprocal temperature, β, through density series based on cluster integral expansions. Previously, power series in β were evaluated by thermodynamic perturbation theory (TPT) using molecular simulation of a reference fluid. The present methodology uses cluster integrals to evaluate coefficients of the density expansion at each individual order of temperature. While Paper I [J. R. Elliott, A. J. Schultz, and D. A. Kofke, J. Chem. Phys. 143, 114110 (2015)] developed this methodology for square well (SW) spheres, the present work extends the methodology to Lennard-Jones (LJ) spheres, where the reference fluid is the Weeks-Chandler-Andersen potential. Comparisons of TPT coefficients computed from cluster integrals to those from molecular simulation show good agreement through third order in β when coefficients are expressed with effective approximants. Notably, the agreement for LJ spheres is much better than for SW spheres although fewer coefficients of the density series (B2–B5) are available than for SW spheres (B2–B6). The coefficients for Bi(β) of the reference fluid are shown to follow a simple relationship to the virial coefficients of hard sphere fluids, corrected for the temperature dependency of the equivalent hard sphere diameter. This lays the foundation for a correlation of the second virial coefficient of LJ spheres B2(β) that extrapolates to infinite order in temperature. This correlation of B2(β) provides a basis for estimating the low density limit of TPT coefficients at all orders in temperature, facilitating a recursive extrapolation formula to estimate TPT coefficients of fourth order and higher over the entire density range. The applicability of the resulting equation of state is demonstrated by computing the thermodynamic properties for LJ spheres and comparing to standard simulation results. [ABSTRACT FROM AUTHOR]

The present paper is the second part of a series of papers aimed at assessing the accuracy of alchemical computational approaches based on nonequilibrium techniques for solvation free energy of organic molecules in the context of molecular dynamics simulations. In Paper I [Procacci, J. Chem. Phys. 151, 144113 (2019)], we dealt with bidirectional estimates of solvation free energies using nonequilibrium approaches. Here, we assess accuracy and precision of unidirectional estimates with the focus on the Gaussian and Jarzynski estimators. We present a very simple methodology to increase the statistics in the work distribution, hence boosting the accuracy and precision of the Jarzynski unidirectional estimates at no extra cost, exploiting the observed decorrelation between the random variables represented by the Lennard-Jones solute-solvent recoupling or decoupling work and by the electrostatic work due to the charging/discharging of the solute in the solvent. [ABSTRACT FROM AUTHOR]

In the context of molecular dynamics simulations, alchemical approaches based on nonequilibrium techniques are recently emerging as a powerful method for the computation of solvation free energy of druglike compounds. Here, we present a rigorous and extensive analysis of the accuracy and precision of the method as a function of the parameters qualifying the nonequilibrium alchemical protocol (e.g., number and length of the nonequilibrium trajectories and decoupling or recoupling alchemical schedule) on a selection of drug-size organic compounds characterized by a nontrivial conformational free energy landscape. The study is organized in two contributions. The first paper includes a detailed description of method and of the conformational behavior of molecular systems. Results are focused on the accuracy and precision bidirectional estimates of solvation free energy, notably those based on the so-called Bennett acceptance ratio. In the second paper, unidirectional estimates for solvation free energy are analyzed in depth. [ABSTRACT FROM AUTHOR]

In continuation with Paper I [S. R. Reddy et al., J. Chem. Phys. 151, 054303 (2019)], the vibronic structure and dynamics of the 1 Σ u + electronic state of C15, C17, C19, and C21 chains in the coupled manifold of 1 Σ u + –1Πg–1Πu– 1 Σ g + electronic states have been investigated in this paper. The model vibronic Hamiltonian developed through extensive ab initio quantum chemistry calculations in Paper I is employed, and first principles nuclear dynamics calculations are carried out to obtain the photoabsorption band of the 1 Σ u + electronic state. Both time-independent and time-dependent quantum mechanical calculations are carried out to precisely locate the vibrational levels, assign them with the progression of vibrational modes, and elucidate the impact of both Renner-Teller and pseudo-Renner-Teller couplings on them. The nonradiative decay of the 1 Σ u + electronic state is studied, and it is found that the decay rate is comparable with the prediction made for them to be qualified as a carrier of diffuse interstellar bands in the literature. The theoretical results are found to be in good accord with the available experimental results. [ABSTRACT FROM AUTHOR]

Exact nonadiabatic quantum evolution preserves many geometric properties of the molecular Hilbert space. In the first paper of this series ["Paper I," S. Choi and J. Vaníček, J. Chem. Phys. 150, 204112 (2019)], we presented numerical integrators of arbitrary-order of accuracy that preserve these geometric properties exactly even in the adiabatic representation, in which the molecular Hamiltonian is not separable into kinetic and potential terms. Here, we focus on the separable Hamiltonian in diabatic representation, where the split-operator algorithm provides a popular alternative because it is explicit and easy to implement, while preserving most geometric invariants. Whereas the standard version has only second-order accuracy, we implemented, in an automated fashion, its recursive symmetric compositions, using the same schemes as in Paper I, and obtained integrators of arbitrary even order that still preserve the geometric properties exactly. Because the automatically generated splitting coefficients are redundant, we reduce the computational cost by pruning these coefficients and lower memory requirements by identifying unique coefficients. The order of convergence and preservation of geometric properties are justified analytically and confirmed numerically on a one-dimensional two-surface model of NaI and a three-dimensional three-surface model of pyrazine. As for efficiency, we find that to reach a convergence error of 10−10, a 600-fold speedup in the case of NaI and a 900-fold speedup in the case of pyrazine are obtained with the higher-order compositions instead of the second-order split-operator algorithm. The pyrazine results suggest that the efficiency gain survives in higher dimensions. [ABSTRACT FROM AUTHOR]

Vapor-phase ammonia, NH3(g), and hydrochloric acid, HCl(g), undergo a series of complex reactions, including nucleation and growth, to form solid ammonium chloride, NH4Cl(s). The counterdiffusional experiment, whereby HCl(g) and NH3(g) diffuse from opposite ends of a tube and react to form spatiotemporally complex patterns, has a rich history of study. In this paper, we combine experimental data, molecular simulations, and analysis and simulations of a partial differential equation model to address the questions of where the first unobserved vapor product NH4Cl(g) and visually observable precipitate NH4Cl(s) form and how these positions depend on experimental parameters. These analyses yield a consistent picture which involves a moving reaction front as well as previously unobserved heterogeneous nucleation, wall nucleation, and homogeneous nucleation. The experiments combined with modeling allow for an estimate of the heterogeneous and homogeneous nucleation thresholds for the vapor-to-solid phase transition. The results, synthesized with the literature on this vapor-to-particle reaction, inform a discussion of the details of the reaction mechanism, including the role of water, which concludes the paper. [ABSTRACT FROM AUTHOR]

Cluster perturbation (CP) theory was developed in Paper I [F. Pawłowski et al., J. Chem. Phys. 150, 134108 (2019)] for a coupled cluster (CC) target state and is extended in this paper to comprehend a cluster linear (CL) target state, for which the embedding of a CC parent state in the target excitation space is described using a linear parametrization. The theory is developed for determining the energy and molecular properties for a CL state. When CP theory is applied to a CL target state, a series of corrections is determined in orders of the CC parent-state similarity-transformed fluctuation potential, where the zeroth-order term is the energy or molecular property of the CC parent state and where the series formally converges to the energy or molecular property of the CL target state. The determination of energies and molecular properties is simpler for a CL state than for a CC state because the CL state is linearly parametrized. The amplitude equations are quadratic for a CL target state, while quartic for a CC target state, and molecular property expressions for a CL target state have the same simple structure as for a configuration interaction state. The linear parametrization introduces non-size-extensive contributions in the energy and molecular property expressions. However, since the linear parametrization describes the embedding of the CC parent state in the target excitation space, the energy and molecular properties for a CL state are weakly size-extensive. For the energy, weak size-extensivity means that non-size-extensive contributions enter in sixth and higher orders in the CP energy series, whereas for molecular properties, weak size-extensivity means that non-size-extensive contributions enter in second and higher orders. Weak size-extensivity therefore has a little or vanishing effect on calculated energies or molecular properties. The determination of the CP energy and molecular property corrections does not require that amplitude or response equations are solved explicitly for the target state and it becomes computationally tractable to use low-order corrections from these series to obtain energies and molecular properties of CL target state quality. For three simple molecules, HF, N2, and CH2, the accuracy of the CL approach for ground-state energies is tested using a parent state including single and double excitations (i.e., the CC singles-and-doubles state, CCSD) and a target state that includes triple excitations. It is found that the size-extensive fifth-order CL energies deviate by less than 0.0001 hartree from the energies of a target CC that includes triple excitations (i.e., the CC singles-doubles-and-triples state, CCSDT). CP theory with a CL target state therefore becomes a very attractive replacement of standard CC theory for high-accuracy energy and molecular property calculations, in which triple and higher excitation levels are considered. [ABSTRACT FROM AUTHOR]