Search

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]

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]

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]

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]

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]

Often, a single correlation function is used to measure the kinetics of a complex system. In contrast, a large set of k-vector modes and their correlation functions are commonly defined for motion in free space. This set can be transformed to the van Hove correlation function, which is the Green's function for molecular diffusion. Here, these ideas are generalized to other observables. A set of correlation functions of nonlinear functions of an observable is used to extract the corresponding Green's function. Although this paper focuses on nonlinear correlation functions of an equilibrium time series, the results are directly connected to other types of nonlinear kinetics, including perturbation–response experiments with strong fields. Generalized modes are defined as the orthogonal polynomials associated with the equilibrium distribution. A matrix of mode-correlation functions can be transformed to the complete, single-time-interval (1D) Green's function. Diagonalizing this matrix finds the eigendecays. To understand the advantages and limitation of this approach, Green's functions are calculated for a number of models of complex dynamics within a Gaussian probability distribution. Examples of non-diffusive motion, rate heterogeneity, and range heterogeneity are examined. General arguments are made that a full set of nonlinear 1D measurements is necessary to extract all the information available in a time series. However, when a process is neither dynamically Gaussian nor Markovian, they are not sufficient. In those cases, additional multidimensional measurements are needed. [ABSTRACT FROM AUTHOR]

Allowing triplet components of individual geminals, spin-contaminated strongly orthogonal geminal wave functions may emerge, which can be ameliorated by spin-projection techniques. Of the latter, half-projection was previously shown to be useful, offering a compromise between the amount of remaining spin-contamination and the violation of size consistency generated by projection. This paper investigates how a half-projected spin-contaminated geminal wave function can be improved by multi-configuration perturbation theory to incorporate dynamical correlation effects. [ABSTRACT FROM AUTHOR]

The generalized quantum master equation (GQME) provides a general and formally exact framework for simulating the reduced dynamics of open quantum systems. The recently introduced modified approach to the GQME (M-GQME) corresponds to a specific implementation of the GQME that is geared toward simulating the dynamics of the electronic reduced density matrix in systems governed by an excitonic Hamiltonian. Such a Hamiltonian, which is often used for describing energy and charge transfer dynamics in complex molecular systems, is given in terms of diabatic electronic states that are coupled to each other and correspond to different nuclear Hamiltonians. Within the M-GQME approach, the effect of the nuclear degrees of freedom on the time evolution of the electronic density matrix is fully captured by a memory kernel superoperator, which can be obtained from short-lived (compared to the time scale of energy/charge transfer) projection-free inputs. In this paper, we test the ability of the M-GQME to predict the energy transfer dynamics within a seven-state benchmark model of the Fenna–Matthews–Olson (FMO) complex, with the short-lived projection-free inputs obtained via the Ehrenfest method. The M-GQME with Ehrenfest-based inputs is shown to yield accurate results across a wide parameter range. It is also found to dramatically outperform the direct application of the Ehrenfest method and to provide better-behaved convergence with respect to memory time in comparison to an alternative implementation of the GQME approach previously applied to the same FMO model. [ABSTRACT FROM AUTHOR]

The Infrared (IR) and Raman spectra of various defects in silicon, containing both oxygen atoms (in the interstitial position, Oi) and a vacancy, are computed at the quantum mechanical level by using a periodic supercell approach based on a hybrid functional (B3LYP), an all-electron Gaussian-type basis set, and the Crystal code. The first of these defects is VO: the oxygen atom, twofold coordinated, saturates the unpaired electrons of two of the four carbon atoms on first neighbors of the vacancy. The two remaining unpaired electrons on the first neighbors of the vacancy can combine to give a triplet (Sz = 1) or a singlet (Sz = 0) state; both states are investigated for the neutral form of the defect, together with the doublet solution, the ground state of the negatively charged defect. Defects containing two, three, and four oxygen atoms, in conjunction with the vacancy V, are also investigated as reported in many experimental papers: VO2 and VOOi (two oxygen atoms inside the vacancy, or one in the vacancy and one in interstitial position between two Si atoms) and VO2Oi and VO22Oi (containing three and four oxygen atoms). This study integrates and complements a recent investigation referring to Oi defects [Gentile et al., J. Chem. Phys. 152, 054502 (2020)]. A general good agreement is observed between the simulated IR spectra and experimental observations referring to VOx (x = 1–4) defects. [ABSTRACT FROM AUTHOR]

The optical surface of extreme ultraviolet (EUV) lithography machines is highly vulnerable to contamination by hydrocarbons, resulting in the formation of carbon deposits that significantly degrade the quality and efficiency of lithography. The dynamic gas lock (DGL) has been proven as an effective approach to alleviate carbon deposition. However, the majority of existing studies on carbon deposition neglect the influence of the DGL. This paper is dedicated to investigating the phenomena of hydrocarbon adsorption, desorption, and cleavage with considering the effects of the DGL. A comprehensive mathematical model of the carbon deposition process is established, and the impact of radiation intensity, temperature, and hydrocarbon types on the depositing rate is considered. The results suggest that the primary cause of carbon deposition is the direct cracking of hydrocarbons induced by photons with a wavelength range between 12.5 and 14.5 nm. Additionally, it has been observed that the carbon deposition rate decreases exponentially as clean gas flow increases when EUV radiation intensity exceeds 50 mW/mm2. Conversely, at low EUV radiation intensity, clean gas flow has little effect on the carbon deposition rate. An effective approach to mitigate carbon deposition is to elevate the temperature of the optical surface and employ light hydrocarbon materials in the EUV process. [ABSTRACT FROM AUTHOR]

Previous studies have shown that NaB6, KB6, and RbB6 adopting Pm 3 ̄ m are superconductors with a relatively high Tc under ambient conditions. In this paper, we conducted systematic structural and related properties research on CsB6 through a genetic evolution algorithm and total energy calculations based on density functional theory between 0 and 20 GPa. Our results reveal a cubic Pm 3 ̄ m CsB6, which is dynamically stable under the pressures we studied. We systematically calculated the formation enthalpies, electronic properties, and superconducting properties of Pm 3 ̄ m MB6 (M = Na, K, Rb, Cs). They all exhibit metallic features, and boron has high contributions to band structures, density of states, and electron–phonon coupling (EPC). The calculated results about the Helmholtz free energy difference of Pm 3 ̄ m CsB6 at 0, 10, and 20 GPa indicate that it is stable upon chemical decomposition (decomposition to simple substances Cs and B) from 0 to 400 K. The phonon density of states indicates that boron atoms occupy the high frequency area. The EPC results show that Pm 3 ̄ m CsB6 is a superconductor with Tc = 11.7 K at 0 GPa, close to NaB6 (13.1 K), KB6 (11.7 K), and RbB6 (11.3 K) at 0 GPa in our work, which indicates that boron atoms play an essential role in superconductivity: vibrations of B6 regular octagons lead to the high Tc of Pm 3 ̄ m MB6. Our work about Pm 3 ̄ m hexaborides provides a supplementary study on the borides of the group IA elements (without Fr and Li) and has an important guiding significance for the experimental synthesis of CsB6. [ABSTRACT FROM AUTHOR]

This paper presents a general second-quantized form of a permutation operator interchanging n pairs of electrons between interacting subsystems in the framework of the symmetry-adapted perturbation theory (SAPT). We detail the procedure for constructing this operator through the consecutive multiplication of single-pair permutation operators. This generalized form of the permutation operator has enabled the derivation of universal formulas for S2n approximations of the exchange energies in the first and second order of the interaction operator. We present expressions for corrections of S4 approximations and assess its efficacy on a selection of systems anticipated to exhibit a slowly converging overlap expansion. Additionally, we outline a method to sum the overlap expansion series to infinity in second-quantization, up to the second order in V. This new approach offers an alternative to the existing formalism based on density-matrix formulations. When combined with a symbolic algebra program for automated derivations, it paves the way for advancements in SAPT theory, particularly for intricate wavefunction theories. [ABSTRACT FROM AUTHOR]

Modern 4-wave mixing spectroscopies are expensive to obtain experimentally and computationally. In certain cases, the unfavorable scaling of quantum dynamics problems can be improved using a generalized quantum master equation (GQME) approach. However, the inclusion of multiple (light–matter) interactions complicates the equation of motion and leads to seemingly unavoidable cubic scaling in time. In this paper, we present a formulation that greatly simplifies and reduces the computational cost of previous work that extended the GQME framework to treat arbitrary numbers of quantum measurements. Specifically, we remove the time derivatives of quantum correlation functions from the modified Mori–Nakajima–Zwanzig framework by switching to a discrete-convolution implementation inspired by the transfer tensor approach. We then demonstrate the method's capabilities by simulating 2D electronic spectra for the excitation-energy-transfer dimer model. In our method, the resolution of data can be arbitrarily coarsened, especially along the t2 axis, which mirrors how the data are obtained experimentally. Even in a modest case, this demands O (1 0 3 ) fewer data points. We are further able to decompose the spectra into one-, two-, and three-time correlations, showing how and when the system enters a Markovian regime where further measurements are unnecessary to predict future spectra and the scaling becomes quadratic. This offers the ability to generate long-time spectra using only short-time data, enabling access to timescales previously beyond the reach of standard methodologies. [ABSTRACT FROM AUTHOR]

This paper reports the energies and charge and spin distributions of the low-lying excited states in singlet and triplet N2V defects in diamond from direct Δ-SCF calculations based on Gaussian orbitals within the B3LYP, PBE0, and HSE06 functionals. They assign the observed absorption at 2.463 eV, first reported by Davies et al. [Proc. R. Soc. London 351, 245 (1976)], to the excitation of a N(sp3) lone-pair electron in the singlet and triplet states, respectively, with estimates of ∼1.1 eV for that of the unpaired electrons, C(sp3). In both cases, the excited states are predicted to be highly local and strongly excitonic with 81% of the C(sp3) and 87% of the N(sp3) excited charges localized at the three C atoms nearest neighbor (nn) to the excitation sites. Also reported are the higher excited gap states of both the N lone pair and C unpaired electron. Calculated excitation energies of the bonding sp3 hybrids of the C atoms nn to the four inner atoms are close to that of the bulk, which indicates that the N2V defect is largely a local defect. The present results are in broad agreement with those reported by Udvarhelyi et al. [Phys. Rev. B 96, 155211 (2017)] from plane wave HSE06 calculations, notably for the N lone pair excitation energy, for which both predict an energy of ∼2.7 eV but with a difference of ∼0.5 eV for the excitation of the unpaired electron. [ABSTRACT FROM AUTHOR]

Many methods to accelerate sampling of molecular configurations are based on the idea that temperature can be used to accelerate rare transitions. These methods typically compute equilibrium properties at a target temperature using reweighting or through Monte Carlo exchanges between replicas at higher temperatures. A recent paper [G. M. Rotskoff and E. Vanden-Eijnden, Phys. Rev. Lett. 122, 150602 (2019)] demonstrated that accurate equilibrium densities of states can also be computed through a nonequilibrium "quench" process, where sampling is performed at a higher temperature to encourage rapid mixing and then quenched to lower energy states with dissipative dynamics. Here, we provide an implementation of the quench dynamics in LAMMPS and evaluate a new formulation of nonequilibrium estimators for the computation of partition functions or free energy surfaces (FESs) of molecular systems. We show that the method is exact for a minimal model of N-independent harmonic springs and use these analytical results to develop heuristics for the amount of quenching required to obtain accurate sampling. We then test the quench approach on alanine dipeptide, where we show that it gives an FES that is accurate near the most stable configurations using the quench approach but disagrees with a reference umbrella sampling calculation in high FE regions. We then show that combining quenching with umbrella sampling allows the efficient calculation of the free energy in all regions. Moreover, by using this combined scheme, we obtain the FES across a range of temperatures at no additional cost, making it much more efficient than standard umbrella sampling if this information is required. Finally, we discuss how this approach can be extended to solute tempering and demonstrate that it is highly accurate for the case of solvated alanine dipeptide without any additional modifications. [ABSTRACT FROM AUTHOR]

The network structure within polymers significantly influences their mechanical properties, including their strength, toughness, and fatigue resistance. All-atom molecular dynamics (AAMD) simulations offer a method to investigate the energy dissipation mechanism within polymers during deformation and fracture; Such an approach is, however, computationally inefficient when used to analyze polymers with complex network structures, such as the common chemically double-networked hydrogels. Alternatively, coarse-grained molecular dynamics (CGMD) models, which reduce the computational degrees of freedom by concentrating a set of adjacent atoms into a coarse-grained bead, can be employed. In CGMD simulations, a coarse-grained force field (CGFF) is a critical factor affecting the simulation accuracy. In this paper, we proposed a data-based method for predicting the CGFF parameters to improve the simulation efficiency of complex cross-linked network in polymers. Here, we utilized a typical chemically double-networked hydrogel as an example. An artificial neural network was selected, and it was trained with the tensile stress–strain data from the CGMD simulations using different CGFF parameters. The CGMD simulations using the predicted CGFF parameters show good agreement with the AAMD simulations and are almost fifty times faster. The data-inspired CGMD model presented here broadens the applicability of molecular dynamics simulations to cross-linked polymers and has the potential to provide insights that will aid the design of polymers with desirable mechanical properties. [ABSTRACT FROM AUTHOR]

Stimulated Raman Adiabatic Passage (STIRAP) and its variants, such as M-type chainwise-STIRAP, allow for efficiently transferring the populations in a multilevel system and have widely been used to prepare molecules in their rovibrational ground state. However, their transfer efficiencies are generally imperfect. The main obstacle is the presence of losses and the requirement to make the dynamics adiabatic. To this end, in the present paper, a new theoretical method is proposed for the efficient and robust creation and detection of deeply bound molecules in three-level Λ-type and five-level M-type systems via "Invariant-based shortcut-to-adiabaticity." In the regime of large detunings, we first reduce the dynamics of three- and five-level molecular systems to those of effective two- and three-level counterparts. By doing so, the major molecular losses from the excited states can be well suppressed. Consequently, the effective two-level counterpart can be directly compatible with two different "Invariant-based Inverse Engineering" protocols; the results show that both protocols give a comparable performance and have a good experimental feasibility. For the effective three-level counterpart, by considering a relation among the four incident pulses, we show that this model can be further generalized to an effective Λ-type one with the simplest resonant coupling. This generalized model permits us to borrow the "Invariant-based Inverse Engineering" protocol from a standard three-level Λ-type system to a five-level M-type system. Numerical calculations show that the weakly bound molecules can be efficiently transferred to their deeply bound states without strong laser pulses, and the stability against parameter variations is well preserved. Finally, the detection of ultracold deeply bound molecules is discussed. [ABSTRACT FROM AUTHOR]

Recent experiments by Brückner et al. [Science 380, 1357 (2023)] have observed an anomalous chain length dependence of the time of near approach of widely separated pairs of genomic elements on transcriptionally active chromosomal DNA. In this paper, I suggest that the anomaly may have its roots in internal friction between neighboring segments on the DNA backbone. The basis for this proposal is a model of chain dynamics formulated in terms of a continuum scaled Brownian walk (sBw) of polymerization index N. The sBw is an extension of the simple Brownian walk model widely used in path integral calculations of polymer properties, differing from it in containing an additional parameter H (the Hurst index) that can be tuned to produce varying degrees of correlation between adjacent monomers. A calculation using the sBw of the mean time τc for chain closure predicts—under the Wilemski–Fixman approximation for diffusion-controlled reactions—that at early times, τc varies as the 2/3 power of N, in close agreement with the findings of the Brückner et al. study. Other scaling relations of that study, including those related to the probability of loop formation and the mean square displacements of terminal monomers, are also satisfactorily accounted for by the model. [ABSTRACT FROM AUTHOR]

A major limitation of amine-based post-combustion carbon capture technology is the necessity to regenerate amines at high temperatures, which dramatically increases operating costs. This paper concludes the effect of solvent choice as a possible route to modify the thermodynamics and kinetics characterizing the involved amine regeneration reactions and discusses whether these modifications can be economically beneficial. We report experimentally benchmarked computational chemistry calculations of monoethanolamine regeneration reactions employing aqueous and non-aqueous solvents with a wide range of dielectric constants. Unlike previous studies, our improved computational chemistry framework could accurately reproduce the right experimental activation energy of zwitterion formation. From the thermodynamics and kinetics of the predicted reactions, the use of non-aqueous solvents with small dielectric constants led to reductions in regeneration Gibbs free energies, activation barriers, and enthalpy changes. This can reduce energy consumption and give an opportunity to run desorption columns at relatively lower temperatures, thus offering the possibility of relying on low-grade waste heat as an energy input. [ABSTRACT FROM AUTHOR]

As an optimal one-dimensional reaction coordinate, the committor function not only describes the probability of a trajectory initiated at a phase space point first reaching the product state before reaching the reactant state but also preserves the kinetics when utilized to run a reduced dynamics model. However, calculating the committor function in high-dimensional systems poses significant challenges. In this paper, within the framework of milestoning, exact expressions for committor functions at two levels of coarse graining are given, including committor functions of phase space point to point (CFPP) and milestone to milestone (CFMM). When combined with transition kernels obtained from trajectory analysis, these expressions can be utilized to accurately and efficiently compute the committor functions. Furthermore, based on the calculated committor functions, an adaptive algorithm is developed to gradually refine the transition state region. Finally, two model examples are employed to assess the accuracy of these different formulations of committor functions. [ABSTRACT FROM AUTHOR]

In this paper, we demonstrate that the weak temperature dependence of the structure factor of supercooled liquids, a defining feature of the glass transition, is a consequence of the averaging of the scattering intensity due to angular averaging. We show that the speckle at individual wavevectors, calculated from a simulated glass former, exhibits a Debye–Waller factor with a sufficiently large temperature dependence to represent a structural order parameter capable of distinguishing liquid from glass. We also extract from the speckle intensities a quantity proportional to the variance of the local restraint, i.e., a direct experimental measure of the amplitude of structural heterogeneity. [ABSTRACT FROM AUTHOR]

Modeling excited state processes in molecular crystals is relevant for several applications. A popular approach for studying excited state molecular crystals is to use cluster models embedded in point charges. In this paper, we compare the performance of several embedding models in predicting excited states and S1–S0 optical gaps for a set of crystals from the X23 molecular crystal database. The performance of atomic charges based on ground or excited states was examined for cluster models, Ewald embedding, and self-consistent approaches. We investigated the impact of various factors, such as the level of theory, basis sets, embedding models, and the level of localization of the excitation. We consider different levels of theory, including time-dependent density functional theory and Tamm–Dancoff approximation (TDA) (DFT functionals: ωB97X-D and PBE0), CC2, complete active space self-consistent field, and CASPT2. We also explore the impact of selection of the QM region, charge leakage, and level of theory for the description of different kinds of excited states. We implemented three schemes based on distance thresholds to overcome overpolarization and charge leakage in molecular crystals. Our findings are compared against experimental data, G0W0-BSE, periodic TDA, and optimally tuned screened range-separated functionals. [ABSTRACT FROM AUTHOR]

In a recent paper, Becke et al. [J. Chem. Phys. 158, 151103 (2023)] presented a novel double hybrid density functional, "DH23," whose terms are based on good physics. Its 12 coefficients were trained on the GMTKN55 (general main-group thermochemistry, kinetics, and noncovalent interactions) chemical database of Goerigk et al. [Phys. Chem. Chem. Phys. 19, 32184 (2017)]. The lowest GMTKN55 "WTMAD2" error to date for any hybrid or double hybrid density functional was obtained (1.76 kcal/mol). Here, we make some revisions to DH23 and test its efficacy on reference data beyond GMTKN55, namely, organometallic reaction energies and barrier heights. The results confirm that DH23 is robust outside its training set. In the process, a slightly smaller GMTKN55 WTMAD2 of 1.73 kcal/mol is achieved. [ABSTRACT FROM AUTHOR]

TurboGenius is an open-source Python package designed to fully control ab initio quantum Monte Carlo (QMC) jobs using a Python script, which allows one to perform high-throughput calculations combined with TurboRVB [Nakano et al. J. Phys. Chem. 152, 204121 (2020)]. This paper provides an overview of the TurboGenius package and showcases several results obtained in a high-throughput mode. For the purpose of performing high-throughput calculations with TurboGenius, we implemented another open-source Python package, TurboWorkflows, that enables one to construct simple workflows using TurboGenius. We demonstrate its effectiveness by performing (1) validations of density functional theory (DFT) and QMC drivers as implemented in the TurboRVB package and (2) benchmarks of Diffusion Monte Carlo (DMC) calculations for several datasets. For (1), we checked inter-package consistencies between TurboRVB and other established quantum chemistry packages. By doing so, we confirmed that DFT energies obtained by PySCF are consistent with those obtained by TurboRVB within the local density approximation (LDA) and that Hartree–Fock (HF) energies obtained by PySCF and Quantum Package are consistent with variational Monte Carlo energies obtained by TurboRVB with the HF wavefunctions. These validation tests constitute a further reliability check of the TurboRVB package. For (2), we benchmarked the atomization energies of the Gaussian-2 set, the binding energies of the S22, A24, and SCAI sets, and the equilibrium lattice parameters of 12 cubic crystals using DMC calculations. We found that, for all compounds analyzed here, the DMC calculations with the LDA nodal surface give satisfactory results, i.e., consistent either with high-level computational or with experimental reference values. [ABSTRACT FROM AUTHOR]

Polymer models describing the statistics of biomolecules under confinement have applications to a wide range of single-molecule experimental techniques and give insight into biologically relevant processes in vivo. In this paper, we determine the transverse position and bending correlation functions for a wormlike chain confined within slits and cylinders (with one and two confined dimensions, respectively) using a mean-field approach that enforces rigid constraints on average. We show the theoretical predictions accurately capture the statistics of a wormlike chain from Monte Carlo simulations in both confining geometries for both weak and strong confinement. We also show that the longitudinal correlation function is accurately computed for a chain confined to a slit and leverages the accuracy of the model to suggest an experimental technique to infer the (often unobservable) transverse statistics from the (directly observable) longitudinal end-to-end distance. [ABSTRACT FROM AUTHOR]

In this paper, a global and full-dimensional potential energy surface at the 2A″ ground state for the Al + O2 → AlO + O reaction was constructed, for the first time, based on extensive electronic structure calculations using the doubly hybrid functional XYG3 and potential energy surface fittings by neural networks. Details of the reaction paths have been analyzed. It was found that both two intermediates, the cyclic-AlO2 and the linear-OAlO, were able to dissociate to the AlO + O products, and the isomerization process between these two intermediates was controlled by conical intersections between two 2A″ states. Ro-vibrational state resolved integral cross sections have also been calculated at collision energies from 1.0 to 10.0 kcal/mol. The results support the harpooning mechanism in this metal-oxidant-involved reaction. [ABSTRACT FROM AUTHOR]

Instanton-based rate theory is a powerful tool that is used to explore tunneling in many-dimensional systems. Yet, it diverges at the so-called "crossover temperature." Using the uniform semiclassical transmission probability of Kemble [Phys. Rev. 48, 549 (1935)], we showed recently that in one dimension, one might derive a uniform semiclassical instanton rate theory, which has no divergence. In this paper, we generalize this uniform theory to many-dimensional systems. The resulting theory uses the same input as in the previous instanton theory, yet does not suffer from the divergence. The application of the uniform theory to dissipative systems is considered and used to revise Wolynes' well-known analytical expression for the rate [P. G. Wolynes, Phys. Rev. Lett. 47, 968 (1981)] so that it does not diverge at the "crossover temperature." [ABSTRACT FROM AUTHOR]

Chemical and photochemical reactivity, as well as supramolecular organization and several other molecular properties, can be modified by strong interactions between light and matter. Theoretical studies of these phenomena require the separation of the Schrödinger equation into different degrees of freedom as in the Born–Oppenheimer approximation. In this paper, we analyze the electron–photon Hamiltonian within the cavity Born–Oppenheimer approximation (CBOA), where the electronic problem is solved for fixed nuclear positions and photonic parameters. In particular, we focus on intermolecular interactions in representative dimer complexes. The CBOA potential energy surfaces are compared with those obtained using a polaritonic approach, where the photonic and electronic degrees of freedom are treated at the same level. This allows us to assess the role of electron–photon correlation and the accuracy of CBOA. [ABSTRACT FROM AUTHOR]

The fewest switches surface hopping (FSSH) method proposed by Tully in 1990 [Tully, J. Chem. Phys. 93, 1061 (1990)]—along with its many later variations—forms the basis for most practical simulations of molecular dynamics with electronic transitions in realistic systems. Despite its popularity, a rigorous formal derivation of the algorithm has yet to be achieved. In this paper, we derive the energy-conserving momentum jumps employed by FSSH from the perspective of quantum trajectory surface hopping (QTSH) [Martens, J. Phys. Chem. A 123, 1110 (2019)]. In the limit of localized nonadiabatic transitions, simple mathematical and physical arguments allow the FSSH algorithm to be derived from first principles. For general processes, the quantum forces characterizing the QTSH method provide accurate results for nonadiabatic dynamics with rigorous energy conservation, at the ensemble level, within the consistency of the underlying stochastic surface hopping without resorting to the artificial momentum rescaling of FSSH. [ABSTRACT FROM AUTHOR]

In Paper I [H. Francisco, A. C. Cancio, and S. B. Trickey, J. Chem. Phys. 159, 214102 (2023)], we gave a regularization of the Tao–Mo exchange functional that removes the order-of-limits problem in the original Tao–Mo form and also eliminates the unphysical behavior introduced by an earlier regularization while essentially preserving compliance with the second-order gradient expansion. The resulting simplified, regularized (sregTM) functional delivers performance on standard molecular and solid state test sets equal to that of the earlier revised, regularized Tao–Mo functional. Here, we address de-orbitalization of that new sregTM into a pure density functional. We summarize the failures of the Mejía-Rodríguez and Trickey de-orbitalization strategy [Phys. Rev. A 96, 052512 (2017)] when used with both versions. We discuss how those failures apparently arise in the so-called z′ indicator function and in substitutes for the reduced density Laplacian in the parent functionals. Then, we show that the sregTM functional can be de-orbitalized somewhat well with a rather peculiarly parameterized version of the previously used deorbitalizer. We discuss, briefly, a de-orbitalization that works in the sense of reproducing error patterns but that apparently succeeds by cancelation of major qualitative errors associated with the de-orbitalized indicator functions α and z, hence, is not recommended. We suggest that the same issue underlies the earlier finding of comparatively mediocre performance of the de-orbitalized Tao–Perdew–Staroverov–Scuseri functional. Our work demonstrates that the intricacy of such two-indicator functionals magnifies the errors introduced by the Mejía-Rodríguez and Trickey de-orbitalization approach in ways that are extremely difficult to analyze and correct. [ABSTRACT FROM AUTHOR]

The revised, regularized Tao–Mo (rregTM) exchange-correlation density functional approximation (DFA) [A. Patra, S. Jana, and P. Samal, J. Chem. Phys. 153, 184112 (2020) and Jana et al., J. Chem. Phys. 155, 024103 (2021)] resolves the order-of-limits problem in the original TM formulation while preserving its valuable essential behaviors. Those include performance on standard thermochemistry and solid data sets that is competitive with that of the most widely explored meta-generalized-gradient-approximation DFAs (SCAN and r2SCAN) while also providing superior performance on elemental solid magnetization. Puzzlingly however, rregTM proved to be intractable for de-orbitalization via the approach of Mejía-Rodríguez and Trickey [Phys. Rev. A 96, 052512 (2017)]. We report investigation that leads to diagnosis of how the regularization in rregTM of the z indicator functions (z = the ratio of the von-Weizsäcker and Kohn–Sham kinetic energy densities) leads to non-physical behavior. We propose a simpler regularization that eliminates those oddities and that can be calibrated to reproduce the good error patterns of rregTM. We denote this version as simplified, regularized Tao–Mo, sregTM. We also show that it is unnecessary to use rregTM correlation with sregTM exchange: Perdew–Burke–Ernzerhof correlation is sufficient. The subsequent paper shows how sregTM enables some progress on de-orbitalization. [ABSTRACT FROM AUTHOR]

Quantum walks are the quantum counterpart of classical random walks and have various applications in quantum information science. Polar molecules have rich internal energy structure and long coherence time and thus are considered as a promising candidate for quantum information processing. In this paper, we propose a theoretical scheme for implementing discrete-time quantum walks on a circle with dipole–dipole coupled SrO molecules. The states of the walker and the coin are encoded in the pendular states of polar molecules induced by an external electric field. We design the optimal microwave pulses for implementing quantum walks on a four-node circle and a three-node circle by multi-target optimal control theory. To reduce the accumulation of decoherence and improve the fidelity, we successfully realize a step of quantum walk with only one optimal pulse. Moreover, we also encode the walker into a three-level molecular qutrit and a four-level molecular ququart and design the corresponding optimal pulses for quantum walks, which can reduce the number of molecules used. It is found that all the quantum walks on a circle in our scheme can be achieved via optimal control fields with high fidelities. Our results could shed some light on the implementation of discrete-time quantum walks and high-dimensional quantum information processing with polar molecules. [ABSTRACT FROM AUTHOR]

In a binary mixture, stripes refer to a one-dimensional periodicity of the composition, namely, a regular alternation of layers filled with particles of mostly one species. We have recently introduced [Munaò et al., Phys. Chem. Chem. Phys. 25, 16227 (2023)] a model that possibly provides the simplest binary mixture endowed with stripe order. The model consists of two species of identical hard spheres with equal concentration, which mutually interact through a square-well potential. In that paper, we have numerically shown that stripes are present in both liquid and solid phases when the attraction range is rather long. Here, we study the phase behavior of the model in terms of a density functional theory capable to account for the existence of stripes in the dense mixture. Our theory is accurate in reproducing the phases of the model, at least insofar as the composition inhomogeneities occur on length scales quite larger than the particle size. Then, using Monte Carlo simulations, we prove the existence of solid stripes even when the square well is much thinner than the particle diameter, making our model more similar to a real colloidal mixture. Finally, when the width of the attractive well is equal to the particle diameter, we observe a different and more complex form of compositional order in the solid, where each species of particle forms a regular porous matrix holding in its holes the other species, witnessing a surprising variety of emergent behaviors for a very basic model of interaction. [ABSTRACT FROM AUTHOR]

In this paper, we investigate the time evolution of quantum coherence—the off-diagonal elements of the density matrix of a multistate quantum system—from the perspective of the Wigner–Moyal formalism. This approach provides an exact phase space representation of quantum mechanics. We consider the coherent evolution of nuclear wavepackets in a molecule with two electronic states. For harmonic potentials, the problem is analytically soluble for both a fully quantum mechanical description and a semiclassical description. We highlight the serious deficiencies of the semiclassical treatment of coherence for general systems and illustrate how even qualitative accuracy requires higher order terms in the Moyal expansion to be included. The model provides an experimentally relevant example of a molecular Schrödinger's cat state. The alive and dead cats of the exact two-state quantum evolution collapse into a "zombie" cat in the semiclassical limit—an averaged behavior, neither alive nor dead, leading to significant errors. The inclusion of the Moyal correction restores a faithful simultaneously alive and dead representation of the cat that is experimentally observable. [ABSTRACT FROM AUTHOR]

The computation of the nuclear quantum dynamics of molecules is challenging, requiring both accuracy and efficiency to be applicable to systems of interest. Recently, theories have been developed for employing time-dependent basis functions (denoted modals) with vibrational coupled cluster theory (TDMVCC). The TDMVCC method was introduced along with a pilot implementation, which illustrated good accuracy in benchmark computations. In this paper, we report an efficient implementation of TDMVCC, covering the case where the wave function and Hamiltonian contain up to two-mode couplings. After a careful regrouping of terms, the wave function can be propagated with a cubic computational scaling with respect to the number of degrees of freedom. We discuss the use of a restricted set of active one-mode basis functions for each mode, as well as two interesting limits: (i) the use of a full active basis where the variational modal determination amounts essentially to the variational determination of a time-dependent reference state for the cluster expansion; and (ii) the use of a single function as an active basis for some degrees of freedom. The latter case defines a hybrid TDMVCC/TDH (time-dependent Hartree) approach that can obtain even lower computational scaling. The resulting computational scaling for hybrid and full TDMVCC[2] is illustrated for polyaromatic hydrocarbons with up to 264 modes. Finally, computations on the internal vibrational redistribution of benzoic acid (39 modes) are used to show the faster convergence of TDMVCC/TDH hybrid computations towards TDMVCC compared to simple neglect of some degrees of freedom. [ABSTRACT FROM AUTHOR]

Solving the Liouville–von-Neumann equation using a density operator provides a more complete picture of dynamical quantum phenomena than by using a wavepacket and solving the Schrödinger equation. As density operators are not restricted to the description of pure states, they can treat both thermalized and open systems. In practice, however, they are rarely used to study molecular systems as the computational resources required are even more prohibitive than those needed for wavepacket dynamics. In this paper, we demonstrate the potential utility of a scheme based on the powerful multi-layer multi-configurational time-dependent Hartree algorithm for propagating multi-dimensional density operators. Studies of two systems using this method are presented at a range of temperatures and including up to 13 degrees of freedom. The first case is single proton transfer in salicylaldimine, while the second is double proton transfer in porphycene. A comparison is also made with the approach of using stochastic wavepackets. [ABSTRACT FROM AUTHOR]

A model system of identical particles interacting via a hard-sphere potential is essential in condensed matter physics; it helps to understand in and out of equilibrium phenomena in complex fluids, such as colloidal dispersions. Yet, most of the fixed time-step algorithms to study the transport properties of those systems have drawbacks due to the mathematical nature of the interparticle potential. Because of this, mapping a hard-sphere potential onto a soft potential has been recently proposed [Báez et al., J. Chem. Phys. 149, 164907 (2018)]. More specifically, using the second virial coefficient criterion, one can set a route to estimate the parameters of the soft potential that accurately reproduces the thermodynamic properties of a monocomponent hard-sphere system. However, real colloidal dispersions are multicomponent or polydisperse, making it important to find an efficient way to extend the potential model for dealing with such kind of many-body systems. In this paper, we report on the extension and applicability of the second virial coefficient criterion to build a description that correctly captures the phenomenology of both multicomponent and polydisperse hard-sphere dispersions. To assess the accuracy of the continuous potentials, we compare the structure of soft polydisperse systems with their hard-core counterpart. We also contrast the structural and thermodynamic properties of soft binary mixtures with those obtained through mean-field approximations and the Ornstein–Zernike equation for the two-component hard-sphere dispersion. [ABSTRACT FROM AUTHOR]

Effect mechanisms of the undercooling degree and the surface configuration on the ice growth characteristics were revealed under micro-droplets icing conditions. Preferential ice crystals appear firstly on the surfaces due to the randomness of icing, and obtain growth advantages to form protruding structures. Protruding structures block the incoming droplets from contacting the substrates, causing voids around the structures. The undercooling degree mainly affects the density and the growth rate of preferential ice crystals. With the increase of undercooling degree, the preferential ice crystals have higher density and growth rate, resulting in stronger growth advantage and higher porosity. The surface configuration affects the growth mode, and the ice layer grows with uniform mode, spreading mode and structure-induced mode on the aluminum, smooth Polytetrafluoroethylene (PTFE) and rough PTFE surface respectively, causing the needle-like, ridge-like and cluster-like ice crystals. The rough structures effectively improve the porosity of the ice layer, which is beneficial for optimizing the icephobic property of the materials. This paper provides important theoretical guidance for the design of subsequent icephobic materials. [ABSTRACT FROM AUTHOR]

Computing the solubility of crystals in a solvent using atomistic simulations is notoriously challenging due to the complexities and convergence issues associated with free-energy methods, as well as the slow equilibration in direct-coexistence simulations. This paper introduces a molecular-dynamics workflow that simplifies and robustly computes the solubility of molecular or ionic crystals. This method is considerably more straightforward than the state-of-the-art, as we have streamlined and optimised each step of the process. Specifically, we calculate the chemical potential of the crystal using the gas-phase molecule as a reference state, and employ the S0 method to determine the concentration dependence of the chemical potential of the solute. We use this workflow to predict the solubilities of sodium chloride in water, urea polymorphs in water, and paracetamol polymorphs in both water and ethanol. Our findings indicate that the predicted solubility is sensitive to the chosen potential energy surface. Furthermore, we note that the harmonic approximation often fails for both molecular crystals and gas molecules at or above room temperature, and that the assumption of an ideal solution becomes less valid for highly soluble substances. [ABSTRACT FROM AUTHOR]