Influence of nuclear quantum effects on the electronic properties of amorphous carbon.

We carry out quantum simulations to study the physical properties of diamond-like amorphous carbon by coupling first-principles molecular dynamics with a quantum thermostat, and we analyze multiple samples representative of different defective sites present in the disordered network. We show that quantum vibronic coupling is critical in determining the electronic properties of the system, in particular its electronic and mobility gaps, while it has a moderate influence on the structural properties. We find that despite localized electronic states near the Fermi level, the quantum nature of the nuclear motion leads to a renormalization of the electronic gap surprisingly similar to that found in crystalline diamond. We also discuss the notable influence of nuclear quantum effects on band-like and variable-hopping mechanisms contributing to electrical conduction. Our calculations indicate that methods often used to evaluate electron-phonon coupling in ordered solids are inaccurate to study the electronic and transport properties of amorphous semiconductors composed of light atoms.

Nuclear Quantum Effects Accelerate Charge Recombination but Boost the Stability of Inorganic Perovskites in Mild Humidity.

Experiments have demonstrated that mild humidity can enhance the stability of the CsPbBr3 perovskite, though the underlying mechanism remains unclear. Utilizing ab initio molecular dynamics, ring polymer molecular dynamics, and non-adiabatic molecular dynamics, our study reveals that nuclear quantum effects (NQEs) play a crucial role in stabilizing the lattice rigidity of the perovskite while simultaneously shortening the charge carrier lifetime. NQEs reduce the extent of geometric disorder and the number of atomic fluctuations, diminish the extent of hole localization, and thereby improve the electron-hole overlap and non-adiabatic coupling. Concurrently, these effects significantly suppress phonon modes and slow decoherence. As a result, these factors collectively accelerate charge recombination by a factor of 1.42 compared to that in scenarios excluding NQEs. The resulting sub-10 ns recombination time scales align remarkably well with experimental findings. This research offers novel insight into how moisture resistance impacts the stability and charge carrier lifetime in all-inorganic perovskites.

The Hydrogen-Bond Continuum in the Salt/Cocrystal Systems of Quinoline and Chloro-Nitrobenzoic Acids.

This study investigates the hydrogen-bond geometry in six two-component solid systems composed of quinoline and chloro-nitrobenzoic acids. New X-ray diffraction studies were conducted using both the conventional independent-atom model and the more recent Hirshfeld atom-refinement method, with the latter providing precise hydrogen-atom positions. The systems can be divided into salts (the hydrogen atom transferred to the quinoline nitrogen), cocrystals (the hydrogen atom retained by the acid), and intermediate structures. Solid-state NMR experiments corroborated the X-ray diffraction-derived H-N distances. DFT calculations, using five functionals including hybrid B3LYP and PBE0, showed varying energy profiles for the hydrogen bonds, with notable differences across functionals. These calculations revealed different preferences for salt or cocrystal structures, depending on the functional used. Path-integral molecular dynamics simulations incorporating nuclear quantum effects demonstrated significant hydrogen-atom delocalization, forming a hydrogen-bond continuum, and provided average N-H distances in excellent agreement with experimental results. This comprehensive experimental and theoretical approach highlights the complexity of multicomponent solids. The study emphasizes that the classification into salts or cocrystals is frequently inadequate, as the hydrogen atom is often significantly delocalized in the hydrogen bond. This insight is crucial for understanding and predicting the behavior of such systems in pharmaceutical applications.

A nuclear configuration interaction approach to study nuclear spin effects: an application to ortho- and para-3 He2 @C60.

We introduce a non-orthogonal configuration interaction approach to investigate nuclear quantum effects on energies and densities of confined fermionic nuclei. The Hamiltonian employed draws parallels between confined systems and many-electron atoms, where effective non-Coulombic potentials represent the interactions of the trapped particles. One advantage of this method is its generality, as it offers the potential to study the nuclear quantum effects of various confined species affected by effective isotropic or anisotropic potentials. As a first application, we analyze the quantum states of two 3 He atoms encapsulated in C60 . At the Hartree-Fock level, we observe the breaking of spin and spatial symmetries. To ensure wavefunctions with the correct symmetries, we mix the broken-symmetry Hartree-Fock states within the non-orthogonal configuration interaction expansion. Our proposed approach predicts singly and triply degenerate ground states for the singlet (para-3 He2 @C60 ) and triplet (ortho-3 He2 @C60 ) nuclear spin configurations, respectively. The ortho-3 He2 @C60 ground state is 5.69 cm-1 higher in energy than the para-3 He2 @C60 ground state. The nuclear densities obtained for these states exhibit the icosahedral symmetry of the C60 embedding potential. Importantly, our calculated energies for the lowest 85 states are in close agreement with perturbation theory results based on a harmonic oscillator plus rigid rotor model of 3 He2 @C60 .

Photosynthesis tunes quantum-mechanical mixing of electronic and vibrational states to steer exciton energy transfer.

Photosynthetic species evolved to protect their light-harvesting apparatus from photoxidative damage driven by intracellular redox conditions or environmental conditions. The Fenna-Matthews-Olson (FMO) pigment-protein complex from green sulfur bacteria exhibits redox-dependent quenching behavior partially due to two internal cysteine residues. Here, we show evidence that a photosynthetic complex exploits the quantum mechanics of vibronic mixing to activate an oxidative photoprotective mechanism. We use two-dimensional electronic spectroscopy (2DES) to capture energy transfer dynamics in wild-type and cysteine-deficient FMO mutant proteins under both reducing and oxidizing conditions. Under reducing conditions, we find equal energy transfer through the exciton 4-1 and 4-2-1 pathways because the exciton 4-1 energy gap is vibronically coupled with a bacteriochlorophyll-a vibrational mode. Under oxidizing conditions, however, the resonance of the exciton 4-1 energy gap is detuned from the vibrational mode, causing excitons to preferentially steer through the indirect 4-2-1 pathway to increase the likelihood of exciton quenching. We use a Redfield model to show that the complex achieves this effect by tuning the site III energy via the redox state of its internal cysteine residues. This result shows how pigment-protein complexes exploit the quantum mechanics of vibronic coupling to steer energy transfer.

Using Instanton Theory to Study Quantum Effects in Photosensitization.

Electronic excitation is usually accomplished using light (photoexcitation) and is a key step in a vast number of important physical and biological processes. However, in instances where photoexcitation is not possible, a photosensitizer can excite the target molecule in a process called photosensitization. Unfortunately, full details of its mechanism are still unknown. This perspective gives an overview of the current understanding of photosensitization and describes how instanton theory can be used to fill the gaps, especially with regard tothe importance of quantum tunnelling effects.

The Fate of the Formic Acid Proton on the Anatase TiO2(101) Surface.

As a prototype adsorption reaction of gas Brønsted acid on oxides, we study the adsorption of formic acid on anatase. We perform infrared spectroscopy measurements of adsorbed HCOOH and HCOOD on TiO2 nanopowders, from 13 K up to room temperature in an ultra-high vacuum chamber. We assign the IR signals via computed spectra from nuclear quantum dynamics simulations using our divide-and-conquer semiclassical ab initio molecular dynamics method. The acid proton forms an extraordinarily short and strong hydrogen bond with the surface oxygen. The strength of this hydrogen bond, that compares to H bonds in ice at high pressures, is at the root of a substantial redshift with respect to the typical free OH stretching frequency, which eludes its straightforward detection.

Attachment of Hydrogen Molecules to Atomic Ions (Na+ , Cl- ): Examination of an Adiabatic Separation of the H2 Rotational Motion.

Interactions between molecular hydrogen and ions are of interest in cluster science, astrochemistry and hydrogen storage. In dynamical simulations, H2 molecules are usually modelled as point particles, an approximation that can fail for anisotropic interactions. Here, we apply an adiabatic separation of the H2 rotational motion to build effective pseudoatom-ion potentials and in turn study the properties of (H2 )n Na+ /Cl- clusters. These interaction potentials are based on high-level ab initio calculations and Improved Lennard-Jones parametrizations, while the subsequent dynamics has been performed by quantum Monte Carlo calculations. By comparisons with simulations explicitly describing the molecular rotations, it is concluded that the present adiabatic model is very adequate. Interestingly, we find differences in the cluster stabilities and coordination shells depending on the spin isomer considered (para- or ortho-H2 ), especially for the anionic clusters.

Photodetachment of Deprotonated R-Mandelic Acid: The Role of Proton Delocalization on the Radical Stability.

The photodetachment and stability of R-Mandelate, the deprotonated form of the R-Mandelic acid, was investigated by observing the neutral species issued from either simple photodetachment or dissociative photodetachment in a cold anions set-up. R-Mandalate has the possibility to form an intramolecular ionic hydrogen-bond between adjacent hydroxyl and carboxylate groups. The potential energy surface along the proton transfer (PT) coordinate between both groups (O- …H+ …- OCO) features a single local minima, with the proton localized on the O- group (OH…- OCO). However, the structure with the proton localized on the - OCO group (O- …HOCO) is also observed because it falls within the extremity of the vibrational wavefunction of the OH…- OCO isomer along the PT coordinate. The stability of the corresponding radicals, produced upon photodetachment, is strongly dependent on the position of the proton in the anion: the radicals produced from the OH…- OCO isomer decarboxylate without barrier, while the radicals produced from the O- …HOCO isomer are stable.

Ring-polymer Molecular Dynamics Simulation for the Adsorption of H2 on Ice Clusters (H2 O)n (n=8, 10, and 12).

In the interstellar medium, the H2 adsorption and desorption on the solid water ice are crucial for chemical and physical processes. We have recently investigated the probabilities of H2 sticking on the (H2 O)8 ice, which has quadrilateral surfaces. We have extended the previous work using classical MD and ring-polymer molecular dynamics (RPMD) simulations to the larger ice clusters, (H2 O)10 and (H2 O)12 , which have pentagonal and hexagonal surfaces, respectively. The H2 sticking probabilities decreased as the temperature increased for both cluster cases, whereas the cluster-size-independent profiles were observed. It is thought that the size independence of the probabilities is qualitatively understood from the similar binding energies for all the three cluster systems. Furthermore, the RPMD sticking probabilities are smaller than the classical ones because of the reduction in the binding energies owing to nuclear quantum effects, such as vibrational quantization.

Microtubules as a potential platform for energy transfer in biological systems: a target for implementing individualized, dynamic variability patterns to improve organ function.

Variability characterizes the complexity of biological systems and is essential for their function. Microtubules (MTs) play a role in structural integrity, cell motility, material transport, and force generation during mitosis, and dynamic instability exemplifies the variability in the proper function of MTs. MTs are a platform for energy transfer in cells. The dynamic instability of MTs manifests itself by the coexistence of growth and shortening, or polymerization and depolymerization. It results from a balance between attractive and repulsive forces between tubulin dimers. The paper reviews the current data on MTs and their potential roles as energy-transfer cellular structures and presents how variability can improve the function of biological systems in an individualized manner. The paper presents the option for targeting MTs to trigger dynamic improvement in cell plasticity, regulate energy transfer, and possibly control quantum effects in biological systems. The described system quantifies MT-dependent variability patterns combined with additional personalized signatures to improve organ function in a subject-tailored manner. The platform can regulate the use of MT-targeting drugs to improve the response to chronic therapies. Ongoing trials test the effects of this platform on various disorders.

Nuclear Quantum Effects on the Charge-Density Wave Transition in NbX2 (X = S, Se).

Understanding the origin of charge-density wave (CDW) instability is important for manipulating novel collective electronic states. Many layered transition metal dichalcogenides (TMDs) share similarity in the structural and electronic instability, giving rise to diverse CDW phases and superconductivity. It is still puzzling that even isostructural and isoelectronic TMDs show distinct CDW features. For instance, bulk NbSe2 exhibits CDW order at low temperature, while bulk NbS2 displays no CDW instability. The CDW transitions in single-layer NbS2 and NbSe2 are also different. In the classic limit, we investigate the electron correlation effects on the dimensionality dependence of the CDW ordering. By performing ab initio path integral molecular dynamics simulations and comparative analyses, we further revealed significant nuclear quantum effects in these systems. Specifically, the quantum motion of sulfur anions significantly reduces the CDW transition temperature in both bulk and single-layer NbS2, resulting in distinct CDW features in the NbS2 and NbSe2 systems.

Characterization and Application of Supramolecular Junctions.

The convergence of supramolecular chemistry and single-molecule electronics offers a new perspective on supramolecular electronics, and provides a new avenue toward understanding and application of intermolecular charge transport at the molecular level. In this review, we will provide an overview of the advances in the characterization technique for the investigation of intermolecular charge transport, and summarize the experimental investigation of several non-covalent interactions, including π-π stacking interactions, hydrogen bonding, host-guest interactions and σ-σ interactions at the single-molecule level. We will also provide a perspective on supramolecular electronics and discuss the potential applications and future challenges.

Theoretical perspectives on non-Born-Oppenheimer effects in chemistry.

The Born-Oppenheimer approximation, which assumes that the electrons respond instantaneously to the motion of the nuclei, breaks down for a wide range of chemical and biological processes. The rate constants of such nonadiabatic processes can be calculated using analytical theories, and the real-time nonequilibrium dynamics can be described using numerical atomistic simulations. The selection of an approach depends on the desired balance between accuracy and efficiency. The computational expense of generating potential energy surfaces on-the-fly often favours the use of approximate, robust and efficient methods such as trajectory surface hopping for large, complex systems. The development of formally exact non-Born-Oppenheimer methods and the exploration of well-defined approximations to such methods are critical for providing benchmarks and preparing for the next generation of faster computers. Thus, the parallel development of rigorous but computationally expensive methods and more approximate but computationally efficient methods is optimal. This Perspective briefly summarizes the available theoretical and computational non-Born-Oppenheimer methods and presents examples illustrating how analytical theories and nonadiabatic dynamics simulations can elucidate the fundamental principles of chemical and biological processes. These examples also highlight how theoretical calculations are able to guide the interpretation of experimental data and provide experimentally testable predictions for nonadiabatic processes. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.

Impact of nuclear quantum effects on the structural inhomogeneity of liquid water.

The 2D Raman-terahertz (THz) response of liquid water is studied in dependence of temperature and isotope substitution ([Formula: see text]O, [Formula: see text]O, and [Formula: see text]O). In either case, a very short-lived (i.e., between 75 and 95 fs) echo is observed that reports on the inhomogeneity of the low-frequency intermolecular modes and hence, on the heterogeneity of the hydrogen bond networks of water. The echo lifetime slows down by about 20% when cooling the liquid from room temperature to the freezing point. Furthermore, the echo lifetime of [Formula: see text]O is [Formula: see text] slower than that of [Formula: see text]O, and both can be mapped on each other by introducing an effective temperature shift of [Formula: see text] K. In contrast, the temperature-dependent echo lifetimes of [Formula: see text]O and [Formula: see text]O are the same within error. [Formula: see text]O and [Formula: see text]O have identical masses, yet [Formula: see text]O is much closer to [Formula: see text]O in terms of nuclear quantum effects. It is, therefore, concluded that the echo is a measure of the structural inhomogeneity of liquid water induced by nuclear quantum effects.

Ab initio thermodynamics of liquid and solid water.

Thermodynamic properties of liquid water as well as hexagonal (Ih) and cubic (Ic) ice are predicted based on density functional theory at the hybrid-functional level, rigorously taking into account quantum nuclear motion, anharmonic fluctuations, and proton disorder. This is made possible by combining advanced free-energy methods and state-of-the-art machine-learning techniques. The ab initio description leads to structural properties in excellent agreement with experiments and reliable estimates of the melting points of light and heavy water. We observe that nuclear-quantum effects contribute a crucial [Formula: see text] to the stability of ice Ih, making it more stable than ice Ic. Our computational approach is general and transferable, providing a comprehensive framework for quantitative predictions of ab initio thermodynamic properties using machine-learning potentials as an intermediate step.

Cosmological particle production: a review.

This article will review quantum particle creation in expanding universes. The emphasis will be on the basic physical principles and on selected applications to cosmological models. The needed formalism of quantum field theory in curved spacetime will be summarized, and applied to the example of scalar particle creation in a spatially flat Universe. Estimates for the creation rate will be given and applied to inflationary cosmology models. Analog models which illustrate the same physical principles and may be experimentally realizable are also discussed.

Nuclear quantum effect with pure anharmonicity and the anomalous thermal expansion of silicon.

Despite the widespread use of silicon in modern technology, its peculiar thermal expansion is not well understood. Adapting harmonic phonons to the specific volume at temperature, the quasiharmonic approximation, has become accepted for simulating the thermal expansion, but has given ambiguous interpretations for microscopic mechanisms. To test atomistic mechanisms, we performed inelastic neutron scattering experiments from 100 K to 1,500 K on a single crystal of silicon to measure the changes in phonon frequencies. Our state-of-the-art ab initio calculations, which fully account for phonon anharmonicity and nuclear quantum effects, reproduced the measured shifts of individual phonons with temperature, whereas quasiharmonic shifts were mostly of the wrong sign. Surprisingly, the accepted quasiharmonic model was found to predict the thermal expansion owing to a large cancellation of contributions from individual phonons.

Deuterium Isotope Effects on Acid-Base Equilibrium of Organic Compounds.

Deuterium isotope effects on acid-base equilibrium have been investigated using a combined path integral and free-energy perturbation simulation method. To understand the origin of the linear free-energy relationship of ΔpKa=pKaD2O-pKaH2O versus pKaH2O, we examined two theoretical models for computing the deuterium isotope effects. In Model 1, only the intrinsic isotope exchange effect of the acid itself in water was included by replacing the titratable protons with deuterons. Here, the dominant contribution is due to the difference in zero-point energy between the two isotopologues. In Model 2, the medium isotope effects are considered, in which the free energy change as a result of replacing H2O by D2O in solute-solvent hydrogen-bonding complexes is determined. Although the average ΔpKa change from Model 1 was found to be in reasonable agreement with the experimental average result, the pKaH2O dependence of the solvent isotope effects is absent. A linear free-energy relationship is obtained by including the medium effect in Model 2, and the main factor is due to solvent isotope effects in the anion-water complexes. The present study highlights the significant roles of both the intrinsic isotope exchange effect and the medium solvent isotope effect.

Naphthazarin Derivatives in the Light of Intra- and Intermolecular Forces.

Our long-term investigations have been devoted the characterization of intramolecular hydrogen bonds in cyclic compounds. Our previous work covers naphthazarin, the parent compound of two systems discussed in the current work: 2,3-dimethylnaphthazarin (1) and 2,3-dimethoxy-6-methylnaphthazarin (2). Intramolecular hydrogen bonds and substituent effects in these compounds were analyzed on the basis of Density Functional Theory (DFT), Møller-Plesset second-order perturbation theory (MP2), Coupled Clusters with Singles and Doubles (CCSD) and Car-Parrinello Molecular Dynamics (CPMD). The simulations were carried out in the gas and crystalline phases. The nuclear quantum effects were incorporated a posteriori using the snapshots taken from ab initio trajectories. Further, they were used to solve a vibrational Schrödinger equation. The proton reaction path was studied using B3LYP, ωB97XD and PBE functionals with a 6-311++G(2d,2p) basis set. Two energy minima (deep and shallow) were found, indicating that the proton transfer phenomena could occur in the electronic ground state. Next, the electronic structure and topology were examined in the molecular and proton transferred (PT) forms. The Atoms In Molecules (AIM) theory was employed for this purpose. It was found that the hydrogen bond is stronger in the proton transferred (PT) forms. In order to estimate the dimers' stabilization and forces responsible for it, the Symmetry-Adapted Perturbation Theory (SAPT) was applied. The energy decomposition revealed that dispersion is the primary factor stabilizing the dimeric forms and crystal structure of both compounds. The CPMD results showed that the proton transfer phenomena occurred in both studied compounds, as well as in both phases. In the case of compound 2, the proton transfer events are more frequent in the solid state, indicating an influence of the environmental effects on the bridged proton dynamics. Finally, the vibrational signatures were computed for both compounds using the CPMD trajectories. The Fourier transformation of the autocorrelation function of atomic velocity was applied to obtain the power spectra. The IR spectra show very broad absorption regions between 700 cm-1-1700 cm-1 and 2300 cm-1-3400 cm-1 in the gas phase and 600 cm-1-1800 cm-1 and 2200 cm-1-3400 cm-1 in the solid state for compound 1. The absorption regions for compound 2 were found as follows: 700 cm-1-1700 cm-1 and 2300 cm-1-3300 cm-1 for the gas phase and one broad absorption region in the solid state between 700 cm-1 and 3100 cm-1. The obtained spectroscopic features confirmed a strong mobility of the bridged protons. The inclusion of nuclear quantum effects showed a stronger delocalization of the bridged protons.