Search

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]

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

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

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

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

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

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]

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]

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]

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

X-ray photoelectron spectroscopy (XPS) is generally used for chemical analysis of surfaces and interfaces. This method involves the analysis of changes in binding energies and peak shapes of elements under consideration. It is also possible to use XPS to study the effect of x-ray radiation on the electrical properties of thin films. We measured the Si 2p peak using x-ray powers of 300 and 150 W on ∼135 nm silicon dioxide (SiO2) thin films grown on both n- and p-type substrates while applying DC or AC external biases. Using the shifts in the binding energy of the Si 2p peak, we calculated the resistances and the capacitances of the SiO2 thin film. The way that the binding energies of the Si 2p peak and the capacitance of the thin film change as a function of the type of Si substrate and the power of the x-ray are explained using band bending. [ABSTRACT FROM AUTHOR]

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

Experimental demonstrations of polarization-selection two-dimensional Vibrational-Electronic (2D VE) and 2D Electronic-Vibrational (2D EV) spectroscopies aim to map the magnitudes and spatial orientations of coupled electronic and vibrational coordinates in complex systems. The realization of that goal depends on our ability to connect spectroscopic observables with molecular structural parameters. In this paper, we use a model Hamiltonian consisting of two anharmonically coupled vibrational modes in electronic ground and excited states with linear and bilinear vibronic coupling terms to simulate polarization-selective 2D EV and 2D VE spectra. We discuss the relationships between the linear vibronic coupling and two-dimensional Huang–Rhys parameters and between the bilinear vibronic coupling term and Duschinsky mixing. We develop a description of the vibronic transition dipoles and explore how the Hamiltonian parameters and non-Condon effects impact their amplitudes and orientations. Using simulated polarization-selective 2D EV and 2D VE spectra, we show how 2D peak positions, amplitudes, and anisotropy can be used to measure parameters of the vibronic Hamiltonian and non-Condon effects. This paper, along with the first in the series, provides the reader with a detailed description of reading, simulating, and analyzing multimode, polarization-selective 2D EV and 2D VE spectra with an emphasis on extracting vibronic coupling parameters from complex spectra. [ABSTRACT FROM AUTHOR]

Two-dimensional Electronic–Vibrational (2D EV) spectroscopy and two-dimensional Vibrational–Electronic (2D VE) spectroscopy are among the newest additions to the coherent multidimensional spectroscopy toolbox, and they are directly sensitive to vibronic couplings. In this first of two papers, the complete orientational response functions are developed for a model system consisting of two coupled anharmonic oscillators and two electronic states in order to simulate polarization-selective 2D EV and 2D VE spectra with arbitrary combinations of linearly polarized electric fields. Here, we propose analytical methods to isolate desired signals within complicated spectra and to extract the relative orientation between vibrational and vibronic dipole moments of the model system using combinations of polarization-selective 2D EV and 2D VE spectral features. Time-dependent peak amplitudes of coherence peaks are also discussed as means for isolating desired signals within the time-domain. This paper serves as a field guide for using polarization-selective 2D EV and 2D VE spectroscopies to map coupled vibronic coordinates on the molecular frame. [ABSTRACT FROM AUTHOR]

In recent years there has been a rapid growth in the development and application of new stochastic methods in electronic structure. These methods are quite diverse, from many-body wave function techniques in real space or determinant space to being used to sum perturbative expansions. This growth has been spurred by the more favorable scaling with the number of electrons and often better parallelization over large numbers of central processing unit (CPU) cores or graphical processing units (GPUs) than for high-end non-stochastic wave function based methods. This special issue of the Journal of Chemical Physics includes 33 papers that describe recent developments and applications in this area. As seen from the articles in the issue, stochastic electronic structure methods are applicable to both molecules and solids and can accurately describe systems with strong electron correlation. This issue was motivated, in part, by the 2019 Telluride Science Research Center workshop on Stochastic Electronic Structure Methods that we organized. Below we briefly describe each of the papers in the special issue, dividing the papers into six subtopics. [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 this paper, we study electrostatic and structural properties between pH-responsive polyelectrolyte brushes by using a strong stretching theory accounting for excluded volume interactions, the density of polyelectrolyte chargeable sites, and the Born energy difference between the inside and outside of the brush layer. In a free energy framework, we obtain self-consistent field equations to determine electrostatic properties between two pH-responsive polyelectrolyte brushes. We elucidate that in the region between two pH-responsive polyelectrolyte brushes, electrostatic potential at the centerline and osmotic pressure increase not only with excluded volume interaction but also with the density of chargeable sites on a polyelectrolyte molecule. Importantly, we clarify that when two pH-responsive polyelectrolyte brushes approach each other, the brush thickness becomes short and that a large excluded volume interaction and a large density of chargeable sites yield the enhanced contract of polyelectrolyte brushes. In addition, we also demonstrate how the influence of such quantities as pH, the number of Kuhn monomers, the density of charged sites, the lateral separation between adjacent polyelectrolyte brushes, and Kuhn length on the electrostatic and structural properties between the two polyelectrolyte brushes is affected by the exclusion volume interaction. Finally, we investigate the influence of the Born energy difference on the thickness of polyelectrolyte brushes and the osmotic pressure between two pH-responsive polyelectrolyte brushes. [ABSTRACT FROM AUTHOR]

Similar to cells, bacteria, and other micro-organisms, synthetic chemically active colloids can harness the energy from their environment through a surface chemical reaction and use the energy to self-propel in fluidic environments. In this paper, we study the chemo-mechanical coupling that leads to the self-propulsion of chemically active colloids. The coupling between chemical reactions and momentum transport is a consequence of Onsager reciprocal relations. They state that the velocity and the surface reaction rate are related to mechanical and chemical affinities through a symmetric matrix. A consequence of Onsager reciprocal relations is that if a chemical reaction drives the motion of the colloid, then an external force generates a reaction rate. Here, we investigate Onsager reciprocal relations for a spherical active colloid that catalyzes a reversible surface chemical reaction between two species. We solve the relevant transport equations using a perturbation expansion and numerical simulations to demonstrate the validity of reciprocal relations around the equilibrium. Our results are consistent with previous studies and highlight the key role of solute advection in preserving the symmetry of the Onsager matrix. Finally, we show that Onsager reciprocal relations break down around a nonequilibrium steady state, which has implications for the thermal fluctuations of the active colloids used in experiments. [ABSTRACT FROM AUTHOR]

This paper reports the first variationally based predictions of the lowest excited state in diamond (Γ25′ → Γ15) in the unrelaxed (optical) and structurally relaxed (thermal) configurations, from direct Δ-self-consistent-field (SCF) calculations based on B3LYP, PBE0, HSE06, and GGA functionals. For the B3LYP functional, which has the best overall performance, the energy of the optical state, 7.27 eV, is within the observed range of (7.2–7.4) eV and is predicted to be insulating, with indirect bandgaps of (5.6–5.8) eV. Mulliken analyses of the excited state wavefunction indicate extensive redistributions of charge and spin resulting in a strongly excitonic state with a central charge of −0.8ǀeǀ surrounded by charges of +0.12ǀeǀ at the four nearest neighbor sites. The thermally relaxed state is predicted to be similarly excitonic, with comparable bandgaps and atomic charges. Calculations of the ground and excited state relaxations lead to a Stokes shift of 0.47 eV and predicted Γ-point luminescence energy of 6.89 eV. Assuming a similar shift at the band edge (X1), an estimate of 5.29 eV is predicted for the luminescence energy, which compares with the observed value of 5.27 eV. Excited state vibrational spectra show marked differences from the ground state, with the introduction of an infrared peak at 1150 cm−1 and a modest shift of 2 cm−1 in the TO(X) Raman mode at 1340 cm−1. Similar calculations of the lowest energy bi- and triexcitons predict these to be bound states in both optical and thermal configurations and plausible precursors to exciton condensation. Estimates of bi- and triexciton luminescence energies predict red shifts with respect to the single exciton line, which are compared to the recently reported values. [ABSTRACT FROM AUTHOR]

Recently, we proposed a new orbital analysis method, natural reaction orbital (NRO), which automatically extracts orbital pairs that characterize electron transfer in reaction processes by singular value decomposition of the first-order orbital response matrix to the nuclear coordinate displacements [Ebisawa et al., Phys. Chem. Chem. Phys. 24, 3532 (2022)]. NRO analysis along the intrinsic reaction coordinate (IRC) for several typical chemical reactions demonstrated that electron transfer occurs mainly in the vicinity of transition states and in regions where the energy profile along the IRC shows shoulder features, allowing the reaction mechanism to be explained in terms of electron motion. However, its application has been limited to single configuration theories such as Hartree–Fock theory and density functional theory. In this work, the concept of NRO is extended to multiconfigurational wavefunctions and formulated as the multiconfiguration NRO (MC-NRO). The MC-NRO method is applicable to various types of electronic structure theories, including multiconfigurational theory and linear response theory, and is expected to be a practical tool for extracting the essential qualitative features of a broad range of chemical reactions, including covalent bond dissociation and chemical reactions in electronically excited states. In this paper, we calculate the IRC for five basic chemical reaction processes at the level of the complete active space self-consistent field theory and discuss the phenomenon of electron transfer by performing MC-NRO analysis along each IRC. Finally, issues and future prospects of the MC-NRO method are discussed. [ABSTRACT FROM AUTHOR]

This paper combines the valence bond block diabatization approach with the idea of orbital breathing. With highly compact wave functions, the breathing orbital valence bond (BOVB) method is applied to investigate several atomic and molecular properties, including the electron affinity of F, the adiabatic and diabatic potential energy curves and the dipole moment curves of the two lowest-lying 1Σ+ states, the electronic coupling curve and the crossing distance of the two diabatic states, and the spectroscopic constants of the ground states for LiF. The configuration selection scheme proposed in this work is quite general, requiring only the selection of several de-excitation and excitation orbitals in a sense like the restricted active space self-consistent field method. Practically, this is also the first time that BOVB results are extrapolated to complete basis set limit. Armed with the chemical intuition provided by valence bond theory, the classic but challenging covalent-ionic interaction in the title molecule is not only conceptually interpreted but also accurately computed. [ABSTRACT FROM AUTHOR]

The multireference second-order perturbation theory (CASPT2) is known to deliver a quantitative description of various complex electronic states. Despite its near-size-consistent nature, the applicability of the CASPT2 method to large, real-life systems is mostly hindered by large computational and storage costs for the two-external tensors, such as two-electron integrals, amplitudes, and residuum. To this end, Menezes and co-workers developed a reduced-scaling CASPT2 scheme by incorporating the local pair-natural orbital (PNO) representation of the many-body wave functions using non-orthonormal projected atomic orbitals (PAOs) into the CASPT theory [F. Menezes et al., J. Chem. Phys. 145, 124115 (2016)]. Alternatively, in this paper, we develop a new PNO-based CASPT2 scheme using the orthonormal localized virtual molecular orbitals (LVMOs) and assess its performance and accuracy in comparison with the conventional PAO-based counterpart. Albeit the compactness, the LVMOs were considered to perform somewhat poorly compared to PAOs in the local correlation framework because they caused enormously large orbital domains. In this work, we show that the size of LVMO domains can be rendered comparable to or even smaller than that of PAOs by the use of the differential overlap integrals for domain construction. Optimality of the MOs from the CASSCF treatment is a key to reducing the LVMO domain size for the multireference case. Due to the augmented Hessian-based localization algorithm, an additional computational cost for obtaining the LVMOs is relatively minor. We demonstrate that the LVMO-based PNO-CASPT2 method is routinely applicable to large, real-life molecules such as Menshutkin SN2 reaction in a single-walled carbon nanotube reaction field. [ABSTRACT FROM AUTHOR]

This paper describes Wick&d, an implementation of the algebra of second-quantized operators normal ordered with respect to general correlated references and the corresponding Wick theorem [D. Mukherjee, Chem. Phys. Lett. 274, 561 (1997) and W. Kutzelnigg and D. Mukherjee, J. Chem. Phys. 107, 432 (1997)]. Wick&d employs a compact representation of operators and a backtracking algorithm to efficiently evaluate Wick contractions. Since Wick&d can handle both fully and partially contracted terms, it can be applied to both projective and Fock-space many-body formalisms. To demonstrate the usefulness of Wick&d, we use it to evaluate the single-reference coupled cluster equations up to octuple excitations and report an automated derivation and implementation of the second-order driven similarity renormalization group multireference perturbation theory. [ABSTRACT FROM AUTHOR]

Astrophysical modeling of processes in environments that are not in local thermal equilibrium requires the knowledge of state-to-state rate coefficients of rovibrational transitions in molecular collisions. These rate coefficients can be obtained from coupled-channel (CC) quantum scattering calculations, which are very demanding, however. Here, we present various approximate but more efficient methods based on the coupled-states approximation (CSA), which neglects the off-diagonal Coriolis coupling in the scattering Hamiltonian in body-fixed coordinates. In particular, we investigated a method called NNCC (nearest-neighbor Coriolis coupling) [Yang et al., J. Chem. Phys. 148, 084101 (2018)] that includes Coriolis coupling to first order. The NNCC method is more demanding than the common CSA method but still much more efficient than full CC calculations, and it is substantially more accurate than CSA. All of this is illustrated by showing state-to-state cross sections and rate coefficients of rovibrational transitions induced in CO2 by collisions with He atoms. It is also shown that a further reduction of CPU time, practically without loss of accuracy, can be obtained by combining the NNCC method with the multi-channel distorted-wave Born approximation that we applied in full CC calculations in a previous paper. [ABSTRACT FROM AUTHOR]

Modeling and understanding properties of materials from first principles require knowledge of the underlying atomistic structure. This entails knowing the individual chemical identity and position of all atoms involved. Obtaining such information for macro-molecules, nano-particles, and clusters and for the surface, interface, and bulk phases of amorphous and solid materials represents a difficult high-dimensional global optimization problem. The rise of machine learning techniques in materials science has, however, led to many compelling developments that may speed up structure searches. The complexity of such new methods has prompted a need for an efficient way of assembling them into global optimization algorithms that can be experimented with. In this paper, we introduce the Atomistic Global Optimization X (AGOX) framework and code as a customizable approach that enables efficient building and testing of global optimization algorithms. A modular way of expressing global optimization algorithms is described, and modern programming practices are used to enable that modularity in the freely available AGOX Python package. A number of examples of global optimization approaches are implemented and analyzed. This ranges from random search and basin-hopping to machine learning aided approaches with on-the-fly learnt surrogate energy landscapes. The methods are applied to problems ranging from supported clusters over surface reconstructions to large carbon clusters and metal-nitride clusters incorporated into graphene sheets. [ABSTRACT FROM AUTHOR]

The aging rate of glasses has traditionally been modeled as a function of temperature, T, and fictive temperature, while density, ρ, is not explicitly included as a parameter. However, this description does not naturally connect to the modern understanding of what governs the relaxation rate in equilibrium. In equilibrium, it is well known that the relaxation rate, γeq, depends on temperature and density. In addition, a large class of systems obeys density scaling, which means the rate specifically depends on the scaling parameter, Γ = e(ρ)/T, where e(ρ) is a system specific function. This paper presents a generalization of the fictive temperature concept in terms of a fictive scaling parameter, Γfic, and a density scaling conjecture for aging glasses in which the aging rate depends on Γ and Γfic. [ABSTRACT FROM AUTHOR]

Spectral line shapes provide a window into the local environment coupled to a quantum transition in the condensed phase. In this paper, we build upon a stochastic model to account for non-stationary background processes produced by broad-band pulsed laser stimulation, as distinguished from those for stationary phonon bath. In particular, we consider the contribution of pair-fluctuations arising from the full bosonic many-body Hamiltonian within a mean-field approximation, treating the coupling to the system as a stochastic noise term. Using the Itô transformation, we consider two limiting cases for our model, which lead to a connection between the observed spectral fluctuations and the spectral density of the environment. In the first case, we consider a Brownian environment and show that this produces spectral dynamics that relax to form dressed excitonic states and recover an Anderson–Kubo-like form for the spectral correlations. In the second case, we assume that the spectrum is Anderson–Kubo like and invert to determine the corresponding background. Using the Jensen inequality, we obtain an upper limit for the spectral density for the background. The results presented here provide the technical tools for applying the stochastic model to a broad range of problems. [ABSTRACT FROM AUTHOR]

Newly developed coupled-cluster (CC) methods enable simulations of ionization potentials and spectral functions of molecular systems in a wide range of energy scales ranging from core-binding to valence. This paper discusses the results obtained with the real-time equation-of-motion CC cumulant (RT-EOM-CC) approach and CC Green's function (CCGF) approaches in applications to the water and water dimer molecules. We compare the ionization potentials obtained with these methods for the valence region with the results obtained with the coupled-cluster with singles, doubles, and perturbative triples formulation as a difference of energies for N and N − 1 electron systems. All methods show good agreement with each other. They also agree well with the experiment with errors usually below 0.1 eV for the ionization potentials. We also analyze unique features of the spectral functions, associated with the position of satellite peaks, obtained with the RT-EOM-CC and CCGF methods employing single and double excitations, as a function of the monomer OH bond length and the proton transfer coordinate in the dimer. Finally, we analyze the impact of the basis set effects on the quality of calculated ionization potentials and find that the basis set effects are less pronounced for the augmented-type sets. [ABSTRACT FROM AUTHOR]