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]

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]

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]

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]

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]