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]

Several new variants of the locally-renormalized coupled-cluster (CC) approaches that account for the effect of triples (LR-CCSD(T)) have been formulated and implemented for arbitrary reference states using the TENSOR CONTRACTION ENGINE functionality, enabling the automatic generation of an efficient parallel code. Deeply rooted in the recently derived numerator–denominator-connected (NDC) expansion for the ground-state energy [K. Kowalski and P. Piecuch, J. Chem. Phys. 122, 074107 (2005)], LR-CCSD(T) approximations use, in analogy to the completely renormalized CCSD(T) (CR-CCSD(T)) approach, the three-body moments in constructing the noniterative corrections to the energies obtained in CC calculations with singles and doubles (CCSD). In contrast to the CR-CCSD(T) method, the LR-CCSD(T) approaches discussed in this paper employ local denominators, which assure the additive separability of the energies in the noninteracting system limit when the localized occupied spin-orbitals are employed in the CCSD and LR-CCSD(T) calculations. As clearly demonstrated on several challenging examples, including breaking the bonds of the F2, N2, and CN molecules, the LR-CCSD(T) approaches are capable of providing a highly accurate description of the entire potential-energy surface (PES), while maintaining the characteristic N7scaling of the ubiquitous CCSD(T) approach. Moreover, as illustrated numerically for the ozone molecule, the LR-CCSD(T) approaches yield highly competitive values for a number of equilibrium properties including bond lengths, angles, and harmonic frequencies. [ABSTRACT FROM AUTHOR]

The paper referred to in the title, in J. Chem. Phys. 84, 1788 (1986), presented results of computer simulations of non-Markovian barrier crossing dynamics, and compared these results with various theoretical predictions. One comment made here is that the observed limited applicability of the Grote–Hynes theory can be explained by singular perturbation theory. The other comment concerns the importance of a second constant of the motion when the non-Markovian friction has a very long relaxation time. [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]

Density-functional-theory-based chemical reactivity indicators are formulated for degenerate and near-degenerate ground states. For degenerate states, the functional derivatives of the energy with respect to the external potential do not exist, and must be replaced by the weaker concept of functional variation. The resultant reactivity indicators depend on the specific perturbation. Because it is sometimes impractical to compute reactivity indicators for a specific perturbation, we consider two special cases: point-charge perturbations and Dirac delta function perturbations. The Dirac delta function perturbations provide upper bounds on the chemical reactivity. Reactivity indicators using the common used 'average of degenerate states approximation' for degenerate states provide a lower bound on the chemical reactivity. Unfortunately, this lower bound is often extremely weak. Approximate formulas for the reactivity indicators within the frontier-molecular-orbital approximation and special cases (two or three degenerate spatial orbitals) are presented in the supplementary material. One remarkable feature that arises in the frontier molecular orbital approximation, and presumably also in the exact theory, is that removing electrons sometimes causes the electron density to increase at the location of a negative (attractive) Dirac delta function perturbation. That is, the energetic response to a reduction in the external potential can increase even when the number of electrons decreases. [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]

A theoretical treatment of the vibrational–vibrational (VV) contribution to the vibrational energy relaxation time T1 of a solute normal mode in a molecular solvent, which extends a previous treatment [S. A. Adelman, R. H. Stote, and R. Muralidhar, J. Chem. Phys. 99, 1320 (1993), henceforth called Paper I] of the vibrational–translational–rotational (VTR) contribution to T1, is outlined and expressions for this VV contribution, valid for the infinitely dilute diatomic solution, are presented. The treatment is based on the formula T1=β-1(ωl), where β(ω) is the friction kernel of the relaxing solute mode and where ωl is its liquid phase frequency. β(ω) is evaluated as the cosine transform of the autocorrelation function 〈F(t)Fĩ〉0v of the fluctuating generalized force exerted by the vibrating solvent on the solute normal mode coordinate conditional that this coordinate is fixed at its equilibrium value. 〈F(t)F〉0v is expressed as a superposition of the rigid solvent autocorrelation function 〈F(t)F〉0 and a correction which accounts for solvent vibrational motion. For diatomic solvents one has 〈F(t)F〉0v= 〈F(t)F〉0+NSMD(t) cos ωDt F(ΩD), where NS=number of solvent molecules, MD(t) is the vibrational force gradient autocorrelation function, ωD and ΩD are solvent molecule liquid phase frequencies, and F(Ω)=1/2hΩ-1 coth[hΩ/2kBT].The Gaussian model is assumed for 〈F(t)F〉0 and MD(t) yielding β(ω) as a superposition of a Gaussian centered at ω=0 which mediates VTR processes and a Gaussian centered at ω=ωD which mediates VV processes. Vector integral expressions for MD(t), ωD, and ΩD are presented which are similar to the expressions for ωl and 〈F(t)F〉0 given in Paper I. These expressions permit the evaluation of the VV contribution to T1 from the atomic masses, bondlengths, vibrational frequencies, and site–site interaction... [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]

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]

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]