Search

Proton-transfer mechanisms and hydration states were investigated in chitin films possessing the functionality of fuel-cell electrolytes. The absolute hydration number per chitin molecule (N) as a function of relative humidity (RH) was determined from the OH stretching bands of H2O molecules, and the proton conductivity was found to enhance above N = 2 (80%RH). The FIR spectrum at 500–900 cm−1 for 20%RH (N < 1) together with first-principles calculations clearly shows that the w1 site has the same hydration strength as the w2 site. The molecular dynamics simulations for N = 2 demonstrate that H2O molecules with tiny fluctuations are localized on w1 and w2, and the hydrogen-bond (HB) network is formed via the CH2OH group of chitin molecules. Shrinkage of the O–O distance (dOO), which synchronizes with the barrier height, is required for proton transfer from H3O+ to adjacent CH2OH groups or H2O molecules. Nevertheless, dOO is hardly modulated for N = 2 because H2O molecules are strongly constrained on w1 and w2, and therefore, the transfer probability becomes small. For N = 3, novel HBs emerged between the additional H2O molecules broadly distributed on the w3 site and H2O molecules on w1 and w2. The transfer probability is enhanced because large fluctuations and diffusions in the whole H2O molecule yield large modulations of dOO. Consequently, long-range proton hopping is driven by the Zundel-type protonated hydrates in the water network. [ABSTRACT FROM AUTHOR]

The behavior of hydration water around familial Parkinson's disease linked mutants of α-synuclein may be linked to the early-onset of Parkinson's disease. For the first time, this study compares the local structure and dynamics of hydration water around different segments of some of the natural mutants of α-synuclein, i.e., E46K, G51D, A30P, and A53E, with that of the wild-type protein through explicit water MD simulations. The results show that the C-terminal segments of the fast aggregating mutants such as E46K and A30P are less exposed to water, while those of the slow aggregating ones such as A53E and G51D are more exposed to water relative to that of the wild-type protein. In addition, the water molecules are found to be more ordered around the C-terminal segment of the A53E and G51D mutants as compared to the wild-type protein. This is due to an increase in the overall charge of α-syn upon A53E and G51D mutations. The translational and rotational motions of water molecules in the hydration shell of the C-terminal segment of A53E and G51D mutants are found to be faster relative to that of the wild-type protein. This study validates the differential hydration environment around the C-terminal segment for the causative and protective mutants of α-synuclein. [ABSTRACT FROM AUTHOR]

Obtaining sub-chemical accuracy (1 kJ mol−1) for reaction energies of medium-sized gas-phase molecules is a longstanding challenge in the field of thermochemical modeling. The perturbative triples correction to coupled-cluster single double triple [CCSD(T)] constitutes an important component of all high-accuracy composite model chemistries that obtain this accuracy but can be a roadblock in the calculation of medium to large systems due to its O ( N 7 ) scaling, particularly in HEAT-like model chemistries that eschew separation of core and valence correlation. This study extends the work of Lesiuk [J. Chem. Phys. 156, 064103 (2022)] with new approximate methods and assesses the accuracy of five different approximations of (T) in the context of a subset of molecules selected from the W4-17 dataset. It is demonstrated that all of these approximate methods can achieve sub-0.1 kJ mol−1 accuracy with respect to canonical, density-fitted (T) contributions with a modest number of projectors. The approximation labeled Z ̃ T appears to offer the best trade-off between cost and accuracy and shows significant promise in an order-of-magnitude reduction in the computational cost of the CCSD(T) component of high-accuracy model chemistries. [ABSTRACT FROM AUTHOR]

In the conventional GW method, the three-point vertex function (Γ) is approximated to unity (Γ ∼ 1). Here, we developed an all-electron first-principles GWΓ method beyond a conventional GW method by considering a first-order three-point vertex function (Γ(1) = 1 + iGGW) in a one-electron self-energy operator. We applied the GWΓ method to simulate the binding energies (BEs) of B1s, C1s, N1s, O1s, and F1s for 19 small-sized molecules. Contrary to the one-shot GW method [or G0W0(LDA)], which underestimates the experimentally determined absolute BEs by about 3.7 eV for B1s, 5.1 eV for C1s, 6.9 eV for N1s, 7.8 eV for O1s, and 5.8 eV for F1s, the GWΓ method successfully reduces these errors by approximately 1–2 eV for all the elements studied here. Notably, the first-order three-point vertex corrections are more significant for heavier elements, following the order of F > O > N > C > B1s. Finally, the computational cost analysis revealed that one term in the GWΓ one-electron self-energy operator, despite being computationally intensive, contributes negligibly (< 0.1 eV) to the C1s, N1s, O1s, and F1s. [ABSTRACT FROM AUTHOR]

Nonadiabatic quantum dynamics are carried out to illustrate the photoionized spectrum of the cyanopropyne (CH3–C≡C–C≡N) as reported in recent experimental measurements [Lamarre et al., J. Mol. Spectrosc. 315, 206 (2015)]. A detailed electronic structure calculation is performed to analyze the topographical details of the first five ionized states, of which three are degenerate states ( X ̃ 2 E , B ̃ 2 E , and C ̃ 2 E) and two are non-degenerate states ( A ̃ 2 A 1 and D ̃ 2 A 1 ). The degenerate E states of the C3V symmetry molecule are prone to Jahn–Teller (JT) instability, and in addition, symmetry allowed A1 − E vibronic coupling, i.e., pseudo-Jahn–Teller (PJT), effects are expected to have a significant impact in the detailed vibronic structure of these electronic states. The JT splittings of X ̃ 2 E and B ̃ 2 E degenerate states are small, whereas it is quite large at three high frequencies in the C ̃ 2 E electronic states. The large energy separation of X ̃ 2 E from the other states and the non-zero PJT coupling of the B ̃ 2 E state with the close-lying A ̃ 2 A 1 state indicate the uncoupled nature of the X ̃ , A ̃ , and B ̃ vibronic bands of C4H3N. The intersection minima of B ̃ and C ̃ states with the D ̃ state nearly coincide with the energetic minimum of D ̃ state. Therefore, the PJT couplings among these states will lead to a strong vibronic interaction to shape the respective band structure. To completely understand the JT and PJT interactions in the photoionized spectrum of C4H3N, the vibronic coupling model Hamiltonian was constructed to perform nuclear dynamics studies for these electronic states. The vibrational progressions in each vibronic band are identified and compared with the available experimental data in the literature. The impacts of JT and PJT effects in the first five ionized states of cyanopropyne are investigated and discussed in detail. [ABSTRACT FROM AUTHOR]

The generation of cold molecules is an important topic in the field of cold atoms and molecules and has received relevant advanced research attention in ultracold chemistry, quantum computation, and quantum metrology. With a high atomic phase space density, optical dipole traps have been widely used to prepare, trap, and study cold molecules. In this work, Rb2 molecules were photoassociated in a magneto-optical trap to obtain a precise rovibrational spectrum, which provided accurate numerical references for the realization of multiple frequency photoassociation. By meeting the harsh requirements of photoassociation in optical dipole traps, the cold molecule photoassociation process was well explored, and different photoassociation resonances were simultaneously addressed in a single optical dipole trap. This method can be universally extended to simultaneously photoassociate cold molecules with different internal states or atomic species in a single optical dipole trap, thus advancing generous cold molecule studies such as cold molecule collision dynamics. [ABSTRACT FROM AUTHOR]

Second Harmonic Scattering (SHS) is a suitable technique to investigate the orientational correlations between molecules. This article explores the organization of different dye molecules adsorbed within the hydrophobic porosity of a hybrid organic–inorganic nanotube. Experimental polarization resolved SHS measurements highlight different orientational orders ranking from highly ordered and rigid organizations to disordered assemblies. Microscopic models of assemblies inside the pores are presented and discussed in the context of orientational correlation between the dye molecules. This work shows that the degree of order in the nanotube cavity follows the molecule's affinity within the porosity. [ABSTRACT FROM AUTHOR]

Investigating molecular excitation induced by collisions requires the prior determination of accurate analytical potential energy surfaces for the colliding partners. For carbon-chain molecules, such as cyanopolyynes, this has been a longstanding challenge, resulting in the absence of rate coefficients for HC5N, HC7N, HC9N, and others, induced by collisions with He. To overcome this bottleneck, we introduce a new approach: the Radial Angular Network with Gradual Expansion (RANGE). This method jointly connects the construction of ab initio interaction potentials with the determination of their analytical forms. We use the HC3N–He molecular complex as a reference to assess the reliability of our method, given that its analytical potential has been derived using various methods. Additionally, we apply the RANGE approach to construct the analytical representation of the interaction potential for HC5N–He and HC7N–He. The analysis of the analytical potentials reveals three systematic trends: (i) the anisotropy increases with the length of the carbon chain, (ii) the number of local minima correlates with the number of carbon atoms, and (iii) the shallowest local minimum is consistently located at ∼30 cm−1 below the dissociation limit of the complex. Using the time-independent quantum mechanical close-coupling formalism, we briefly estimate the propensity rules governing the excitation of HC3N, HC5N, and HC7N induced by collisions with He. Consequently, the three collisional systems exhibit the same propensity rule, favoring Δj = 2 transitions. [ABSTRACT FROM AUTHOR]

We apply the analytically solvable model of two electrons in two orbitals to diradical molecules, characterized by two unpaired electrons. The effect of doubly occupied and empty orbitals is taken into account by means of random phase approximation (RPA). We show that in the static limit, the direct RPA leads to the renormalization of the parameters of the two-orbital model. We test our model by comparing its predictions for singlet–triplet splitting with the results of several multi-reference methods for a set of thirteen molecules. We find that for this set, the static RPA results are close to those of the NEVPT2 method with two orbitals and two electrons in the active space. [ABSTRACT FROM AUTHOR]

We derive an analytic expression of the non-equilibrium Fermi's golden rule (NE-FGR) expression for a Holstein–Tavis–Cumming Hamiltonian, a universal model for many molecules collectively coupled to the optical cavity. These NE-FGR expressions capture the full-time-dependent behavior of the rate constant for transitions from polariton states to dark states. The rate is shown to be reduced to the well-known frequency domain-based equilibrium Fermi's golden rule (E-FGR) expression in the equilibrium and collective limit and is shown to retain the same scaling with the number of sites in non-equilibrium and non-collective cases. We use these NE-FGR to perform population dynamics with a time-non-local and time-local quantum master equation and obtain accurate population dynamics from the initially occupied upper or lower polariton states. Furthermore, NE-FGR significantly improves the accuracy of the population dynamics when starting from the lower polariton compared to the E-FGR theory, highlighting the importance of the non-Markovian behavior and the short-time transient behavior in the transition rate constant. [ABSTRACT FROM AUTHOR]

Intrinsically disordered proteins (IDPs) are prevalent participants in liquid–liquid phase separation due to their inherent potential for promoting multivalent binding. Understanding the underlying mechanisms of phase separation is challenging, as phase separation is a complex process, involving numerous molecules and various types of interactions. Here, we used a simplified coarse-grained model of IDPs to investigate the thermodynamic stability of the dense phase, conformational properties of IDPs, chain dynamics, and kinetic rates of forming condensates. We focused on the IDP system, in which the oppositely charged IDPs are maximally segregated, inherently possessing a high propensity for phase separation. By varying interaction strengths, salt concentrations, and temperatures, we observed that IDPs in the dense phase exhibited highly conserved conformational characteristics, which are more extended than those in the dilute phase. Although the chain motions and global conformational dynamics of IDPs in the condensates are slow due to the high viscosity, local chain flexibility at the short timescales is largely preserved with respect to that at the free state. Strikingly, we observed a non-monotonic relationship between interaction strengths and kinetic rates for forming condensates. As strong interactions of IDPs result in high stable condensates, our results suggest that the thermodynamics and kinetics of phase separation are decoupled and optimized by the speed-stability balance through underlying molecular interactions. Our findings contribute to the molecular-level understanding of phase separation and offer valuable insights into the developments of engineering strategies for precise regulation of biomolecular condensates. [ABSTRACT FROM AUTHOR]

In this work, the adsorption structure of deuterated water on the stepped platinum surface is studied under an ultra-high vacuum by using heterodyne-detected sum-frequency generation spectroscopy. On a pristine Pt(553), D2O molecules adsorbed at the step sites act as hydrogen bond (H-bond) donors to the adjacent terrace sites. This ensures the net D-down orientation at the terrace sites away from the steps. In particular, the pre-adsorption of oxygen atoms at the step sites significantly alters the D-down configuration. The oxygen pre-adsorption leads to a spontaneous dissociation of the post-adsorbed water molecules at the step to form hydroxyl (OD) species. Since the hydroxyl at the step acts as a strong H-bond acceptor, D2O at the terrace no longer maintains the D-down configuration and adopts flat-lying configurations, significantly reducing the number of D-down molecules at the terrace. Density-functional theoretical calculations support these pictures. This work demonstrates the critical role of steps in controlling the net orientation of the interfacial water and provides an important reference for future considerations of the reactions at electrochemical interfaces. [ABSTRACT FROM AUTHOR]

We present a comprehensive computational study describing the excited state dynamics and consequent photostability of amino-cyclohexenone (ACyO), the central template of mycosporine systems, widely recognized for their photoprotection of aquatic species. Photoexcitation to the first excited electronic state (S1, 1nπ*) of ACyO is considered an optically dark transition, while photoexcitation to the second excited electronic state (S21ππ*) is an optically bright 1ππ* transition and largely responsible for UV absorption properties of this molecule. We show that following initial photoexcitation to S2, ACyO relaxes via two competing deactivation mechanisms, each mediated by an S1/S0 conical intersection, which directs the excited state population to the electronic ground state (S0). Our ab initio computational results are supported with nonadiabatic dynamics simulation results, yielding an excited state lifetime of ∼280 fs for this system in vacuo. These results explain the inherent photostability of this core structure, commonplace in a wide range of microorganisms. [ABSTRACT FROM AUTHOR]

Motivated by experimental reports on chirality induced spin selectivity, we investigate a minimal model that allows us to calculate the charge and spin conductances through helical molecules analytically. The spin–orbit interaction is assumed to be non-vanishing on the molecule and negligible in the reservoirs (leads). The band structure of the molecule features four helical modes with spin-momentum locking that are analogous of edge-currents in the quantum spin Hall effect. While charge is conserved and therefore the charge current is independent of where it is measured—reservoirs or molecule—our detailed calculations reveal that the spin currents in the left and right leads are equal in magnitudes but with opposite signs (in linear response). We predict that transport currents flowing through helical molecules are accompanied by a spin accumulation in the contact region with the same magnetization direction for source and drain. Furthermore, we predict that the spin-conductance can be extracted directly from measuring the (quasi-static) spin accumulation—rather than the spin current itself, which is very challenging to obtain experimentally. [ABSTRACT FROM AUTHOR]

We report a quantum mechanical method for calculating the momentum distributions of constituent atoms of polyatomic molecules in rotational–vibrational eigenstates. Application of the present theory to triatomic molecules in the rovibrational ground state revealed that oscillatory changes appear on the proton momentum distribution in the nonlinear H2O molecule, while no such modulation is present in the case of an oxygen atom in the linear CO2 molecule. The atomic momentum distributions were analyzed in detail by means of a rigid rotor model, and it was found that the oscillation originates from the quantum-mechanical delocalization of the target atom with respect to the other atoms. [ABSTRACT FROM AUTHOR]

For floppy molecules showing internal rotations and/or large amplitude motions, curvilinear internal coordinates are known to be superior to rectilinear normal coordinates within vibrational structure calculations. Due to the myriad definitions of internal coordinates, automated and efficient potential energy surface generators necessitate a high degree of flexibility, supporting the properties arising from these coordinates. Within this work, an approach to deal with these challenges is presented, including key elements, such as the selection of appropriate fit functions, the exploitation of symmetry, the positioning of grid points, or elongation limits for different coordinates. These elements are tested for five definitions of curvilinear coordinates, with three of them being generated in an automated manner. Calculations for semi-rigid molecules, namely H2O, H2CO, CH2F2, and H2CNH, demonstrate the general functionality of the implemented algorithms. Additional calculations for the HOPO molecule highlight the benefits of these curvilinear coordinates for systems with large amplitude motions. This new implementation allowed us to compare the performance of these different coordinate systems with respect to the convergence of the underlying expansion of the potential energy surface and subsequent vibrational configuration interaction calculations. [ABSTRACT FROM AUTHOR]

RNA pseudoknots are RNA molecules with specialized three-dimensional structures that play important roles in various biological processes. To understand the functions and mechanisms of pseudoknots, it is essential to elucidate their structures and folding pathways. The most fundamental step in RNA folding is the opening and closing of a base pair. The effect of flexible loops on the base pair in pseudoknots remains unclear. In this work, we use molecular dynamics simulations and Markov state model to study the configurations, thermodynamic and kinetic of single base pair in pseudoknots. We find that the presence of the loop leads to a trap state. In addition, the rate-limiting step for the formation of base pair is the disruption of the trap state, rather than the open state to the closed state, which is quite different from the previous studies on non-pseudoknot RNA. For the thermodynamic parameters in pseudoknots, we find that the entropy difference upon opening the base pair between this simulation and the nearest-neighbor model results from the different entropy of different lengths of loop in solution. The thermodynamic parameters of the stack in pseudoknot are close to the nearest-neighbor parameters. The bases on the loop have different distribution patterns in different states, and the slow transition states of the loop are determined by the orientation of the bases. [ABSTRACT FROM AUTHOR]

Sizable glass formers feature numerous unique properties and potential applications, but many questions regarding their glass transition dynamics have not been resolved yet. Here, we have analyzed structural relaxation times measured as a function of temperature and pressure in combination with the equation of state obtained from pressure–volume–temperature measurements. Despite evidence from previous dielectric studies indicating a remarkable sensitivity of supercooled dynamics to compression, and contrary to intuition, our results demonstrated the proof for the almost equivalent importance of thermal energy and free volume fluctuations in controlling reorientation dynamics of sizable molecules. The found scaling exponent γ = 3.0 and Ev/Ep ratio of 0.6 were typical for glass-forming materials with relaxation dynamics determined by both effects with a minor advantage of thermal fluctuations involvement. It shows that the high values of key parameters characterizing the sensitivity of the glass transition dynamics to pressure changes, i.e., activation volume ΔV and dTg/dP, are not a valid premise for a remarkable contribution of volume to glass transition dynamics. [ABSTRACT FROM AUTHOR]

In this work, we study the energy flow through anharmonic coupling of vibrational modes after excitation of gas-phase fluorobenzene with a multi-THz pump. We show that to predict the efficiency of anharmonic energy transfer, simple models that only include the anharmonic coupling coefficients and motion of modes at their resonant frequency are not adequate. The full motion of each mode is needed, including the time while the mode is being driven by the pump pulse, because all the frequencies present in the multi-THz pump contribute to the excitation of the non-resonantly excited vibrational modes. Additionally, the model gives us the insight that modes with either A1 or B2 symmetry are more actively involved in anharmonic coupling because these modes have more symmetry-allowed energy transfer pathways. [ABSTRACT FROM AUTHOR]

Fully relativistic calculations, primarily at the 4-component coupled-cluster singles and doubles with perturbative triples [CCSD(T)] level of theory with the Dirac–Coulomb (DC) Hamiltonian, have been carried out for the superheavy (E125)F molecule using large Gaussian basis sets. The electronic ground state is determined to have an [Og]8s25g16f3 configuration on E125 with an Ω = 6 ground state and an 8p electron largely donated to F. A Mulliken population analysis indicates that the ground state is mainly ionic with a partial charge of +0.79 on E125 and a single sigma bond involving the F 2p and E125 8p spinors. The occupied g spinor is not involved in the bonding. With the largest basis set used in this work, the (0 K) dissociation energy was calculated at the DC-CCSD(T) level of theory to be 7.02 eV. Analogous calculations were also carried out for the E125 atom, both the neutral and its cation. The lowest energy electron configuration of E125+, [Og] 8 s 1 / 2 2 5 g 7 / 2 1 6 f 5 / 2 3 with a J = 6 ground state, was found to be similar to that in (E125)F, while the neutral E125 atom has an [Og] 8 s 1 / 2 2 5 g 7 / 2 1 6 f 5 / 2 2 7 d 3 / 2 1 8 p 1 / 2 1 ground state electron configuration with a J = 17/2 ground state. The ionization energy (IE) of E125 is reported for the first time and is calculated to be 4.70 eV at the DC-CCSD(T) level of theory. Non-relativistic calculations were also carried out on the E125 atom and the (E125)F molecule. The non-relativistic ground state of the E125 atom was calculated to have a 5g5 ground state with an IE of just 3.4 eV. The net effect of relativity on (E125)F is to stabilize its bonding. [ABSTRACT FROM AUTHOR]

A summary of the development of Monte Carlo statistical mechanics simulations for the computation of free energies of hydration of organic molecules is followed by presentation of results with the latest version of the optimized potentials for liquid simulations–all atom force field and the TIP4P water model. Scaling of the Lennard-Jones interactions between water, oxygen, and carbon atoms by a factor of 1.25 is found to improve the accuracy of free energies of hydration for 50 prototypical organic molecules from a mean unsigned error of 1.0–1.2 to 0.4 kcal/mol. [ABSTRACT FROM AUTHOR]

Electronic friction and Langevin dynamics is a popular mixed quantum–classical method for simulating the nonadiabatic dynamics of molecules interacting with metal surfaces, as it can be computationally more efficient than fully quantum approaches. In this work, we extend the theory of electronic friction within the hierarchical equations of motion formalism to models with a position-dependent metal–molecule coupling. We show that the addition of a position-dependent metal–molecule coupling adds new contributions to the electronic friction and other forces, which are highly relevant for many physical processes. Our expressions for the electronic forces within the Langevin equation are valid both in and out of equilibrium and for molecular models containing strong interactions. We demonstrate the approach by applying it to different models of interest. [ABSTRACT FROM AUTHOR]

The chemiresistive effect of an oxide significantly influences its electrical properties, which depend greatly on the interactions between the ambient gas molecules and the solid surface, including the gas adsorption and charge transfer still challenging to be clarified. In this work, we investigate the chemiresistive effect of the p-type delafossite CuScO2 microsheets by comparing their responses to various gaseous alcohols, which increase with an approximately linear relationship with the length of straight carbon chains from methanol to n-hexanol. A new mechanism is proposed to elucidate such a dramatic trend of observed chemiresistive change based on the first-principles calculations and test results. The increasing carbon chain length modulates the adsorption configuration and provides supplementary routes for electron transfer, which is assumed to account for the observed chemiresistive effect. This work may provide a novel perspective for the investigation and development of more advanced functional oxides for electrical applications. [ABSTRACT FROM AUTHOR]

We investigate the vertical electron attachment energy (VAE) of 1-methyl-4-nitroimidazole, a model radiosensitizer, employing quantum mechanics/molecular mechanics (QM/MM) and QM/polarized continuum (QM/PCM) solvation models. We considered the solvent-excluded surface (QM/PCM-SES) and Van der Waals (QM/PCM-VDW) cavities within the PCM framework, the electrostatic embedding QM/MM (EE-QM/MM) model, and the self-consistent sequential QM/MM polarizable electrostatic embedding (scPEE-S-QM/MM) model. Due to slow VAE convergence concerning the number of QM solvent molecules, full QM calculations prove inefficient. Ensemble averages in these calculations do not align with VAEs computed for the representative solute–solvent configuration. QM/MM and QM/PCM calculations show agreement with each other for sufficiently large QM regions, although the QM/PCM-VDW model exhibits artifacts linked to the cavity. QM/MM models demonstrate good agreement between ensemble averages and VAEs calculated with the representative configuration. Notably, the VAE computed with the scPEE-S-QM/MM model achieves faster convergence concerning the number of QM water molecules compared to the EE-QM/MM model, attributed to enhanced efficiency from MM charge polarization in the scPEE-S-QM/MM approach. This emphasizes the importance of QM/classical models with accurate solute–solvent and solvent–solvent mutual polarization for obtaining converged VAEs at a reasonable computational cost. The full-QM approach is very inefficient, while the microsolvation model is inaccurate. Computational savings in QM/MM models result from electrostatic embedding and the representative configuration, with the scPEE-S-QM/MM approach emerging as an efficient tool for describing bulk-solvated anions within the QM/MM framework. Its potential extends to improving transient anion state descriptions in biomolecules and radiosensitizers, especially given the frequent employment of microsolvation models. [ABSTRACT FROM AUTHOR]

We investigate an efficient spatial enantioseparation method of chiral molecules in cyclic three-level systems coupled with three optical fields using optically induced forces. When the overall phase differs by π between two enantiomers, significant variations in the magnitude and direction of the optically induced forces are observed. The manipulation of the center of mass of chiral molecules in optical fields can be achieved through the induced gauge force, primarily generated from the variations in the chirality-dependent scalar potentials created by the three inhomogeneous laser fields. By appropriately configuring the system, we can completely separate the slow spatial and fast inner dynamics, making instantaneous eigenstates of the inner Hamiltonian independent of the transverse profiles of the laser beams. Compared to previous methods, which required adiabatic conditions to be satisfied, the proposed method overcomes the limitations of the adiabatic approximation by utilizing a specific system configuration. This allows for increased flexibility in the transverse profiles of the laser beams and relaxes the constraints on the velocity of chiral molecules, leading to significantly greater spatial separations achievable across a broader range of parameters. [ABSTRACT FROM AUTHOR]

We present a four-component relativistic unitary coupled cluster method for atoms and molecules. We have used commutator-based non-perturbative approximation using the "Bernoulli expansion" to derive an approximation to the relativistic unitary coupled cluster method. The performance of the full quadratic unitary coupled-cluster singles and doubles method (qUCCSD), as well as a perturbative approximation variant (UCC3), has been reported for both energies and properties. It can be seen that both methods give results comparable to those of the standard relativistic coupled cluster method. The qUCCSD method shows better agreement with experimental results due to the better inclusion of relaxation effects. The relativistic UCC3 and qUCCSD methods can simulate the spin-forbidden transition with easy access to transition properties. A natural spinor-based scheme to reduce the computational cost of relativistic UCC3 and qUCCSD methods has been discussed. [ABSTRACT FROM AUTHOR]

The theory of polymer dynamics describes the intermediate scattering function for a polymer molecule in terms of relaxation modes defined by normal coordinates for the corresponding coarse-grained model. However, due to the difficulty of defining the normal coordinates for arbitrary polymer molecules, it is generally challenging to express the intermediate scattering function for a polymer molecule in terms of relaxation modes. To overcome this challenge, we propose a general method to calculate the intermediate scattering function for a polymer molecule on the basis of a relaxation mode analysis approach [Takano and Miyashita, J. Phys. Soc. Jpn. 64, 3688 (1995)]. In the proposed method, relaxation modes defined by eigenfunctions in a Markov process are evaluated on the basis of the simulation results for a polymer molecule and used to calculate the intermediate scattering function for that molecule. To demonstrate the effectiveness of the present method, we simulate the dynamics of a linear polymer molecule in a dilute solution and apply it to the calculation of the intermediate scattering function for the polymer molecule. The evaluation results regarding the relaxation modes reasonably describe the intermediate scattering function on the length scale of the radius of gyration of the polymer molecule. Accordingly, we examine the contributions of the pure relaxation and oscillatory relaxation processes to the entire intermediate scattering function. [ABSTRACT FROM AUTHOR]

Cholesterol (Chol) plays a crucial role in shaping the intricate physicochemical attributes of biomembranes, exerting a considerable influence on water molecules proximal to the membrane interface. In this study, we conducted molecular dynamics simulations on the bilayers of two lipid species, dipalmitoylphosphatidylcholine (DPPC) and palmitoyl sphingomyelin; they are distinct with respect to the structures of the hydrogen-bond (H-bond) acceptors. Our investigation focuses on the dynamic properties and H-bonds of water molecules in the lipid-membrane systems, with a particular emphasis on the influence of Chol at varying temperatures. Notably, in the gel phase at 303 K, the presence of Chol extends the lifetimes of H-bonds of the oxygen atoms acting as H-bond acceptors within DPPC with water molecules by a factor of 1.5–2.5. In the liquid-crystalline phase at 323 K, on the other hand, H-bonding dynamics with lipid membranes remain largely unaffected by Chol. This observed shift in H-bonding states serves as a crucial key to unraveling the subtle control mechanisms governing water dynamics in lipid-membrane systems. [ABSTRACT FROM AUTHOR]

Magnetometer cell wall coat molecules play an important role in preserving the lifetime of pumped alkali metal atoms for use in magnetometers that are capable of measuring very small magnetic fields. The goal of this study is to help rationalize the design of the cell coat molecules. Rubidium-87 is studied in terms of its interaction with three template cell coat molecules: ethane, ethene, and methyltrichlorosilane (MeTS). Ab initio electronic structure methods are applied to investigate the effect that the coat molecules have on the 2S ground state and 2P excited state of 87Rb. We find that, from the ab initio results, the three template molecules have differing effects, with MeTS having the largest effect on the ground state and ethane or ethene having the largest effect on the non-degenerate excited states. [ABSTRACT FROM AUTHOR]

We show that optimal control of the electron dynamics is able to prepare molecular ground states, within chemical accuracy, with evolution times approaching the bounds imposed by quantum mechanics. We propose a specific parameterization of the molecular evolution only in terms of interaction already present in the molecular Hamiltonian. Thus, the proposed method solely utilizes quantum simulation routines, retaining their favorable scalings. Due to the intimate relationships between variational quantum algorithms and optimal control, we compare, when possible, our results with state-of-the-art methods in the literature. We found that the number of parameters needed to reach chemical accuracy and algorithmic scaling is in line with compact adaptive strategies to build variational Ansätze. The algorithm, which is also suitable for quantum simulators, is implemented by emulating a digital quantum processor (up to 16 qubits) and tested on different molecules and geometries spanning different degrees of electron correlation. [ABSTRACT FROM AUTHOR]

The control of movement and orientation of gas-phase molecules has become the focus of many research areas in molecular physics. Here, ND3 molecules are polarized in a segmented, curved electrostatic guide and adiabatically aligned inside a rotatable mass spectrometer (MS). Alignment is probed by photoionization using a linearly polarized laser. Rotation of the polarization at fixed MS orientation has the same effect as the rotation of the MS at fixed polarization, proving that the molecular alignment adiabatically follows the MS axis. Polarization-dependent ion signals reveal state-specific populations and allow for a quantification of the aligned sample in the space-fixed reference frame. [ABSTRACT FROM AUTHOR]

Dimensions and molecular structures play pivotal roles in the principle of heat conduction. The dimensional characteristics of a solution within nanoscale systems depend on the degrees of confinement. However, the influence of such variations on heat transfer remains inadequately understood. Here, we perform quasi-one-dimensional non-equilibrium molecular dynamics simulations to calculate the thermal conductivity of water molecules confined in carbon nanotubes. The structure of water molecules is determined depending on the nanotube radius, forming a single-file, a single-layer, and a double-layer structure, corresponding to an increasing radius order. We reveal that the thermal conductivity of liquid water has a sublinear dependency on nanotube length exclusively when water molecules form a single file. A stronger confinement leads to behavioral and structural characteristics closely resembling a one-dimensional nature. Moreover, single-layer-structured water molecules exhibit enhanced thermal conductivity. We elucidate that this is due to the increase in the local water density and the absence of transitions to another layer, which typically occurs in systems with double-layer water structures within relatively large radius nanotubes. [ABSTRACT FROM AUTHOR]

A novel mixed quantum–classical approach to simulating nonadiabatic dynamics of molecules at metal surfaces is presented. The method combines the numerically exact hierarchical equations of motion approach for the quantum electronic degrees of freedom with Langevin dynamics for the classical degrees of freedom, namely, low-frequency vibrational modes within the molecule. The approach extends previous mixed quantum–classical methods based on Langevin equations to models containing strong electron–electron or quantum electronic–vibrational interactions, while maintaining a nonperturbative and non-Markovian treatment of the molecule–metal coupling. To demonstrate the approach, nonequilibrium transport observables are calculated for a molecular nanojunction containing strong interactions. [ABSTRACT FROM AUTHOR]

We have determined the bond dissociation energy of O2 by measuring fully state-to-state resolved threshold fragment yield spectra in the XUV energy region, O 2 X 3 Σ g − , N ″ , J ″ → O ( P J 3 ) + O ( S 1 o 3 ) / O ( S 2 o 5 ). Our results have yielded a bond dissociation energy value of 41 269.19 ± 0.10 cm−1, which is consistent with previous measurements but exhibits a significantly lower uncertainty, approximately five times smaller. It is noteworthy that this study is the first to simultaneously achieve fine structure state resolution for the parent O2 molecule and spin–orbit state resolution for the O(3PJ) fragments in the measurement of O2 bond dissociation energy. As a result, our findings have established a solid foundation for the obtained data. [ABSTRACT FROM AUTHOR]

Molecular chirality has traditionally been viewed as a binary property where a molecule is classified as either chiral or achiral, yet in recent decades, mathematical methods for quantifying chirality have been explored. Here, we use toy molecular systems to systematically compare the performance of two state-of-the-art chirality measures: (1) the Continuous Chirality Measure (CCM) and (2) the Chirality Characteristic (χ). We find that both methods exhibit qualitatively similar behavior when applied to simple molecular systems such as a four-site molecule or the polymer double-helix, but we show that the CCM may be more suitable for evaluating the chirality of arbitrary molecules or abstract structures such as normal vibrational modes. We discuss a range of considerations for applying these methods to molecular systems in general, and we provide links to user-friendly codes for both methods. We aim for this paper to serve as a concise resource for scientists attempting to familiarize themselves with these chirality measures or attempting to implement chirality measures in their own work. [ABSTRACT FROM AUTHOR]

In the rich ice polymorphism landscape, ice XVII, metastable at ambient pressure and at temperatures below 130 K, is surely one of the most interesting from both fundamental and technological perspectives due to its porosity, i.e., its capability to repeatedly absorb and desorb molecular hydrogen by dosing the gas at pressures even below the ambient one. Here, owing to this exceptional key feature, we investigate the roto-vibrational dynamics of the H2 molecules trapped in the fully deuterated ice XVII structure. Making use of the high-resolution and brilliance of the TOSCA neutron vibrational spectrometer, combined with high-resolution Raman data, we are able to efficiently distinguish the center-of-mass translational bands from the rotational ones and to study them as a function of the guest filling of the ice structure, unraveling a peculiar behavior for the confined particle in a low-dimensional system. Moreover, we also report the study of the microscopic dynamics of confined nitrogen and oxygen, which are the most abundant molecular species in the atmosphere and are of paramount interest for technological applications. Finally, we show that the ice XVII porosity is a unique feature, especially in the low pressure regime, within the emptied-hydrate phases discovered to date. [ABSTRACT FROM AUTHOR]

A high-resolution absorption spectrum of the S1–S0 transition of free-base phthalocyanine was observed and analyzed with improved reliability. The spectrum, with a partially resolved rotational structure, was obtained by using the buffer-gas cooling technique and a single-mode tunable laser. Our new analysis reveals that the S 1 ← S 0 0 0 0 band belongs to the a-type transition, where the electronic transition moment aligns parallel to the NH–HN direction, allowing the assignment of the S1 state to 1B3u. These results agree with a prior study using supersonic expansion and are well supported by theoretical calculations. Interestingly, the rotational constant B in the S1 state, which is often smaller than that in the ground state for typical molecules, was found to be slightly larger than that in the S01Ag state. This suggests a change in the character of π bonds with the electronic excitation. [ABSTRACT FROM AUTHOR]

The notions of ionicity and covalency of chemical bonds, effective atomic charges, and decomposition of the cohesive energy into ionic and covalent terms are fundamental yet elusive. For example, different approaches give different values of atomic charges. Pursuing the goal of formulating a universal approach based on firm physical grounds (first-principles or non-empirical), we develop a formalism based on Wannier functions with atomic orbital symmetry and capable of defining these notions and giving numerically robust results that are in excellent agreement with traditional chemical thinking. Unexpectedly, in diamond-like boron phosphide (BP), we find charges of +0.68 on phosphorus and −0.68 on boron atoms, and this anomaly is explained by the Zintl–Klemm nature of this compound. We present a simple model that includes energies of the highest occupied cationic and lowest unoccupied anionic atomic orbitals, coordination numbers, and strength of interatomic orbital overlap. This model captures the essential physics of bonding and accurately reproduces all our results, including anomalous BP. [ABSTRACT FROM AUTHOR]

Resonant electron capture by aza and diaza derivatives of phenanthrene (7,8-benzoquinoline and 1,10-phenanthroline) and anthracene (acridine and phenazine) at incident free electron energies (Ee) in the range of 0–15 eV was studied. All compounds except 7,8-benzoquinoline form long-lived molecular ions (M−) at thermal electron energies (Ee ∼ 0 eV). Acridine and phenazine also form such ions at epithermal electron energies up to Ee = 1.5–2.5 eV. The lifetimes (τa) of M− with respect to electron autodetachment are proportional to the extent of aza-substitution and increase on going from molecules with bent geometry of the fused rings (azaphenanthrenes) to linear isomers (azaanthracenes). These regularities are due to an increase in the adiabatic electron affinities (EAa) of the molecules. The EAa values of the molecules under study were comprehensively assessed based on a comparative analysis of the measured τa values using the Rice–Ramsperger–Kassel–Marcus theory, the electronic structure analysis using the molecular orbital approach, as well as the density functional calculations of the total energy differences between the molecules and anions. The only fragmentation channel of M− ions from the compounds studied is abstraction of hydrogen atoms. When studying [M–H]− ions, electron autodetachment processes were observed, the τa values were measured, and the appearance energies were determined. A comparative analysis of the gas-phase acidity of the molecules and the EAa values of the [M–H]· radicals revealed their proportionality to the EAa values of the parent molecules. [ABSTRACT FROM AUTHOR]

This article is a response to criticisms raised by Hirano et al. regarding research on sum frequency vibrational spectroscopy (SFVS) of benzene at interfaces. The authors address five points raised in the comment, including the local field effect, symmetry breaking of benzene, and field variations at interfaces. They argue that Hirano et al. underestimated the interfacial EQ SFVS and did not calculate the absolute amplitude of SFVS, which undermines the credibility of their research. The authors' computational SFVS aligns well with experimental results in terms of vibrational frequencies, amplitude, and line shapes, but there are discrepancies in vibrational frequencies and spectral line shapes between their results and the experimental data. [Extracted from the article]

Using molecular dynamics simulations, we show that a molecule of moderately active antifreeze protein (type III AFP, QAE HPLC-12 isoform) is able to interact with ice in an indirect manner. This interaction occurs between the ice binding site (IBS) of the AFP III molecule and the surface of ice, and it is mediated by liquid water, which separates these surfaces. As a result, the AFP III molecule positions itself at a specific orientation and distance relative to the surface of ice, which enables the effective binding (via hydrogen bonds) of the molecule with the nascent ice surface. Our results show that the final adsorption of the AFP III molecule on the surface of ice is not achieved by chaotic diffusion movements, but it is preceded by a remote, water-mediated interaction between the IBS and the surface of ice. The key factor that determines the existence of this interaction is the ability of water molecules to spontaneously form large, high-volume aggregates that can be anchored to both the IBS of the AFP molecule and the surface of ice. The results presented in this work for AFP III are in full agreement with the ones obtained by us previously for hyperactive CfAFP, which indicates that the mechanism of the remote interaction of these molecules with ice remains unchanged despite significant differences in the molecular structure of their ice binding sites. For that reason, we can expect that also other types of AFPs interact with the ice surface according to an analogous mechanism. [ABSTRACT FROM AUTHOR]

We develop and demonstrate how to use the Graphical Unitary Group Approach (GUGA)-based MRCISD with Core–Valence Separation (CVS) approximation to compute the core-excited states. First, perform a normal Self-Consistent-Field (SCF) or valence MCSCF calculation to optimize the molecular orbitals. Second, rotate the optimized target core orbitals and append to the active space, form an extended CVS active space, and perform a CVS-MCSCF calculation for core-excited states. Finally, construct the CVS-MRCISD expansion space and perform a CVS-MRCISD calculation to optimize the CI coefficients based on the variational method. The CVS approximation with GUGA-based methods can be implemented by flexible truncation of the Distinct Row Table. Eliminating the valence-excited configurations from the CVS-MRCISD expansion space can prevent variational collapse in the Davidson iteration diagonalization. The accuracy of the CVS-MRCISD scheme was investigated for excitation energies and compared with that of the CVS-MCSCF and CVS-CASPT2 methods using the same active space. The results show that CVS-MRCISD is capable of reproducing well-matched vertical core excitation energies that are consistent with experiments by combining large basis sets and a rational reference space. The calculation results also highlight the fact that the dynamic correlation between electrons makes an undeniable contribution in core-excited states. [ABSTRACT FROM AUTHOR]

We investigate the effect of a cavity on nonlinear two-photon transitions of a molecular system and we analyze how such an effect depends on the cavity quality factor, the field enhancement, and the possibility of dephasing. We find that the molecular response to strong light fields in a cavity with a variable quality factor can be understood as arising from a balance between (i) the ability of the cavity to enhance the field of an external probe and promote multiphoton transitions more easily and (ii) the fact that the strict selection rules on multiphoton transitions in a cavity support only one resonant frequency within the excitation range. Although our simulations use a classical level description of the radiation field (i.e., we solve Maxwell–Bloch or Maxwell–Liouville equations within the Ehrenfest approximation for the field–molecule interaction), based on experience with this level of approximation in the past studies of plasmonic and polaritonic systems, we believe that our results are valid over a wide range of external probing. [ABSTRACT FROM AUTHOR]

We propose a Brownian ratchet for the unidirectional transport of stimuli-responsive molecules confined in a series of asymmetric geometries. It relies on repetitive cycles of aggregation and dispersion, which cause significant changes in molecular distribution within the confining geometry and enable the Brownian motion of the molecules to be ratcheted in a specific direction. To demonstrate the feasibility of the proposed Brownian ratchet, we conducted Brownian dynamics simulations where stimuli-responsive molecules were repeatedly aggregated and dispersed in a series of truncated conical tubes by altering intermolecular interactions. These simulations demonstrated the unidirectional transport of the molecules, indicating the efficacy of the proposed Brownian ratchet. Furthermore, we found that it becomes more effective with higher concentrations of molecules. This study suggests that, through the deliberate control of molecular assembly and disassembly by stimuli-responsive intermolecular interactions, it is possible to achieve directional and controlled molecular transport in various nanoscale applications. [ABSTRACT FROM AUTHOR]

The interplay of molecules gives rise to collective phenomena absent in a single molecule. Many examples of collective phenomena have been reported as their knowledge is essential for understanding the behavior of matter. Here, we consider molecules sufficiently separated from each other to not form chemical bonds. If these molecules are excited, e.g., by a weak laser, can they concertedly relax by emitting a single high-energy photon possessing the total energy of all the relaxing molecules? We show that this concerted emission process is indeed possible. We estimate its probability and analyze its dependence on molecular properties, intermolecular distances, and relative orientations of the molecules. A numerical example on two pyridine molecules is given. The concerted emission found is a fundamental process expected to be operative in gas phase and clusters. Its true relevance lies in its intimate relationship to concerted emission of virtual photons and thus to collective energy transfer ionizing neighboring systems. The estimated rates and examples discussed of this collective intermolecular Coulombic decay shed much light on recent puzzling experiments. [ABSTRACT FROM AUTHOR]

The Effective Relativistic Coupling by Asymptotic Representation (ERCAR) approach is a method to generate fully coupled diabatic potential energy surfaces (PESs) including relativistic effects, especially spin–orbit coupling. The spin–orbit coupling of a full molecule is determined only by the atomic states of selected relativistically treated atoms. The full molecular coupling effect is obtained by a diabatization with respect to asymptotic states, resulting in the correct geometry dependence of the spin–orbit effect. The ERCAR approach has been developed over the last decade and initially only for molecules with a single relativistic atom. This work presents its extension to molecules with more than a single relativistic atom using the iodine molecule as a proof-of-principle example. The theory for the general multiple atomic ERCAR approach is given. In this case, the diabatic basis is defined at the asymptote where all relativistic atoms are separated from the remaining molecular fragment. The effective spin–orbit operator is then a sum of spin–orbit operators acting on isolated relativistic atoms. PESs for the iodine molecule are developed within the new approach and it is shown that the resulting fine structure states are in good agreement with spin–orbit ab initio calculations. [ABSTRACT FROM AUTHOR]

The effect of a locally excited state on charge transfer symmetry breaking (SBCT) in excited quadrupolar molecules in solutions has been studied. The interaction of a locally excited state and two zwitterionic states is found to either increase or decrease the degree of SBCT depending on the molecular parameters. A strategy on how to adjust the molecular parameters to control the extent of SBCT is presented. The influence of level degeneracy on SBCT is identified and discussed in detail. The level degeneracy is shown to lead to the existence of a hidden dipole moment in excited quadrupolar molecules. Its manifestations in SBCT are analyzed. The main conclusions are consistent with the available experimental data. [ABSTRACT FROM AUTHOR]

Semilocal density-functional approximations (DFAs), including the state-of-the-art SCAN functional, are plagued by the self-interaction error (SIE). While this error is explicitly defined only for one-electron systems, it has inspired the self-interaction correction method proposed by Perdew and Zunger (PZ-SIC), which has shown promise in mitigating the many-electron SIE. However, the PZ-SIC method is known for its significant numerical instability. In this study, we introduce a novel constraint that facilitates self-consistent localization of the SIC orbitals in the spirit of Edmiston–Ruedenberg orbitals [Rev. Mod. Phys. 35, 457 (1963)]. Our practical implementation within the all-electron numeric atom-centered orbitals code FHI-aims guarantees efficient and stable convergence of the self-consistent PZ-SIC equations for both molecules and solids. We further demonstrate that our PZ-SIC approach effectively mitigates the SIE in the meta-generalized gradient approximation SCAN functional, significantly improving the accuracy for ionization potentials, charge-transfer energies, and bandgaps for a diverse selection of molecules and solids. However, our PZ-SIC method does have its limitations. It cannot improve the already accurate SCAN results for properties such as cohesive energies, lattice constants, and bulk modulus in our test sets. This highlights the need for new-generation DFAs with more comprehensive applicability. [ABSTRACT FROM AUTHOR]

Core–shell gold–silver cuboidal nanoparticles were produced, with either concave or straight facets. Their incubation with a low concentration of chiral l-glutathione (GSH) biomolecules was found to produce near UV plasmonic extinction and induced circular dichroism (CD) peaks. The effect is sensitive to the silver shell thickness. The GSH molecules were found to cause redistribution of silver in the shell, removing silver atoms from edges/corners and re-depositing them at the nanocuboid facets, probably through some redox and complexation processes between the silver and thiol group of the GSH. Other thiolated chiral biomolecules (and drug molecules) did not show this effect. The emerging near UV surface plasmon resonance is a silver slab resonance, which might also possess some multipolar resonance nature. The concave-shaped nanocuboids exhibited stronger induced plasmonic CD relative to the nanocuboids with straight facets. [ABSTRACT FROM AUTHOR]

This study investigated the dissociation after the Xe 4d Auger decay of weak-bonding XeF2 molecules by multielectron–ion coincidence spectroscopy using a magnetic bottle electron spectrometer. Fragmentations from the XeF22+ states were clarified in the Auger spectra coincident with individual ion species. It was observed that the two-hole population led by the Auger decay was not directly inherited during the fragmentation of XeF22+. Furthermore, the dissociations of XeF23+ states produced by the double Auger decay were investigated. [ABSTRACT FROM AUTHOR]