Search

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]

We have analyzed the effects of the spreads in experimental parameters on the reliability of speeds and angular distributions extracted from a generic surface-scattering experiment based on planar laser-induced fluorescence detection. The numerical model assumes a pulsed beam of projectile molecules is directed at a surface. The spatial distribution of the scattered products is detected by imaging the laser-induced fluorescence excited by a thin, pulsed sheet of laser light. Monte Carlo sampling is used to select from realistic distributions of the experimental parameters. The key parameter is found to be the molecular-beam diameter, expressed as a ratio to the measurement distance from the point of impact. Measured angular distributions are negligibly distorted when this ratio is <∼10%. Measured most-probable speeds are more tolerant, being undistorted when it is <∼20%. In contrast, the spread of speeds or of corresponding arrival times in the incident molecular beam has only very minor systematic effects. The thickness of the laser sheet is also unimportant within realistic practical limits. These conclusions are broadly applicable to experiments of this general type. In addition, we have analyzed the specific set of parameters designed to match the experiments on OH scattering from a liquid perfluoropolyether (PFPE) surface in the Paper I [Roman et al., J. Chem. Phys. 158, 244704 (2023)]. This reveals that the detailed form of the molecular-beam profile is important, particularly on apparent angular distributions, for geometric reasons that we explain. Empirical factors have been derived to correct for these effects. [ABSTRACT FROM AUTHOR]

This is an I easier approach to implement into computer code i , because these equations [Eq. (95) of the paper] are simpler than the corresponding EOMs for { I i SB I n i sb , I i SB I n i sb }. CLARIFICATION ON THE NUMERICAL ALGORITHM USED TO PROPAGATE DYNAMICS We want to clarify the numerical algorithm we used to propagate the EOMs and generate all numerical results presented in the paper. In the above expressions, to compute HT ht , we use Eq. (C2) of the paper. [Extracted from the article]

A preceding paper [M. Dhar, J. A. Dickinson, and M. A. Berg, J. Chem. Phys. 159, 054110 (2023)] shows how to remove additive noise from an experimental time series, allowing both the equilibrium distribution of the system and its Green's function to be recovered. The approach is based on nonlinear-correlation functions and is fully nonparametric: no initial model of the system or of the noise is needed. However, single-molecule spectroscopy often produces time series with either photon or photon-counting noise. Unlike additive noise, photon noise is signal-size correlated and quantized. Photon counting adds the potential for bias. This paper extends noise-corrected-correlation methods to these cases and tests them on synthetic datasets. Neither signal-size correlation nor quantization is a significant complication. Analysis of the sampling error yields guidelines for the data quality needed to recover the properties of a system with a given complexity. We show that bias in photon-counting data can be corrected, even at the high count rates needed to optimize the time resolution. Using all these results, we discuss the factors that limit the time resolution of single-molecule spectroscopy and the conditions that would be needed to push measurements into the submicrosecond region. [ABSTRACT FROM AUTHOR]

This paper focuses on phase and aggregation behavior for linear chains composed of blocks of hydrophilic and hydrophobic segments. Phase and conformational transitions of patterned chains are relevant for understanding liquid–liquid separation of biomolecular condensates, which play a prominent role in cellular biophysics and for surfactant and polymer applications. Previous studies of simple models for multiblock chains have shown that, depending on the sequence pattern and chain length, such systems can fall into one of two categories: displaying either phase separation or aggregation into finite-size clusters. The key new result of this paper is that both formation of finite-size aggregates and phase separation can be observed for certain chain architectures at appropriate conditions of temperature and concentration. For such systems, a bulk dense liquid condenses from a dilute phase that already contains multi-chain finite-size aggregates. The computational approach used in this study involves several distinct steps using histogram-reweighting grand canonical Monte Carlo simulations, which are described in some level of detail. [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]

This article is a response to a comment made by Professor Liu regarding a previously published paper. The response addresses questions raised by Professor Liu and focuses on three key aspects of the validity of the paper. It explains the use of different contractions in the construction of the QED Hamiltonians and clarifies the commutation relations used in the calculations. The article also highlights that the formulation of the molecular orbital method described in the paper is independent of the ordering of the operators and can derive expressions for various perturbation theories. The response concludes by stating that alternative criteria for the QED Hamiltonians are not ruled out and could be explored in future research. [Extracted from the article]

This document is an erratum for a previous paper titled "Modified Fermi's golden rule rate expressions." The erratum corrects an error in equation (41) by adding a missing factor of 1/2. The correction does not change any major conclusions or implications of the paper, but it does require the correction of figures 1-4, which show results of model calculations. The corrected figures and their captions are provided in the erratum. [Extracted from the article]

This work follows a companion article, which will be referred to as Paper I [Campeggio et al., J. Chem. Phys. 158, 244104 (2023)] in which a quantum-stochastic Liouville equation for the description of the quantum–classical dynamics of a molecule in a dissipative bath has been formulated in curvilinear internal coordinates. In such an approach, the coordinates of the system are separated into three subsets: the quantum coordinates, the classical relevant nuclear degrees of freedom, and the classical irrelevant (bath) coordinates. The equation has been derived in natural internal coordinates, which are bond lengths, bond angles, and dihedral angles. The resulting equation needs to be parameterized. In particular, one needs to compute the potential energy surfaces, the friction tensor, and the rate constants for the nonradiative jumps among the quantum states. While standard methods exist for the calculation of energy and dissipative properties, an efficient evaluation of the transition rates needs to be developed. In this paper, an approximated treatment is introduced, which leads to a simple explicit formula with a single adjustable parameter. Such an approximated expression is compared with the exact calculation of transition rates obtained via molecular dynamics simulations. To make such a comparison possible, a simple sandbox system has been used, with two quantum states and a single internal coordinate (together with its conjugate momentum). Results show that the adjustable parameter, which is an effective decoherence time, can be parameterized from the effective relaxation times of the autocorrelation functions of the conjugated momenta of the relevant nuclear coordinates. [ABSTRACT FROM AUTHOR]

Solvation dynamics in ionic liquids show features that are often associated with supercooled liquids, including "stretched" nonexponential relaxation. To better understand the mechanism behind the stretching, the nonlinear mode-correlation methods proposed in Paper I [S. R. Hodge and M. A. Berg, J. Chem. Phys. 155, 024122 (2021)] are applied to a simulation of a prototypical ionic liquid. A full Green's function is recovered. In addition, specific tests for non-Gaussian dynamics are made. No deviations from Gaussian dynamics are found. This finding is incompatible with rate heterogeneity as a cause of the nonexponential relaxation and appears to be in conflict with an earlier multidimensional analysis of the same data. Although this conflict is not resolved here, this work does demonstrate the practicality of mode-correlation analysis in the face of finite datasets and calculations. [ABSTRACT FROM AUTHOR]

The study and modeling of water exchange in complex media using different applications of diffusion and relaxation magnetic resonance (MR) have been of interest in recent years. Most models attempt to describe this process using a first order kinetics expression, which is appropriate to describe chemical exchange; however, it may not be suitable to describe diffusion-driven exchange since it has no direct relationship to diffusion dynamics of water molecules. In this paper, these limitations are addressed through a more general exchange expression that does consider such important properties. This exchange fraction expression features a multi-exponential recovery at short times and a mono-exponential decay at long times, both of which are not captured by the first order kinetics expression. Furthermore, simplified exchange expressions containing partial information of the analyzed system's diffusion and relaxation processes and geometry are proposed, which can potentially be employed in already established estimation protocols. Finally, exchange fractions estimated from simulated MR data and derived here were compared, showing qualitative similarities but quantitative differences, suggesting that the features of the derived exchange fraction in this paper can be partially recovered by employing an existing estimation framework. [ABSTRACT FROM AUTHOR]

Biased and accelerated molecular simulations (BAMS) are widely used tools to observe relevant molecular phenomena occurring on time scales inaccessible to standard molecular dynamics, but evaluation of the physical time scales involved in the processes is not directly possible from them. For this reason, the problem of recovering dynamics from such kinds of simulations is the object of very active research due to the relevant theoretical and practical implications of dynamics on the properties of both natural and synthetic molecular systems. In a recent paper [A. Rapallo et al., J. Comput. Chem. 42, 586–599 (2021)], it has been shown how the coupling of BAMS (which destroys the dynamics but allows to calculate average properties) with Extended Diffusion Theory (EDT) (which requires input appropriate equilibrium averages calculated over the BAMS trajectories) allows to effectively use the Smoluchowski equation to calculate the orientational time correlation function of the head–tail unit vector defined over a peptide in water solution. Orientational relaxation of this vector is the result of the coupling of internal molecular motions with overall molecular rotation, and it was very well described by correlation functions expressed in terms of weighted sums of suitable time-exponentially decaying functions, in agreement with a Brownian diffusive regime. However, situations occur where exponentially decaying functions are no longer appropriate to capture the actual dynamical behavior, which exhibits persistent long time correlations, compatible with the so called subdiffusive regimes. In this paper, a generalization of EDT will be given, exploiting a fractional Smoluchowski equation (FEDT) to capture the non-exponential character observed in the relaxation of intramolecular distances and molecular radius of gyration, whose dynamics depend on internal molecular motions only. The calculation methods, proper to EDT, are adapted to implement the generalization of the theory, and the resulting algorithm confirms FEDT as a tool of practical value in recovering dynamics from BAMS, to be used in general situations, involving both regular and anomalous diffusion regimes. [ABSTRACT FROM AUTHOR]

Despite the importance of rhodium complexes in catalysis, and the favorable 100% natural abundance of the spin-1/2 103Rh nucleus, there are few reports of 103Rh nuclear magnetic resonance (NMR) parameters in the literature. In part, this is the consequence of the very low gyromagnetic ratio of 103Rh and its dismal NMR sensitivity. In a previous paper [Harbor-Collins et al., J. Chem. Phys. 159, 104 307 (2023)], we demonstrated an NMR methodology for 1H-enhanced 103Rh NMR and demonstrated an application to the 103Rh NMR of the dirhodium formate paddlewheel complex. In this paper, we employ selective 18O labeling to break the magnetic equivalence of the 103Rh spin pair of dirhodium formate. This allows the estimation of the 103Rh–103Rh spin–spin coupling and provides access to the 103Rh singlet state. We present the first measurement of a 18O-induced 103Rh secondary isotope shift as well as the first instance of singlet order generated in a 103Rh spin pair. The field-dependence of 103Rh singlet relaxation is measured by field-cycling NMR experiments. [ABSTRACT FROM AUTHOR]

In this paper, we explore the molecular basis of combining photodynamic therapy (PDT), a light-triggered targeted anticancer therapy, with the traditional chemotherapeutic properties of the well-known cytotoxic agent gemcitabine. A photosensitizer prerequisite is significant absorption of biocompatible light in the visible/near IR range, ideally between 600 and 1000 nm. We use highly accurate multiconfigurational CASSCF/MS-CASPT2/MM and TD-DFT methodologies to determine the absorption properties of a series of gemcitabine derivatives with the goal of red-shifting the UV absorption band toward the visible region and facilitating triplet state population. The choice of the substitutions and, thus, the rational design is based on important biochemical criteria and on derivatives whose synthesis is reported in the literature. The modifications tackled in this paper consist of: (i) substitution of the oxygen atom at O2 position with heavier atoms (O → S and O → Se) to red shift the absorption band and increase the spin–orbit coupling, (ii) addition of a lipophilic chain at the N7 position to enhance transport into cancer cells and slow down gemcitabine metabolism, and (iii) attachment of aromatic systems at C5 position to enhance red shift further. Results indicate that the combination of these three chemical modifications markedly shifts the absorption spectrum toward the 500 nm region and beyond and drastically increases spin–orbit coupling values, two key PDT requirements. The obtained theoretical predictions encourage biological studies to further develop this anticancer approach. [ABSTRACT FROM AUTHOR]

Comparing a discrete molecular spectrum to a continuous molecular spectrum in a quantitative manner is a challenging problem, for example, when attempting to fit a theoretical stick spectrum to a continuous spectrum. In this paper, the use of computational optimal transport is investigated for such a problem. In the optimal transport literature, the comparison of a discrete and a continuous spectrum is referred to as semi-discrete optimal transport and is a situation where a metric such as least-squares may be difficult to define except under special conditions. The merits of an optimal transport approach for this problem are investigated using the transport distance defined for the semi-discrete case. A tutorial on semi-discrete optimal transport for molecular spectra is included in this paper, and several well-chosen synthetic spectra are investigated to demonstrate the utility of computational optimal transport for the semi-discrete case. Among several types of investigations, we include calculations showing how the frequency resolution of the continuous spectrum affects the transport distance between a discrete and a continuous spectrum. We also use the transport distance to measure the distance between a continuous experimental electronic absorption spectrum of SO2 and a theoretical stick spectrum for the same system. The comparison of the theoretical and experimental SO2 spectra also allows us to suggest a theoretical value for the band origin that is closer to the observed band origin than previous theoretical values. [ABSTRACT FROM AUTHOR]

Lithium–ion batteries are key elements of electric vehicles and energy storage systems, and their accurate State of Charge (SOC) estimation is momentous for battery energy management, safe operation, and extended service life. In this paper, the Multi-Kernel Relevance Vector Machine (MKRVM) and Whale Optimization Algorithm (WOA) are used to estimate the SOC of lithium–ion batteries under different operating conditions. In order to better learn and estimate the battery SOC, MKRVM is used to establish a model to estimate lithium–ion battery SOC. WOA is used to automatically adjust and optimize weights and kernel parameters of MKRVM to improve estimation accuracy. The proposed model is validated with three lithium–ion batteries under different operating conditions. In contrast to other optimization algorithms, WOA has a better optimization effect and can estimate the SOC more accurately. In contrast to the single kernel function, the proposed multi-kernel function greatly improves the precision of the SOC estimation model. In contrast to the traditional method, the WOA-MKRVM has a higher precision of SOC estimation. [ABSTRACT FROM AUTHOR]

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

Jiang and co-workers use high resolution STM to investigate the reaction and self-assembly of (3,6-dibromo-9,10-phenanthrenequinone, or DBPQ) molecules on Ag (100) and Ag (110) surfaces in order to understand the mechanism of bottom-up assembly on surfaces.[31] They show that, through the inclusion of multiple functional groups within a precursor molecule, it becomes possible to fabricate new low-dimensional materials with unique chemical, physical, and electronic properties. Herbst and Fransson consider the core-valence separation approximation that is often used in the calculation of core-level spectra.[5] They show how to quantify the errors in this approximation, thereby opening the door to error-quantified predictions relevant to x-ray spectroscopy. 153(16), 164108 (2020).10.1063/5.0019557 5 M. F. Herbst and T. Fransson, "Quantifying the error of the core-valence separation approximation", J. Chem. Phys. Zhu and co-workers tackle this problem for a model system containing a 2D semiconductor heterojunction and show convincingly the efficient hot electron transfer from photoexcited MoTe SB 2 sb to WS SB 2 sb .[30] This finding provides important insight into the competition between hot electron cooling and transfer at 2D semiconductor interfaces and suggests an intriguing possibility for the exploration of hot electron devices. [Extracted from the article]

Due to the ubiquity and importance of water, water dimer has been intensively studied. Computing the (ro-)vibrational spectrum of water dimer is challenging. The potential has eight wells separated by low barriers, which makes harmonic approximations of limited utility. A variational approach is imperative, but difficult because there are 12 coupled vibrational coordinates. In this paper, we use a product contracted basis whose functions are products of intramolecular and intermolecular functions computed using an iterative eigensolver. An intermediate matrix F facilitates calculating matrix elements. Using F, it is possible to do calculations on a general potential without storing the potential on the full quadrature grid. We find that surprisingly many intermolecular functions are required. This is due to the importance of coupling between inter- and intra-molecular coordinates. The full G16 symmetry of water dimer is exploited. We calculate, for the first time, monomer excited stretch states and compare P(1) transition frequencies with their experimental counterparts. We also compare with experimental vibrational shifts and tunneling splittings. Surprisingly, we find that the largest tunneling splitting, which does not involve the interchange of the two monomers, is smaller in the asymmetric stretch excited state than in the ground state. Differences between levels we compute and those obtained with a [6+6]D adiabatic approximation [Leforestier et al. J. Chem. Phys. 137 014305 (2012)] are ∼ 0.6 cm−1 for states without monomer excitation, ∼ 4 cm−1 for monomer excited bend states, and as large as ∼ 10 cm−1 for monomer excited stretch states. [ABSTRACT FROM AUTHOR]

Until quite recently, in almost all papers on crystal nucleation in glass-forming substances, it was assumed that nucleation proceeds in a completely relaxed supercooled liquid and, hence, at constant values of the critical parameters determining the nucleation rate for any given set of temperature, pressure, and composition. Here, we analyze the validity of this hypothesis for a model system by studying nucleation in a lithium silicate glass treated for very long times (up to 250 days) in deeply supercooled states, reaching 60 K below the laboratory glass transition temperature, Tg. At all temperatures in the considered range, T < Tg, we observed an enormous difference between the experimental number of nucleated crystals, N(t), and its theoretically expected value computed by assuming the metastable state of the relaxing glass has been reached. Analyzing the origin of this discrepancy, we confirmed that the key parameters determining the nucleation rates change with time as a result of the glass relaxation process. Finally, we demonstrate that, for temperatures below 683 K, this particular glass almost fully crystallizes prior to reaching the ultimate steady-state nucleation regime (e.g., at 663 K, it would take 176 years for the glass to reach 99% crystallization, while 2600 years would be needed for complete relaxation). This comprehensive study proves that structural relaxation strongly affects crystal nucleation in deeply supercooled states at temperatures well below Tg; hence, this phenomenon has to be accounted for in any crystal nucleation model. [ABSTRACT FROM AUTHOR]

The first paper of this series [J. Chem. Phys. 158, 034103 (2023)] demonstrated that excess entropy scaling holds for both fine-grained and corresponding coarse-grained (CG) systems. Despite its universality, a more exact determination of the scaling relationship was not possible due to the semi-empirical nature. In this second paper, an analytical excess entropy scaling relation is derived for bottom-up CG systems. At the single-site CG resolution, effective hard sphere systems are constructed that yield near-identical dynamical properties as the target CG systems by taking advantage of how hard sphere dynamics and excess entropy can be analytically expressed in terms of the liquid packing fraction. Inspired by classical equilibrium perturbation theories and recent advances in constructing hard sphere models for predicting activated dynamics of supercooled liquids, we propose a new approach for understanding the diffusion of molecular liquids in the normal regime using hard sphere reference fluids. The proposed "fluctuation matching" is designed to have the same amplitude of long wavelength density fluctuations (dimensionless compressibility) as the CG system. Utilizing the Enskog theory to derive an expression for hard sphere diffusion coefficients, a bridge between the CG dynamics and excess entropy is then established. The CG diffusion coefficient can be roughly estimated using various equations of the state, and an accurate prediction of accelerated CG dynamics at different temperatures is also possible in advance of running any CG simulation. By introducing another layer of coarsening, these findings provide a more rigorous method to assess excess entropy scaling and understand the accelerated CG dynamics of molecular fluids. [ABSTRACT FROM AUTHOR]

This paper reports the derivation and implementation of the electric dipole-magnetic dipole and electric dipole-electric quadrupole polarizability tensors at the density functional theory level with periodic boundary conditions (DFT-PBC). These tensors are combined to evaluate the Buckingham/Dunn tensor that describes the optical rotation (OR) in oriented chiral systems. We describe several aspects of the derivation of the equations and present test calculations that verify the correctness of the tensor formulation and their implementation. The results show that the full OR tensor is completely origin invariant as for molecules and that PBC calculations match molecular cluster calculations on 1D chains. A preliminary investigation on the choice of density functional, basis set, and gauge indicates a similar dependence as for molecules: the functional is the primary factor that determines the OR magnitude, followed by the basis set and to a much smaller extent the choice of gauge. However, diffuse functions may be problematic for PBC calculations even if they are necessary for the molecular case. A comparison with experimental data of OR for the tartaric acid crystal shows reasonable agreement given the level of theory employed. The development presented in this paper offers the opportunity to simulate the OR of chiral crystalline materials with general-purpose DFT-PBC methods, which, in turn, may help to understand the role of intermolecular interactions on this sensitive electronic property. [ABSTRACT FROM AUTHOR]

Two-dimensional transition metal dichalcogenides and semiconductor metal oxides have shown great potential in photocatalysis. However, their stability and efficiency need to be further improved. In this paper, porous ZnO nanorods with high specific surface area were prepared from metal-organic framework ZIF-8 by a simple hydrothermal method. A MoS2/ZnO composite was constructed by loading MoS2 onto the surface of porous ZnO nanorods. Compared with ZnO materials prepared by other methods, MoS2/ZnO prepared in this paper exhibits superior photocatalytic performance. The enhanced photocatalytic activity of the MoS2/ZnO composite can be attributed to the formation of heterojunctions and strong interaction between them, which greatly facilitate the separation of electrons and holes at the contact interface. In addition, due to the wide absorption region of the visible spectrum, MoS2 can greatly broaden the light absorption range of the material after the formation of the composite material, increase the utilization rate of visible light, and reduce the combination of electrons and holes. This study provides a new way to prepare cheap and efficient photocatalysts. [ABSTRACT FROM AUTHOR]

Fluorescence-encoded infrared (FEIR) spectroscopy is an emerging technique for performing vibrational spectroscopy in solution with detection sensitivity down to single molecules. FEIR experiments use ultrashort pulses to excite a fluorescent molecule's vibrational and electronic transitions in a sequential, time-resolved manner, and are therefore sensitive to intervening vibrational dynamics on the ground state, vibronic coupling, and the relative orientation of vibrational and electronic transition dipole moments. This series of papers presents a theoretical treatment of FEIR spectroscopy that describes these phenomena and examines their manifestation in experimental data. This first paper develops a nonlinear response function description of Fourier-transform FEIR experiments for a two-level electronic system coupled to multiple vibrations, which is then applied to interpret experimental measurements in the second paper [L. Whaley-Mayda et al., J. Chem. Phys. 159, 194202 (2023)]. Vibrational coherence between pairs of modes produce oscillatory features that interfere with the vibrations' population response in a manner dependent on the relative signs of their respective Franck–Condon wavefunction overlaps, leading to time-dependent distortions in FEIR spectra. The orientational response of population and coherence contributions are analyzed and the ability of polarization-dependent experiments to extract relative transition dipole angles is discussed. Overall, this work presents a framework for understanding the full spectroscopic information content of FEIR measurements to aid data interpretation and inform optimal experimental design. [ABSTRACT FROM AUTHOR]

Developing fluorescence-encoded infrared (FEIR) vibrational spectroscopy for single-molecule applications requires a detailed understanding of how the molecular response and external experimental parameters manifest in the detected signals. In Paper I [L. Whaley-Mayda, A. Guha, and A. Tokmakoff, J. Chem. Phys. 159, 194201 (2023)] we introduced a nonlinear response function theory to describe vibrational dynamics, vibronic coupling, and transition dipole orientation in FEIR experiments with ultrashort pulses. In this second paper, we apply the theory to investigate the role of intermode vibrational coherence, the orientation of vibrational and electronic transition dipoles, and the effects of finite pulse durations in experimental measurements. We focus on measurements at early encoding delays—where signal sizes are largest and therefore of most value for single-molecule experiments, but where many of these phenomena are most pronounced and can complicate the appearance of data. We compare experiments on coumarin dyes with finite-pulse response function simulations to explain the time-dependent behavior of FEIR spectra. The role of the orientational response is explored by analyzing polarization-dependent experiments and their ability to resolve relative dipole angles in the molecular frame. This work serves to demonstrate the molecular information content of FEIR experiments, and develop insight and guidelines for their interpretation. [ABSTRACT FROM AUTHOR]

Double Resonance is a powerful spectroscopic method that unambiguously assigns the rigorous quantum numbers of one state of a transition. However, there is often ambiguity as to the branch (ΔJ) of that transition. Spectroscopists have resolved this ambiguity by using the dependence of the double resonance intensity on the relative polarization directions of pump and probe radiation. However, published theoretical predictions for this ratio are based upon a weak (i.e., non-saturating) field approximation. This paper presents theoretical predictions for these intensity ratios for cases where the pump field is strongly saturating in the two limits of transitions dominated by homogeneous or of inhomogeneous broadening. Saturation reduces but does not eliminate the magnitude of the polarization effect (driving the intensity ratio closer to unity) even with strong pump saturation. For the case of an inhomogeneously broadened line, such as when Doppler broadened linewidth dominates over the power-broadened homogeneous line width, a large fraction of the low pump power polarization anisotropy remains. This paper reports predicted polarization ratios for both linear and circular pump and probe field polarizations. The present predictions are compared with experimental measurements on CH4 ground state → ν3 → 3ν3 transitions recently reported by de Oliveira et al.63 and these are in better agreement than with the weak field predictions. [ABSTRACT FROM AUTHOR]

This paper presents an exploration of a specific type of a generalized multicomponent solution model, which appears to be first given by Saulov in the current explicit form. The assumptions of the underlying theory and a brief derivation of the main equation have been provided preliminarily for completeness and notational consistency. The resulting formulae for the Gibbs free energy of mixing and the chemical potentials are multivariate polynomials with physically meaningful coefficients and the mole fractions of the components as variables. With one additional assumption about the relative magnitudes of the solvent–solute and solute–solute interaction exchange energies, combining rules were obtained that express the mixed coefficients of the polynomial in terms of its pure coefficients. This was done by exploiting the mathematical structure of the asymmetric form of the solvent chemical potential equation. The combining rules allow one to calculate the thermodynamic properties of the solvent with multiple solutes from binary mixture data only (i.e., each solute with the solvent), and hence, are of practical importance. Furthermore, a connection was established between the osmotic virial coefficients derived in this work and the original osmotic virial coefficients of Hill found by employing a different procedure, illustrating the equivalency of what appears to be two different theories. A validation of the combining rules derived here has been provided in a separate paper where they were successfully used to predict the freezing points of ternary salt solutions of water. [ABSTRACT FROM AUTHOR]

This paper series aims to establish a complete correspondence between fine-grained (FG) and coarse-grained (CG) dynamics by way of excess entropy scaling (introduced in Paper I). While Paper II successfully captured translational motions in CG systems using a hard sphere mapping, the absence of rotational motions in single-site CG models introduces differences between FG and CG dynamics. In this third paper, our objective is to faithfully recover atomistic diffusion coefficients from CG dynamics by incorporating rotational dynamics. By extracting FG rotational diffusion, we unravel, for the first time reported to our knowledge, a universality in excess entropy scaling between the rotational and translational diffusion. Once the missing rotational dynamics are integrated into the CG translational dynamics, an effective translation-rotation coupling becomes essential. We propose two different approaches for estimating this coupling parameter: the rough hard sphere theory with acentric factor (temperature-independent) or the rough Lennard-Jones model with CG attractions (temperature-dependent). Altogether, we demonstrate that FG diffusion coefficients can be recovered from CG diffusion coefficients by (1) incorporating "entropy-free" rotational diffusion with translation-rotation coupling and (2) recapturing the missing entropy. Our findings shed light on the fundamental relationship between FG and CG dynamics in molecular fluids. [ABSTRACT FROM AUTHOR]

Protein folding in confined media has attracted wide attention over the past 15 years due to its importance to both in vivo and in vitro applications. It is generally believed that protein stability increases by decreasing the size of the confining medium, if the medium’s walls are repulsive, and that the maximum folding temperature in confinement is in a pore whose size D0 is only slightly larger than the smallest dimension of a protein’s folded state. Until recently, the stability of proteins in pores with a size very close to that of the folded state has not received the attention it deserves. In a previous paper [L. Javidpour and M. Sahimi, J. Chem. Phys. 135, 125101 (2011)], we showed that, contrary to the current theoretical predictions, the maximum folding temperature occurs in larger pores for smaller α-helices. Moreover, in very tight pores, the free energy surface becomes rough, giving rise to a new barrier for protein folding close to the unfolded state. In contrast to unbounded domains, in small nanopores proteins with an α-helical native state that contain the β structures are entropically stabilized implying that folding rates decrease notably and that the free energy surface becomes rougher. In view of the potential significance of such results to interpretation of many sets of experimental data that could not be explained by the current theories, particularly the reported anomalously low rates of folding and the importance of entropic effects on proteins’ misfolded states in highly confined environments, we address the following question in the present paper: To what extent the geometry of a confined medium affects the stability and folding rates of proteins? Using millisecond-long molecular dynamics simulations, we study the problem in three types of confining media, namely, cylindrical and slit pores and spherical cavities. Most importantly, we find that the prediction of the previous theories that the dependence of the maximum folding temperature Tf on the size D of a confined medium occurs in larger media for larger proteins is correct only in spherical geometry, whereas the opposite is true in the two other geometries that we study. Also studied is the effect of the strength of the interaction between the confined media’s walls and the proteins. If the walls are only weakly or moderately attractive, a complex behavior emerges that depends on the size of the confining medium. [ABSTRACT FROM AUTHOR]

In Paper I, the performances of pre-screening (PS), extended PS (EPS), and cumulant (CU) approximations to the fourth-order density matrix were examined in the context of second-order N-electron valence state perturbation theory (NEVPT2). It has been found that the CU, PS, and even EPS approximations with loose thresholds may introduce intruder states. In the present work, the origin of these "false intruder" states introduced by approximated density matrices is discussed. Canonical NEVPT2 implementations employ a rank reduction trick. By analyzing its residual error, we find that the omission of the rank reduction leads to a more stable multireference perturbation theory for incomplete active space reference wave functions. Such a full rank (FR)-NEVPT2 formulation is equivalent to the conventional NEVPT2 method for the complete active space self-consistent field/complete active space configuration interaction reference wave function. A major drawback of the FR-NEVPT2 formulation is the necessity of the fifth-order density matrix. To avoid the construction of the high-order density matrices, the combination of the FR-NEVPT2 with the CU approximation is studied. However, we find that the CU approximation remains problematic as it still introduces intruder states. The question of how to robustly and efficiently perform internally contracted multireference perturbation theories with approximate densities remains a challenging field of investigation. [ABSTRACT FROM AUTHOR]

This paper proposes the first linear scaling implementation for the domain decomposition approach of the polarizable continuum model (ddPCM) for the computation of the solvation energy and forces. The ddPCM-equation consists of a (non-local) integral equation on the van der Waals or solvent accessible surface of the solute's cavity resulting in a dense solution matrix, and, in turn, one matrix–vector multiplication has a quadratic arithmetic complexity with respect to the number of atoms of the solute molecule. The use of spherical harmonics as basis functions makes it natural to employ the fast multipole method (FMM) in order to provide an asymptotically linear scaling method. In this paper, we employ the FMM in a non-uniform manner with a clusterization based on a recursive inertial bisection. We present some numerical tests illustrating the accuracy and scaling of our implementation. [ABSTRACT FROM AUTHOR]

Room temperature ionic liquids (RTILs) are molten salts consisting entirely of ions and have over the past decades gained increased interest due to their high potential in applications. These structurally complex systems often display multiple relaxation modes in the response functions at lower frequencies, hinting to complex underlying mechanisms. While the existence of these multimodal spectra in the shear mechanical, dielectric, and light scattering response of RTILs has been confirmed multiple times, controversy still surrounds the origin. This paper, therefore, aims to provide additional insights into the multimodal spectra seen in RTILs by presenting new shear mechanical results on seven different RTILs: Pyr1n-TFSI with n = 4, 6, and 8; Pyr18-TFSI mixed with Li-TFSI in two high concentrations; and Cn-mim-BF4 with n = 3 and 8. Dynamic depolarized light scattering was also measured on one of the Pyr18-TFSI Li-salt mixtures. These specific cases were analyzed in detail and put into a bigger perspective together with an overview of the literature. Recent literature offers two specific explanations for the origin of the multimodal shear mechanical spectra: (1) cation–anion time scale separation or (2) combined cation–anion relaxation in addition to a dynamic signal from mesoscale aggregates at lower frequencies. However, neither of these two pictures can consistently explain all the results on different ionic liquids. Instead, we conclude that the origin of the multimodal spectrum is system specific. This underlines the complexity of this class of liquids and shows that great care must be taken when making general conclusions based on specific cases. [ABSTRACT FROM AUTHOR]

The doping of CdS quantum dots (QDs) with Cu(I) disrupts electron–hole correlation due to hole trapping by the dopant ion, post-photoexcitation. The present paper examines the effect of such disruption on the rate of photoinduced electron transfer (PET) from the QDs to methyl viologen (MV2+), with implications in their photocatalytic activity. A significantly greater efficiency of PL quenching by MV2+ is observed for the doped QDs than for the undoped ones. Interestingly, the Stern–Volmer plots constructed using PL intensities exhibit an upward curvature for both the cases, while the PL lifetimes remain unaffected. This observation is rationalized by considering the adsorption of the quencher on the surface of the QDs and ultrafast PET post-photoexcitation. Ultrafast transient absorption experiments confirm a faster electron transfer for the doped QDs. It is also realized that the transient absorption experiment yields a more accurate estimate of the binding constant of the quencher with the QDs, than the PL experiment. [ABSTRACT FROM AUTHOR]

Fermi's golden rule (FGR) offers an empirical framework for understanding the dynamics of spin–lattice relaxation in magnetic molecules, encompassing mechanisms like direct (one-phonon) and Raman (two-phonon) processes. These principles effectively model experimental longitudinal relaxation rates, denoted as T 1 − 1 . However, under scenarios of increased coupling strength and nonlinear spin–lattice interactions, FGR's applicability may diminish. This paper numerically evaluates the exact spin–lattice relaxation rate kernels, employing the extended dissipaton equation of motion formalism. Our calculations reveal that when quadratic spin–lattice coupling is considered, the rate kernels exhibit a free induction decay-like feature, and the damping rates depend on the interaction strength. We observe that the temperature dependence predicted by FGR significantly deviates from the exact results since FGR ignores the higher order effects and the non-Markovian nature of spin–lattice relaxation. Our methods can be easily extended to study other systems with nonlinear spin–lattice interactions and provide valuable insights into the temperature dependence of T1 in molecular qubits when the coupling is strong. [ABSTRACT FROM AUTHOR]

The biochemical activity inside a cell has recently been suggested to act as a source of hydrodynamic fluctuations that can speed up or slow down enzyme catalysis [Tripathi et al., Commun. Phys. 5, 101 (2022).] The idea has been tested against and largely corroborated by simulations of activated barrier crossing in a simple fluid in the presence of thermal and athermal noise. The present paper attempts a wholly analytic solution to the same noise-driven barrier crossing problem but generalizes it to include viscoelastic memory effects of the kind likely to be present in cellular interiors. A calculation of the model's barrier crossing rate, using Kramers' flux-over-population formalism, reveals that in relation to the case where athermal noise is absent, athermal noise always accelerates barrier crossing, though the extent of enhancement depends on the duration τ0 over which the noise acts. More importantly, there exists a critical τ0—determined by the properties of the medium—at which Kramers' theory breaks down and, on approach to which, the rate grows significantly. The possibility of such a giant enhancement is potentially open to experimental validation using optically trapped nanoparticles in viscoelastic media that are acted on by externally imposed colored noise. [ABSTRACT FROM AUTHOR]

This paper presents the first implementation of a coupling between advanced wavefunction theories and molecular density functional theory (MDFT). This method enables the modeling of solvent effect into quantum mechanical (QM) calculations by incorporating an electrostatic potential generated by solvent charges into the electronic Hamiltonian. Solvent charges are deduced from the spatially and angularly dependent solvent particle density. Such a density is obtained through the minimization of the functional associated with the molecular mechanics (MM) Hamiltonian describing the interaction between the fluid particles. The introduced QM/MDFT framework belongs to QM/MM family of methods, but its originality lies in the use of MDFT as the MM solver, offering two main advantages. First, its functional formulation makes it competitive with respect to sampling-based molecular mechanics. Second, it preserves a molecular-level description lost in macroscopic continuum approaches. The excited state properties of water and formaldehyde molecules solvated into water have been computed at the selected configuration interaction (SCI) level. The excitation energies and dipole moments have been compared with experimental data and previous theoretical work. A key finding is that using the Hartree–Fock method to describe the solute allows for predicting the solvent charge around the ground state with sufficient precision for the subsequent SCI calculations of excited states. This significantly reduces the computational cost of the described procedure, paving the way for the study of more complex molecules. [ABSTRACT FROM AUTHOR]

Electron transfer reactions play an essential role in many chemical and biological processes. Fermi's golden rule (GR), which assumes that the coupling between electronic states is small, has formed the foundation of electron transfer rate theory; however, in short range electron/energy transfer reactions, this coupling can become very large, and, therefore, Fermi's GR fails to make even qualitatively accurate rate predictions. In this paper, I present a simple modified GR theory to describe electron transfer in the Marcus inverted regime at arbitrarily large electronic coupling strengths. This theory is based on an optimal global rotation of the diabatic states, which makes it compatible with existing methods for calculating GR rates that can account for nuclear quantum effects with anharmonic potentials. Furthermore, the optimal GR (OGR) theory can also be combined with analytic theories for non-adiabatic rates, such as Marcus theory and Marcus–Levich–Jortner theory, offering clear physical insights into strong electronic coupling effects in non-adiabatic processes. OGR theory is also tested on a large set of spin-boson models and an anharmonic model against exact quantum dynamics calculations, where it performs well, correctly predicting rate turnover at large coupling strengths. Finally, an example application to a boron-dipyrromethane–anthracene photosensitizer reveals that strong coupling effects inhibit excited state charge recombination in this system, reducing the rate of this process by a factor of 4. Overall, OGR theory offers a new approach to calculating electron transfer rates at strong couplings, offering new physical insights into a range of non-adiabatic processes. [ABSTRACT FROM AUTHOR]

Two-dimensional (2D) Janus transition metal dichalcogenides MXY manifest novel physical properties owing to the breaking of out-of-plane mirror symmetry. Recently, the 2H phase of MoSH has been demonstrated to possess intrinsic superconductivity, whereas the 1T phase exhibits a charge density waves state. In this paper, we have systematically studied the stability and electron–phonon interaction characteristics of MoSLi. Our results have shown that both the 2H and 1T phases of MoSLi are stable, as indicated by the phonon spectrum and the ab initio molecular dynamics. However, the 1T phase exhibits an electron–phonon coupling constant that is twice as large as that of the 2H phase. In contrast to MoSH, the 1T phase of MoSLi exhibits intrinsic superconductivity. By employing the ab initio anisotropic Migdal-Eliashberg formalism, we have revealed the two-gap superconducting nature of 1T-MoSLi, with a transition temperature (Tc) of 14.8 K. The detailed analysis indicates that the superconductivity in 1T-MoSLi primarily originates from the interplay between the vibration of the phonon modes in the low-frequency region and the dz2 orbital. These findings provide a fresh perspective on superconductivity within Janus structures. [ABSTRACT FROM AUTHOR]

Vibrational coherences in ultrafast pump–probe (PP) and 2D electronic spectroscopy (2DES) provide insights into the excited state dynamics of molecules. Femtosecond coherence spectra and 2D beat maps yield information about displacements of excited state surfaces for key vibrational modes. Half-broadband 2DES uses a PP configuration with a white light continuum probe to extend the detection range and resolve vibrational coherences in the excited state absorption (ESA). However, the interpretation of these spectra is difficult as they are strongly dependent on the spectrum of the pump laser and the relative displacement of the excited states along the vibrational coordinates. We demonstrate the impact of these convoluting factors for a model based upon cresyl violet. A careful consideration of the position of the pump spectrum can be a powerful tool in resolving the ESA coherences to gain insights into excited state displacements. This paper also highlights the need for caution in considering the spectral window of the pulse when interpreting these spectra. [ABSTRACT FROM AUTHOR]

It is significant to investigate the calcium carbonate (CaCO3) precipitation mechanism during the carbon capture process; nevertheless, CaCO3 precipitation is not clearly understood yet. Understanding the carbonation mechanism at the atomic level can contribute to the mineralization capture and utilization of carbon dioxide, as well as the development of new cementitious materials with high-performance. There are many factors, such as temperature and CO2 concentration, that can influence the carbonation reaction. In order to achieve better carbonation efficiency, the reaction conditions of carbonation should be fully verified. Therefore, based on molecular dynamics simulations, this paper investigates the atomic-scale mechanism of carbonation. We investigate the effect of carbonation factors, including temperature and concentration, on the kinetics of carbonation (polymerization rate and activation energy), the early nucleation of calcium carbonate, etc. Then, we analyze the local stresses of atoms to reveal the driving force of early stage carbonate nucleation and the reasons for the evolution of polymerization rate and activation energy. Results show that the higher the calcium concentration or temperature, the higher the polymerization rate of calcium carbonate. In addition, the activation energies of the carbonation reaction increase with the decrease in calcium concentrations. [ABSTRACT FROM AUTHOR]

Endohedrally doped clusters form a large category of cage clusters, with unique structures, diverse elemental compositions, and highly tunable electronic structures and physisochemical properties. They have been widely achieved in laboratory and may serve as functional building blocks for assembling new supermolecular structures and devices. In this paper, for the first time, we disclosed the luminescence properties of endohedrally doped group-IV clusters by time-dependent density functional theory calculations. A total of 64 cage clusters have been explored in terms of stability, emission wavelength, and the energy difference between the first excited singlet and triplet states. The key geometric and electronic factors governing the photophysical properties of these cage clusters were unveiled, to provide crucial insights for crafting atomically precise nanoclusters for optical and optoelectronic applications. [ABSTRACT FROM AUTHOR]

Quantum Dot Cellular Automata (QCA) is an emerging trend in the field of nanoelectronics, and computing can be regarded as an alternative to the traditional complementary metal–oxide–semiconductor technology. The paper is devoted to the study of the key functional properties of the cells for molecular QCA based on mixed valence molecules. The theoretical results for the heat dissipation under the conditions of the fast nonadiabatic switching event and cell–cell response function are obtained in the framework of the quantum-mechanical vibronic approach. These results are parameterized using the previous reliable ab initio calculations performed for oxidized norbornadiene and its polycyclic derivatives with variable lengths of the bridge. The comparative analysis of the dimeric and bidimeric molecular cells composed of these compounds is given. It is underlined that the conditions of a strong non-linear response and a low heat release are contradictory. However, despite this problem, a parametric regime is proposed, which provides a low heat release in combination with a strong nonlinear response of the working cell to the electric field induced by the polarized driver cell. [ABSTRACT FROM AUTHOR]

This paper evaluates the accuracy of the Hermitian form of the downfolding procedure using the double unitary coupled cluster (DUCC) ansatz on the benchmark systems of linear chains of hydrogen atoms, H6 and H8. The computational infrastructure employs the occupation-number-representation codes to construct the matrix representation of arbitrary second-quantized operators, allowing for the exact representation of exponentials of various operators. The tests demonstrate that external amplitudes from standard single-reference coupled cluster methods that sufficiently describe external (out-of-active-space) correlations reliably parameterize the Hermitian downfolded effective Hamiltonians in the DUCC formalism. The results show that this approach can overcome the problems associated with losing the variational character of corresponding energies in the corresponding SR-CC theories. [ABSTRACT FROM AUTHOR]

The superadiabatic dynamical density functional theory (superadiabatic-DDFT) is a promising new method for the study of colloidal systems out-of-equilibrium. Within this approach, the viscous forces arising from interparticle interactions are accounted for in a natural way by explicitly treating the dynamics of the two-body correlations. For bulk systems subject to spatially homogeneous shear, we use the superadiabatic-DDFT framework to calculate the steady-state pair distribution function and the corresponding viscosity for low values of the shear-rate. We then consider a variant of the central approximation underlying this superadiabatic theory and obtain an inhomogeneous generalization of a rheological bulk theory due to Russel and Gast. This paper thus establishes for the first time a connection between DDFT approaches, formulated to treat inhomogeneous systems, and existing work addressing nonequilibrium microstructure and rheology in bulk colloidal suspensions. [ABSTRACT FROM AUTHOR]

We propose a new collocation multi-configuration time-dependent Hartree (MCTDH) method. It reduces point-set error by using more points than basis functions. Collocation makes it possible to use MCTDH with a general potential energy surface without computing any integrals. The collocation points are associated with a basis larger than the basis used to represent wavefunctions. Both bases are obtained from a direct product basis built from single-particle functions by imposing a pruning condition. The collocation points are those on a sparse grid. Heretofore, collocation MCTDH calculations with more points than basis functions have only been possible if both the collocation grid and the basis set are direct products. In this paper, we exploit a new pseudo-inverse to use both more points than basis functions and a pruned basis and grid. We demonstrate that, for a calculation of the lowest 50 vibrational states (energy levels and wavefunctions) of CH2NH, errors can be reduced by two orders of magnitude by increasing the number of points, without increasing the basis size. This is true also when unrefined time-independent points are used. [ABSTRACT FROM AUTHOR]

In the previous paper [C. W. Kim and I. Franco, J. Chem. Phys. 160, 214111-1–214111-13 (2024)], we developed a theory called MQME-D, which allows us to decompose the overall energy dissipation process in open quantum system dynamics into contributions by individual components of the bath when the subsystem dynamics is governed by a Markovian quantum master equation (MQME). Here, we contrast the predictions of MQME-D against the numerically exact results obtained by combining hierarchical equations of motion (HEOM) with a recently reported protocol for monitoring the statistics of the bath. Overall, MQME-D accurately captures the contributions of specific bath components to the overall dissipation while greatly reducing the computational cost compared to exact computations using HEOM. The computations show that MQME-D exhibits errors originating from its inherent Markov approximation. We demonstrate that its accuracy can be significantly increased by incorporating non-Markovianity by exploiting time scale separations (TSS) in different components of the bath. Our work demonstrates that MQME-D combined with TSS can be reliably used to understand how energy is dissipated in realistic open quantum system dynamics. [ABSTRACT FROM AUTHOR]