Search

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]

A complete set of four polarization-dependent differential cross sections in the reactions of Cl + aligned-CHD3(v1 = 1, |jK) → DCl(v = 0) + CHD2(v1 = 1) is reported here for two different, rotationally polarized states with j = 1: specifically the |jK = |10 state prepared via the R(0) excitation and the |1 ± 1 state via Q(1). In stark contrast to the complicated situation of the HCl(v) + CD3(v = 0) channel reported in Paper-I, the stereo-requirement of this isotopic channel for both polarized reactants appears quite straightforward and consistent with a direct rebound mechanism. The extent of steric effects is moderate and relatively smaller than the alternative H-atom abstraction channel. All major findings reported here can qualitatively be understood by first noting that the present reaction invokes abstracting a D-atom, which is the spectator in the IR-excitation process. Next, it is recognized that the directional properties of two polarized states of CHD3(v1 = 1, |jK) should manifest primarily in the IR-excited C-H bond, leaving secondary imprints in the unexcited CD3-moiety. The stereo-specificity of the DCl + CHD2 product channel is further reduced by the fact that the abstraction can occur with any one of the three spatially distinct D-atoms. [ABSTRACT FROM AUTHOR]

We obtain the phase diagram of the hard core lattice gas with third nearest neighbor exclusion on the triangular lattice using Monte Carlo simulations that are based on a rejection-free flat histogram algorithm. In a recent paper [Darjani et al., J. Chem. Phys. 151, 104702 (2019)], it was claimed that the lattice gas with third nearest neighbor exclusion undergoes two phase transitions with increasing density with the phase at intermediate densities exhibiting hexatic order with continuously varying exponents. Although a hexatic phase is expected when the exclusion range is large, it has not been seen earlier in hard core lattice gases with short range exclusion. In this paper, by numerically determining the entropies for all densities, we show that there is only a single phase transition in the system between a low-density fluid phase and a high density ordered sublattice phase and that a hexatic phase is absent. The transition is shown to be first order in nature, and the critical parameters are determined accurately. [ABSTRACT FROM AUTHOR]

Lyotropic liquid crystal phases (LCPs) are widely studied for diverse applications, including protein crystallization and drug delivery. The structure and properties of LCPs vary widely depending on the composition, concentration, temperature, pH, and pressure. High-throughput structural characterization approaches, such as small-angle x-ray scattering (SAXS), are important to cover meaningfully large compositional spaces. However, high-throughput LCP phase analysis for SAXS data is currently lacking, particularly for patterns of multiphase mixtures. In this paper, we develop semi-automated software for high throughput LCP phase identification from SAXS data. We validate the accuracy and time-savings of this software on a total of 668 SAXS patterns for the LCPs of the amphiphile hexadecyltrimethylammonium bromide (CTAB) in 53 acidic or basic ionic liquid derived solvents, within a temperature range of 25–75 °C. The solvents were derived from stoichiometric ethylammonium nitrate (EAN) or ethanolammonium nitrate (EtAN) by adding water to vary the ionicity, and adding precursor ions of ethylamine, ethanolamine, and nitric acid to vary the pH. The thermal stability ranges and lattice parameters for CTAB-based LCPs obtained from the semi-automated analysis showed equivalent accuracy to manual analysis, the results of which were previously published. A time comparison of 40 CTAB systems demonstrated that the automated phase identification procedure was more than 20 times faster than manual analysis. Moreover, the high throughput identification procedure was also applied to 300 unpublished scattering patterns of sodium dodecyl-sulfate in the same EAN and EtAN based solvents in this study, to construct phase diagrams that exhibit phase transitions from micellar, to hexagonal, cubic, and lamellar LCPs. The accuracy and significantly low analysis time of the high throughput identification procedure validates a new, rapid, unrestricted analytical method for the determination of LCPs. [ABSTRACT FROM AUTHOR]

The ability of biomolecules to exert forces on their surroundings or resist compression from the environment is essential in a variety of biologically relevant contexts. For filaments in the low-temperature limit and under a constant compressive force, Euler buckling theory predicts a sudden transition from a compressed state to a bent state in these slender rods. In this paper, we use a mean-field theory to show that if a semiflexible chain is compressed at a finite temperature with a fixed end-to-end distance (permitting fluctuations in the compressive forces), it exhibits a continuous phase transition to a buckled state at a critical level of compression. We determine a quantitatively accurate prediction of the transverse position distribution function of the midpoint of the chain that indicates this transition. We find that the mean compressive forces are non-monotonic as the extension of the filament varies, consistent with the observation that strongly buckled filaments are less able to bear an external load. We also find that for the fixed extension (isometric) ensemble, the buckling transition does not coincide with the local minimum of the mean force (in contrast to Euler buckling). We also show that the theory is highly sensitive to fluctuations in length in two dimensions and the buckling transition can still be accurately recovered by accounting for those fluctuations. These predictions may be useful in understanding the behavior of filamentous biomolecules compressed by fluctuating forces, relevant in a variety of biological contexts. [ABSTRACT FROM AUTHOR]

Vapor-phase ammonia, NH3(g), and hydrochloric acid, HCl(g), undergo a series of complex reactions, including nucleation and growth, to form solid ammonium chloride, NH4Cl(s). The counterdiffusional experiment, whereby HCl(g) and NH3(g) diffuse from opposite ends of a tube and react to form spatiotemporally complex patterns, has a rich history of study. In this paper, we combine experimental data, molecular simulations, and analysis and simulations of a partial differential equation model to address the questions of where the first unobserved vapor product NH4Cl(g) and visually observable precipitate NH4Cl(s) form and how these positions depend on experimental parameters. These analyses yield a consistent picture which involves a moving reaction front as well as previously unobserved heterogeneous nucleation, wall nucleation, and homogeneous nucleation. The experiments combined with modeling allow for an estimate of the heterogeneous and homogeneous nucleation thresholds for the vapor-to-solid phase transition. The results, synthesized with the literature on this vapor-to-particle reaction, inform a discussion of the details of the reaction mechanism, including the role of water, which concludes the paper. [ABSTRACT FROM AUTHOR]

General transformation expressions of the second-order non-adiabatic Hamiltonian of the atomic nuclei, including the kinetic-energy correction terms, are derived upon the change from laboratory-fixed Cartesian coordinates to general curvilinear coordinate systems commonly used in rovibrational computations. The kinetic-energy or so-called "mass-correction" tensor elements are computed with the stochastic variational method and floating explicitly correlated Gaussian functions for the H 2 + molecular ion in its ground electronic state. {Further numerical applications for the 4 He 2 + molecular ion are presented in the forthcoming paper, Paper II [E. Mátyus, J. Chem. Phys. 149, 194112 (2018)]}. The general, curvilinear non-adiabatic kinetic energy operator expressions are used in the examples, and non-adiabatic rovibrational energies and corrections are determined by solving the rovibrational Schrödinger equation including the diagonal Born–Oppenheimer as well as the mass-tensor corrections. [ABSTRACT FROM AUTHOR]

The purpose of this paper is to investigate several analytical methods of solving first passage (FP) problem for the Rouse model, a simplest model of a polymer chain. We show that this problem has to be treated as a multi-dimensional Kramers' problem, which presents rich and unexpected behavior. We first perform direct and forward-flux sampling (FFS) simulations and measure the mean first-passage time t(z) for the free end to reach a certain distance z away from the origin. The results show that the mean FP time is getting faster if the Rouse chain is represented by more beads. Two scaling regimes of t(z) are observed, with transition between them varying as a function of chain length. We use these simulation results to test two theoretical approaches. One is a well known asymptotic theory valid in the limit of zero temperature. We show that this limit corresponds to fully extended chain when each chain segment is stretched, which is not particularly realistic. A new theory based on the well known Freidlin-Wentzell theory is proposed, where dynamics is projected onto the minimal action path. The new theory predicts both scaling regimes correctly, but fails to get the correct numerical prefactor in the first regime. Combining our theory with the FFS simulations leads us to a simple analytical expression valid for all extensions and chain lengths. One of the applications of polymer FP problem occurs in the context of branched polymer rheology. In this paper, we consider the arm-retraction mechanism in the tube model, which maps exactly on the model we have solved. The results are compared to the Milner-McLeish theory without constraint release, which is found to overestimate FP time by a factor of 10 or more. [ABSTRACT FROM AUTHOR]

The rocksalt structure of ZnO has a very promising bandgap for optoelectronic applications. Unfortunately, this high-pressure phase is unstable under ambient conditions. This paper presents experimental results for rocksalt-type ZnO/MgO superlattices and theoretical considerations of the critical thickness of MgxZn1−xO layers. The correlations between the layer/spacer thickness ratio, elastic strain, chemical composition, and critical thickness are analyzed. The Matthews and Blakeslee model is revisited to find analytic conditions for the critical layer thickness resulting in phase transition. Our analysis shows that due to the decrease in misfit stresses below some critical limit, the growth of multiple quantum wells composed of rocksalt ZnO layers and MgO spacers is possible only for very large layer/spacer thickness ratios. [ABSTRACT FROM AUTHOR]

It has been recently shown that the Tao–Mo (TM) [J. Tao and Y. Mo, Phys. Rev. Lett. 117, 073001 (2016)] semilocal exchange–correlation energy functional suffers from the order-of-limit problem, which affects the functional performance in phase transition pressures [Furness et al., J. Chem. Phys. 152, 244112 (2020)]. The root of the order-of-limit problem of the TM functional is inherent within the interpolation function, which acts as a switch between the compact density and the slowly varying density. This paper proposes a different switch function that avoids the order-of-limit problem and correctly interpolates the compact density and the slowly varying fourth-order density correction. By circumventing the order-of-limit problem, the proposed form enhances the applicability of the original TM functional on the diverse nature of solid-state properties. Our conclusion is ensured by examining the functional in predicting properties related to general-purpose solids, quantum chemistry, and phase transition pressure. Besides, we discuss the connection between the order-of-limit problem, phase transition pressure, and bandgap of solids. [ABSTRACT FROM AUTHOR]

Chemical reaction networks offer a natural nonlinear generalization of linear Markov jump processes on a finite state-space. In this paper, we analyze the dynamical large deviations of such models, starting from their microscopic version, the chemical master equation. By taking a large-volume limit, we show that those systems can be described by a path integral formalism over a Lagrangian functional of concentrations and chemical fluxes. This Lagrangian is dual to a Hamiltonian, whose trajectories correspond to the most likely evolution of the system given its boundary conditions. The same can be done for a system biased on time-averaged concentrations and currents, yielding a biased Hamiltonian whose trajectories are optimal paths conditioned on those observables. The appropriate boundary conditions turn out to be mixed so that, in the long time limit, those trajectories converge to well-defined attractors. We are then able to identify the largest value that the Hamiltonian takes over those attractors with the scaled cumulant generating function of our observables, providing a nonlinear equivalent to the well-known Donsker-Varadhan formula for jump processes. On that basis, we prove that chemical reaction networks that are deterministically multistable generically undergo first-order dynamical phase transitions in the vicinity of zero bias. We illustrate that fact through a simple bistable model called the Schlögl model, as well as multistable and unstable generalizations of it, and we make a few surprising observations regarding the stability of deterministic fixed points and the breaking of ergodicity in the large-volume limit. [ABSTRACT FROM AUTHOR]

Solid-state nudged elastic band (SSNEB) methods can be used for finding solid-solid transition paths when solids are subjected to external stress fields. However, previous SSNEB methods may lead to inaccurate barriers and deviated reaction paths for transitions under stress and finite deformation due to an inaccurate evaluation of the external work contributions in enthalpies. In this paper, a finite deformation nudged elastic band (FD-NEB) method is formulated for finding transition paths of solids under finite deformation. Applications of FD-NEB to a phase transition of silicon from the diamond phase to the β-tin phase under uniaxial compression are presented. The results are compared with those from the generalized solid-state nudged elastic band method. [ABSTRACT FROM AUTHOR]

This paper is the first application of the tensor-train (TT) cross approximation procedure for potential energy surface fitting. In order to reduce the complexity, we combine the TT-approach with another technique recently introduced in the field of numerical analysis: an affine transformation of Cartesian coordinates into the active subspaces where the PES function has the most variability. The numerical experiments for the water molecule and for the nitrous acid molecule confirm the efficiency of this approach. [ABSTRACT FROM AUTHOR]

This paper determines the thermodynamic phase diagram of the EXP system of particles interacting by the purely repulsive exponential pair potential. The solid phase is face-centered cubic (fcc) at low densities and pressures. At higher densities and pressures, the solid phase is body-centered cubic (bcc) with a re-entrant liquid phase at the highest pressures simulated. The investigation first identifies the phase diagram at zero temperature at which the following four crystal structures are considered: fcc, bcc, hexagonal close packed, and cubic diamond. There is a T = 0 phase transition at pressure 2.651 × 10−3 with the thermodynamically stable structure being fcc below and bcc above this pressure. The densities of the two crystal structures at the phase transition are 1.7469 × 10−2 (fcc) and 1.7471 × 10−2 (bcc). At finite temperatures, the fcc–bcc, fcc-liquid, and bcc-liquid coexistence lines are determined by numerical integration of the Clausius–Clapeyron equation and validated by interface-pinning simulations at selected state points. The bcc-fcc phase transition is a weak first-order transition. The liquid-fcc–bcc triple point, which is determined by the interface-pinning method, has temperature 5.9 × 10−5 and pressure 2.5 × 10−6; the triple-point densities are 1.556 × 10−3 (liquid), 1.583 × 10−3 (bcc), and 1.587 × 10−3 (fcc). [ABSTRACT FROM AUTHOR]

We demonstrate the utility of an unsupervised machine learning tool for the detection of phase transitions in off-lattice systems. We focus on the application of principal component analysis (PCA) to detect the freezing transitions of two-dimensional hard-disk and three-dimensional hard-sphere systems as well as liquid-gas phase separation in a patchy colloid model. As we demonstrate, PCA autonomously discovers order-parameter-like quantities that report on phase transitions, mitigating the need for a priori construction or identification of a suitable order parameter—thus streamlining the routine analysis of phase behavior. In a companion paper, we further develop the method established here to explore the detection of phase transitions in various model systems controlled by compositional demixing, liquid crystalline ordering, and non-equilibrium active forces. [ABSTRACT FROM AUTHOR]

The solid-solid and melting transitions that occur in Lennard-Jones LJn clusters have been both fascinating and challenging for the computational physics community over the last several decades. A number of attempts to extend these studies to the quantum case have also been made. Particularly interesting is the exploration of the parallel between the thermally induced and quantum-induced transitions. Yet, both numerically accurate and systematic studies of the latter are still lacking. In this paper, we apply the diffusion Monte Carlo method to the especially difficult case of LJ38. Starting with the truncated octahedral global minimum configuration, as the de Boer quantum delocalization parameterΛincreases, the system undergoes two consecutive solid-solid transitions, switching to anti-Mackay configurations. At sufficiently large values of Λ, the cluster is completely "melted" which is manifested by delocalization of the ground statewavefunction over a very large number of minima that represent several structural motifs. [ABSTRACT FROM AUTHOR]

Nature employs lipids to construct nanostructured membranes that self-assemble in an aqueous environment to separate the cell interior from the exterior environment. Membrane composition changes among species and according to environmental conditions, which allows organisms to occupy a wide variety of different habitats. Lipid bilayers are phase-change materials that exhibit strong thermotropic and lyotropic phase behavior in an aqueous environment, which may also cause thermal rectification. Among different types of lipids, archaeal lipids are of great interest due to their ability to withstand extreme conditions. In this paper, nonequilibrium molecular dynamics simulations were employed to study the nanostructures and thermal properties of different archaeols and to investigate thermal rectification effects in asymmetric archaeal membranes. In particular, we are interested in understanding the role of bridged phytanyl chains and cyclopentane groups in controlling the phase transition temperature and heat flow across the membrane. Our results indicate that the bridged phytanyl chains decrease the molecular packing of lipids, whereas the existence of cyclopentane rings on the tail groups increases the molecular packing by enhancing the interactions between isoprenoid chains. We found that macrocyclic archaeols have the highest thermal conductivity, whereas macrocyclic archaeols with two cyclopentane rings have the lowest. The effect of the temperature on the variation of thermal conductivity was found to be progressive. Our results further indicate that small thermal rectification effects occur in asymmetric archaeol bilayer membranes at around 25 K temperature gradient. The calculated thermal rectification factor was around 0.09 which is in the range of rectification factor obtained experimentally for nanostructures such as carbon nanotubes (0.07). Such phenomena may be of biological significance and could also be optimized for use in various engineering applications. [ABSTRACT FROM AUTHOR]

This paper is broadly concerned with the dynamics of a polymer confined to a rectangular slit of width D and deformed by a planar elongational flow of strength ẏ. It is interested, more specifically, in the nature of the coil-stretch transition that such polymers undergo when the flow strength ẏ is varied, and in the degree to which this transition is affected by the presence of restrictive boundaries. These issues are explored within the framework of a finitely extensible Rouse model that includes pre-averaged surface-mediated hydrodynamic interactions. Calculations of the chain's steady-state fractional extension x using this model suggest that different modes of relaxation (which are characterized by an integer p) exert different levels of control on the coil-stretch transition. In particular, the location of the transition (as identified from the graph of x versus the Weissenberg number Wi, a dimensionless parameter defined by the product of ẏ and the time constant τp of a relaxation mode p) is found to vary with the choice of τp. In particular, when τ1 is used in the definition of Wi, the x vs. Wi data for different D lie on a single curve, but when τ3 is used instead (with τ3 > τ1) the corresponding data lie on distinct curves. These findings are in close qualitative agreement with a number of experimental results on confinement effects on DNA stretching in electric fields. Similar D-dependent trends are seen in our calculated force vs. Wi data, but force vs. x data are essentially D-independent and lie on a single curve. [ABSTRACT FROM AUTHOR]

The laser induced fluorescence (LIF) spectrum of the à ²Δ - X ²Π transition was obtained for SiCN generated by laser ablation under supersonic free jet expansion. The vibrational structures of the dispersed fluorescence (DF) spectra from single vibronic levels (SVL's) were analyzed with consideration of the Renner-Teller (R-T) interaction. Analysis of the pure bending (ν2) structure by a perturbation approach including R-T, anharmonicity, spin-orbit (SO), and Herzberg-Teller (H-T) interactions indicated considerably different spin splitting for the μ and ? levels of the X ²Π state of SiCN, in contrast to identical spin splitting for general species derived from the perturbation approach, where μ and ? specify the lower and upper levels, respectively, separated by R-T. Further analysis of the vibrational structure including R-T, anharmonicity, SO, H-T, Fermi, and Sears interactions was carried out via a direct diagonalization procedure, where Sears resonance is a second-order interaction combined from SO and H-T interactions with ΔK = ±1, Δ∑ = ±1, and ΔP = 0, and where P is a quantum number, P = K + ∑. The later numerical analysis reproduced the observed structure, not only the pure ν2 structure but also the combination structure of the ν2 and the Si-CN stretching (ν3) modes. As an example, the analysis demonstrates Sears resonance between vibronic levels, (01¹0) ? ∑(+) and (0200) μ Π½, with ΔK = ±1 and ΔP = 0. On the basis of coefficients of their eigen vectors derived from the numerical analysis, it is interpreted as an almost one-to-one mixing between the two levels. The mixing coefficients of the two vibronic levels agree with those obtained from computational studies. The numerical analysis also indicates that some of the vibronic levels show chaotic characteristics in view of the two-dimensional harmonic oscillator (2D-HO) basis which is used as the basis function in the present numerical analysis; i.e., the eigen vectors for some of the observed levels have several components of the basis, and we have not been able to give precise vibronic assignments for the levels, but just vibronically labeled, referring the largest component in their vectors. (To emphasize this situation, we do not use the word "assignment," but prefer to use "label" as the meaning of just "label," but not "assign," throughout this paper.) The latter shows that the vibronic labels of the levels are meaningless, and the P quantum number and the order of their eigen states in the P matrix block derived in the numerical analysis only characterize the vibronic levels. Comparing the constants obtained for all of the interactions with those of species showing Sears resonance and studied previously, it is found that none of them are strong, but are moderate. It is thus concluded that the chaotic appearance is not derived by any strong interaction, but is induced by complex and accidental proximities of the vibronic levels caused by the moderate interactions. [ABSTRACT FROM AUTHOR]

When a CHD3 molecule is pumped to the C-H stretching-excited state by absorbing a linearly polarized infrared (IR) photon via the R(0) branch of the v1 = 1 ← 0 transition, the rotational angular momentum j of the prepared state |jK = |10 preferentially lies in a plane perpendicular to the IR polarization axis εIR. By way of contrast, when the Q(1) branch is used, the state of |jK = |1 ± 1 is prepared with j aligned along the direction of εIR. Reported here is a detailed study of the title reaction by actively controlling the collision geometries under these two IR-excitation schemes at collision energy Ec = 8.6 kcal mol−1, using a crossed molecular beam, product imaging approach. We found that under the R(0) excitation, the polarization-dependent differential cross sections for the HCl(v = 0) + CD3(00) channel can largely be understood by invoking dual reaction mechanisms. The forward-scattered products are most likely mediated by a time-delayed resonance mechanism--as the formation of the HCl(v = 1) + CD3(00) channel, whereas the backward/sideways scattered products are governed by a direct abstraction mechanism. Compared to the previous results at lower Ec of 3.8 kcal mol−1, the sighting of opening-up the attack angle at the transition state of the direct pathway is proposed. Results under the Q(1) excitation are, however, perplexing and bear no obvious correlation to the corresponding ones for the R(0) excitation, defying simple intuitive interpretation. Possible reasons are put forward, which call for theoretical investigations for deeper insights. The results on the alternative isotope channel, DCl + CHD2, will be reported in the following paper. [ABSTRACT FROM AUTHOR]

Restricting the number of attractive physical "bonds" that can form between particles in a fluid suppresses the usual demixing phase transition to very low particle concentrations, allowing for the formation of open, percolated, and homogeneous states, aptly called equilibrium or "empty" gels. Most demonstrations of this concept have directly limited the microscopic particle valence via anisotropic (patchy) attractions; however, an alternative macroscopic valence limitation would be desirable for greater experimental tunability and responsiveness. One possibility, explored in this paper, is to employ primary particles with attractions mediated via a secondary species of linking particles. In such a system, the linker-to-primary particle ratio serves as a macroscopic control parameter for the average microscopic valence. We show that the phase behavior of such a system predicted by Wertheim's first order perturbation theory is consistent with equilibrium gel formation: the primary particle concentrations corresponding to the two-phase demixing transition are significantly suppressed at both low and high linker-to-primary particle ratios. Extensive molecular dynamics simulations validate these theoretical predictions but also reveal the presence of loops of bonded particles, which are neglected in the theory. Such loops cause densification and inhibit percolation, and hence the range of viable empty gel state conditions is somewhat reduced relative to theWertheim theory predictions. [ABSTRACT FROM AUTHOR]

The Mössbauer effect was used to study the structural transitions in a PbZr0.72Sn0.28O3 single crystal. Two kinds of quadrupole splittings were registered and connected with two different environments of the Sn ion occupying the center of SnO6 octahedra. What is responsible for the existence of these two environments is the disorder in tilts of oxygen octahedra and antiparallel shifts of pairs of lead ions. Both disorders decide on structural transitions which the PbZr0.72Sn0.28O3 single crystal undergoes. The two kinds of quadruple splittings that have been observed do not disappear at temperatures far above phase transitions. This indicates that the structure of the paraelectric phase is locally non-centrosymmetric and confirms pre-transitional effects previously reported for Sn doped PbZrO3 single crystals. [ABSTRACT FROM AUTHOR]

This paper is devoted to bulk-mediated surface diffusion of a particle which can diffuse both on a flat surface and in the bulk layer above the surface. It is assumed that the particle is on the surface initially (at t = 0) and at time t, while in between it may escape from the surface and come back any number of times. We propose a new approach to the problem, which reduces its solution to that of a two-state problem of the particle transitions between the surface and the bulk layer, focusing on the cumulative residence times spent by the particle in the two states. These times are random variables, the sum of which is equal to the total observation time t. The advantage of the proposed approach is that it allows for a simple exact analytical solution for the double Laplace transform of the conditional probability density of the cumulative residence time spent on the surface by the particle observed for time t. This solution is used to find the Laplace transform of the particle mean square displacement and to analyze the peculiarities of its time behavior over the entire range of time. We also establish a relation between the double Laplace transform of the conditional probability density and the Fourier-Laplace transform of the particle propagator over the surface. The proposed approach treats the cases of both finite and infinite bulk layer thicknesses (where bulk-mediated surface diffusion is normal and anomalous at asymptotically long times, respectively) on equal footing. [ABSTRACT FROM AUTHOR]

The first-order liquid-liquid phase transition in supercooled Si is revisited by long-time first-principle molecular dynamics simulations. As the focus of the present paper, its nature is revealed by analyzing the inherent structures of low-density liquid (LDL) and high-density liquid (HDL). Our results show that it is a transition between a sp³-hybridization LDL and a white-tin-like HDL. This uncovers the origin of the semimetal-metal transition accompanying it and also proves that HDL is the metastable extension of high temperature equilibrium liquid into the supercooled regime. The pressuretemperature diagram of supercooled Si thus can be regarded in some respects as shifted reflection of its crystalline phase diagram. [ABSTRACT FROM AUTHOR]

In this paper, we present an improved method for obtaining unbiased estimates of the free energy difference between two thermodynamic states using the work distribution measured in nonequilibrium driven experiments connecting these states. The method is based on the assumption that any observed work distribution is given by a mixture of Gaussian distributions, whose normal components are identical in either direction of the nonequilibrium process, with weights regulated by the Crooks theorem. Using the prototypical example for the driven unfolding/folding of deca-alanine, we show that the predicted behavior of the forward and reverse work distributions, assuming a combination of only two Gaussian components with Crooks derived weights, explains surprisingly well the striking asymmetry in the observed distributions at fast pulling speeds. The proposed methodology opens the way for a perfectly parallel implementation of Jarzynski-based free energy calculations in complex systems. [ABSTRACT FROM AUTHOR]

This paper explores the driving mechanisms for structural transitions of atomic clusters between oblate and prolate isomers. We employ the hyperspherical coordinates to investigate structural dynamics of a seven-atom cluster at a coarse-grained level in terms of the dynamics of three gyration radii and three principal axes, which characterize overall mass distributions of the cluster. Dynamics of gyration radii is governed by two kinds of forces. One is the potential force originating from the interactions between atoms. The other is the dynamical forces called the internal centrifugal forces, which originate from twisting and shearing motions of the system. The internal centrifugal force arising from twisting motions has an effect of breaking the symmetry between two gyration radii. As a result, in an oblate isomer, activation of the internal centrifugal force that has the effect of breaking the symmetry between the two largest gyration radii is crucial in triggering structural transitions into prolate isomers. In a prolate isomer, on the other hand, activation of the internal centrifugal force that has the effect of breaking the symmetry between the two smallest gyration radii is crucial in triggering structural transitions into oblate isomers. Activation of a twisting motion that switches the movement patterns of three principal axes is also important for the onset of structural transitions between oblate and prolate isomers. Based on these trigger mechanisms, we finally show that selective activations of specific gyration radii and twisting motions, depending on the isomer of the cluster, can effectively induce structural transitions of the cluster. The results presented here could provide further insights into the control of molecular reactions. [ABSTRACT FROM AUTHOR]

Oxygen transport is central to many properties of oxides such as stoichiometric changes, phase transformation, and ionic conductivity. In this paper, we report a mechanism for oxygen transport in uranium dioxide (UO2) in which the kinetics is mediated by defect clustering dynamics. In particular, the kinetic Monte Carlo method has been used to investigate the kinetics of oxygen transport in UO2 under the condition of creation and annihilation of oxygen vacancies and interstitials as well as oxygen interstitial clustering, with variable off-stoichiometry and temperature conditions. It is found that in hypo-stoichiometric U O2-x, oxygen transport is well described by the vacancy diffusion mechanism while in hyper-stoichiometric UO2+x, oxygen interstitial cluster diffusion contributes significantly to oxygen transport kinetics, particularly at high temperatures and high off-stoichiometry levels. It is also found that di-interstitial clusters and single interstitials play dominant roles in oxygen diffusion while other larger clusters have negligible contributions. However, the formation, coalescence, and dissociation of these larger clusters indirectly affects the overall oxygen diffusion due to their interactions with mono and di-interstitials, thus providing an explanation of the experimental observation of saturation or even drop of oxygen diffusivity at high off-stoichiometry. [ABSTRACT FROM AUTHOR]

Hybrid systems consisting of organic molecules at inorganic semiconductor surfaces are gaining increasing importance as thin film devices for optoelectronics. The efficiency of such devices strongly depends on the collective behavior of the adsorbed molecules. In the present paper, we propose a novel, coarse-grained model addressing the condensed phases of a representative hybrid system, that is, para-sexiphenyl (6P) at zinc-oxide (ZnO). Within our model, intermolecular interactions are represented via a Gay-Berne potential (describing steric and van-der-Waals interactions) combined with the electrostatic potential between two linear quadrupoles. Similarly, the molecule-substrate interactions include a coupling between a linear molecular quadrupole to the electric field generated by the line charges characterizing ZnO(10-10). To validate our approach, we perform equilibrium Monte Carlo simulations, where the lateral positions are fixed to a 2D lattice, while the rotational degrees of freedom are continuous. We use these simulations to investigate orientational ordering in the condensed state. We reproduce various experimentally observed features such as the alignment of individual molecules with the line charges on the surface, the formation of a standing uniaxial phase with a herringbone structure, as well as the formation of a lying nematic phase. [ABSTRACT FROM AUTHOR]

We have updated the isotopically invariant Dunham fit of O2 with newly reported literature transitions to derive (1) the energy levels, partition sums, band-by-band molecular constants, and RKR poten-tials for the X³Σg-, a¹Δg+, and b¹Σg+ states of the six O2 isotopologues: 16O16O, 16O17O, 16O18O, 17O17O, 17O18O, and 18O18O; (2) Franck-Condon factors for their a¹Δg- X³Σg-, b¹Σg+ - X³Σg-, and a¹Δg -- b¹Σg+ band systems. This new spectroscopic parameterization characterizes all known transitions within and between the X³Σg--, a¹Δg, and b¹Σg+ states within experimental uncertainty and can be used for accurate predictions of as yet unmeasured transitions. All of these results are necessary to provide a consistent linelist of all transitions which will be reported in a followup paper. [ABSTRACT FROM AUTHOR]

The theory of optical transitions developed in Barford and Marcus ["Theory of optical transitions in conjugated polymers. I. Ideal systems," J. Chem. Phys. 141, 164101 (2014)] for linear, ordered polymer chains is extended in this paper to model conformationally disordered systems. Our key result is that in the Born-Oppenheimer regime the emission intensities are proportional to S(1)/IPR, where S(1) is the Huang-Rhys parameter for a monomer. IPR is the average inverse participation ratio for the emitting species, i.e., local exciton ground states (LEGSs). Since the spatial coherence of LEGSs determines the spatial extent of chromophores, the significance of this result is that it directly relates experimental observables to chromophore sizes (where IPR is half the mean chromophore size in monomer units). This result is independent of the chromophore shape, because of the Born-Oppenheimer factorization of the many body wavefunction. We verify this prediction by density matrix renormalization group (DMRG) calculations of the Frenkel-Holstein model in the adiabatic limit for both linear, disordered chains and for coiled, ordered chains. We also model optical spectra for poly(p-phenylene) and poly(p-phenylene-vinylene) oligomers and polymers. For oligomers, we solve the fully quantized Frenkel-Holstein model via the DMRG method. For polymers, we use the much simpler method of solving the one-particle Frenkel model and employ the Born-Oppenheimer expressions relating the effective Franck-Condon factor of a chromophore to its inverse participation ratio. We show that increased disorder decreases chromophore sizes and increases the inhomogeneous broadening, but has a non-monotonic effect on transition energies. We also show that as planarizing the polymer chain increases the exciton band width, it causes the chromophore sizes to increase, the transition energies to decrease, and the broadening to decrease. Finally, we show that the absorption spectra are more broadened than the emission spectra and that the broadening of the absorption spectra increases as the chains become more coiled. This is primarily because absorption occurs to both LEGSs and quasi-extended exciton states (QEESs), and QEES acquire increased intensity as chromophores bend, while emission only occurs from LEGSs. [ABSTRACT FROM AUTHOR]

Commensurate-incommensurate (C-IC) transition of krypton molecular layer on graphite received much attention in recent decades in theoretical and experimental researches. However, there still exists a possibility of generalization of the phenomenon from thermodynamic viewpoint on the basis of accurate molecular simulation. Recently, a new technique was developed for analysis of two-dimensional (2D) phase transitions in systems involving a crystalline phase, which is based on accounting for the effect of temperature and the chemical potential on the lattice constant of the 2D layer using the Gibbs-Duhem equation [E. A. Ustinov, J. Chem. Phys. 140, 074706 (2014)]. The technique has allowed for determination of phase diagrams of 2D argon layers on the uniform surface and in slit pores. This paper extends the developed methodology on systems accounting for the periodic modulation of the substrate potential. The main advantage of the developed approach is that it provides highly accurate evaluation of the chemical potential of crystalline layers, which allows reliable determination of temperature and other parameters of various 2D phase transitions. Applicability of the methodology is demonstrated on the kryptongraphite system. Analysis of phase diagram of the krypton molecular layer, thermodynamic functions of coexisting phases, and a method of prediction of adsorption isotherms is considered accounting for a compression of the graphite due to the krypton-carbon interaction. The temperature and heat of C-IC transition has been reliably determined for the gas-solid and solid-solid system. [ABSTRACT FROM AUTHOR]

The purpose of this paper is to introduce a quasiparticle-based multi-reference coupled-cluster (MRCC) approach. The quasiparticles are introduced via a unitary transformation which allows us to represent a complete active space reference function and other elements of an orthonormal multi-reference (MR) basis in a determinant-like form. The quasiparticle creation and annihilation operators satisfy the fermion anti-commutation relations. On the basis of these quasiparticles, a generalization of the normal-ordered operator products for the MR case can be introduced as an alternative to the approach of Mukherjee and Kutzelnigg [Recent Prog. Many-Body Theor. 4, 127 (1995); J. Chem. Phys. 107, 432 (1997)]. Based on the new normal ordering any quasiparticle-based theory can be formulated using the well-known diagram techniques. Beyond the general quasiparticle framework we also present a possible realization of the unitary transformation. The suggested transformation has an exponential form where the parameters, holding exclusively active indices, are defined in a form similar to the wave operator of the unitary coupled-cluster approach. The definition of our quasiparticle-based MRCC approach strictly follows the form of the single-reference coupled-cluster method and retains several of its beneficial properties. Test results for small systems are presented using a pilot implementation of the new approach and compared to those obtained by other MR methods. [ABSTRACT FROM AUTHOR]

Franck-Condon vibrational overlap integrals for the Ã1Au—X̃1Σg+> transition in acetylene have been calculated in full dimension in the harmonic normal mode basis. The calculation uses the method of generating functions first developed for polyatomic Franck-Condon factors by Sharp and Rosenstock [J. Chem. Phys. 41(11), 3453-3463 (1964)], and previously applied to acetylene by Watson [J. Mol. Spectrosc. 207(2), 276-284 (2001)] in a reduced-dimension calculation. Because the transition involves a large change in the equilibrium geometry of the electronic states, two different types of corrections to the coordinate transformation are considered to first order: corrections for axis-switching between the Cartesian molecular frames and corrections for the curvilinear nature of the normal modes at large amplitude. The angular factor in the wavefunction for the out-of-plane component of the trans bending mode, v"4, is treated as a rotation, which results in an Eckart constraint on the polar coordinates of the bending modes. To simplify the calculation, the other degenerate bending mode, v"5, is integrated in the Cartesian basis and later transformed to the constrained polar coordinate basis, restoring the conventional v and l quantum numbers. An updated Ã-state harmonic force field obtained recently in the R. W. Field research group is evaluated. The results for transitions involving the gerade vibrational modes are in qualitative agreement with experiment. Calculated results for transitions involving ungerade modes are presented in Paper II of this series [G. B. Park, J. H. Baraban, and R. W. Field, "Full dimensional Franck-Condon factors for the acetylene Ã1Au—X̃1Σg+ transition. II. Vibrational overlap factors for levels involving excitation in ungerade modes," J. Chem. Phys. 141, 134305 (2014)]. [ABSTRACT FROM AUTHOR]

Although the volume transition of the polyelectrolyte gel has been studied for decades, less study on the finite size effects of the mobile ions has been conducted. In the present paper, Tanaka's classical theory of polyelectrolyte gel is extended to the cases of mobile ions of finite volume. In the salt free limit, the theoretical results show that the discontinuous volume transition of the polyelectrolyte gel will become a continuous one for counterionic size larger than a critical value. When a significant amount of salt is added, the critical value for the volume transition increases as a result of electrostatic screening. An increase in salt concentration can also make the polyelectrolyte gel in poor solvent collapse. Poorer solvent is needed to trigger the salt-induced collapse in polyelectrolyte gel with larger mobile ions than that with smaller ones. The effects of ionic size on the critical points and phase diagram of the volume transition are also discussed. The theoretical results suggest that the swelling behavior of polyelectrolyte gel might be tuned with salt of different volumes. [ABSTRACT FROM AUTHOR]

We outline how principal component analysis can be applied to particle configuration data to detect a variety of phase transitions in off-lattice systems, both in and out of equilibrium. Specifically, we discuss its application to study (1) the nonequilibrium random organization (RandOrg) model that exhibits a phase transition from quiescent to steady-state behavior as a function of density, (2) orientationally and positionally driven equilibrium phase transitions for hard ellipses, and (3) a compositionally driven demixing transition in the non-additive binary Widom-Rowlinson mixture. [ABSTRACT FROM AUTHOR]

We report the structure of transient fluctuations in the liquid phase of a two-dimensional system that exhibits several ordered phases with different symmetries. The density-temperature phase diagram of the system studied, composed of particles with a repulsive shouldered soft-core pair interaction, has regions with stable liquid and hexatic phases, a square solid phase, two separate hexagonal solid phases, and a quasi-crystalline phase with 12-fold symmetry. We have examined the character of the structured fluctuations by computing the same-time aperture cross correlation function of particle configurations in several fluid regions near to and far from phase transition lines. The two primary goals of our study are (1) determination if the spectrum of structures of the fluctuations in the liquid is broader than or limited to the motifs exhibited by the ordered phases supported by the system and (2) determination of the density domains in the liquid that support particular transient structured fluctuations. In the system studied, along a low-temperature isotherm in the temperature-density plane that intersects all the ordered phases we find that the liquid phase exhibits structured fluctuations with hexagonal symmetry near both liquid-hexatic transition lines. Along the same isotherm and in the stable liquid between the lower density hexatic-to-liquid and the higher density liquid-to-square solid transitions, we find that transient hexagonal ordered fluctuations dominate the liquid region near the hexatic-to-liquid transition and square ordered fluctuations dominate the liquid region near the liquid-to square solid transition, but both of these structured fluctuations occur at all densities between these transition lines. At a higher temperature, at phase points in the liquid above, but close to the density maximum of an underlying transition, there are ordered fluctuations that can be correlated with the structure of the lower temperature phase. Although it is expected that very close to a liquid-ordered phase boundary a structured fluctuation in the liquid will have the same symmetry as the ordered phase, it is not obvious that structured fluctuations in thermodynamic states deep in the liquid phase will be similarly restricted. The most striking result of our calculations is that no evidence is found in the liquid phase for structured fluctuations with other symmetries than those of the ordered phases of the system. [ABSTRACT FROM AUTHOR]

Crystallization of amorphous materials by thermal annealing has been investigated for numerous applications in the fields of nanotechnology, such as thin-film transistors and thermoelectric devices. The phase transition and shape evolution of amorphous germanium (Ge) and Ag@Ge core–shell nanoparticles with average diameters of 10 and 12 nm, respectively, were investigated by high-energy electron beam irradiation and in situ heating within a transmission electron microscope. The transition of a single Ge amorphous nanoparticle to the crystalline diamond cubic structure at the atomic scale was clearly demonstrated. Depending on the heating temperature, a hollow Ge structure can be maintained or transformed into a solid Ge nanocrystal through a diffusive process during the amorphous to crystalline phase transition. Selected area diffraction patterns were obtained to confirm the crystallization process. In addition, the thermal stability of Ag@Ge core–shell nanoparticles with an average core of 7.4 and a 2.1 nm Ge shell was studied by applying the same beam conditions and temperatures. The results show that at a moderate temperature (e.g., 385 °C), the amorphous Ge shell can completely crystallize while maintaining the well-defined core–shell structure, while at a high temperature (e.g., 545 °C), the high thermal energy enables a freely diffusive process of both Ag and Ge atoms on the carbon support film and leads to transformation into a phase segregated Ag–Ge Janus nanoparticle with a clear interface between the Ag and Ge domains. This study provides a protocol as well as insight into the thermal stability and strain relief mechanism of complex nanostructures at the single nanoparticle level with atomic resolution. [ABSTRACT FROM AUTHOR]

As an important phase-change material, GeTe has many high-pressure phases as well, but its phase transitions under pressure are still lack of clarity. It is challenging to identify high-pressure GeTe crystal structures owing to the phase coexistence in a wide pressure range and the reversibility of phase transitions. Hence, first-principles calculations are required to provide further information in addition to limited experimental characterizations. In this work, a new orthorhombic Cmca GeTe high-pressure phase has been predicted via the CALYPSO method as the most energetically favorable phase in the pressure range between ∼30 and ∼38.5 GPa, which would update the GeTe high-pressure phase transition sequence. The crystal structure of the Cmca phase is composed of alternate stacking puckered layers of Ge six-membered rings and Te four-membered rings along the b direction. The high density of states near the Fermi level and delocalization of electrons from the two-dimensional electron localization function indicate a strong metallic property of the Cmca phase. Electron–phonon coupling calculations indicate that the Cmca phase is superconductive below ∼4.2 K at 35 GPa. The simulated x-ray diffraction pattern of the Cmca phase implies that this phase might coexist with the Pnma-boat phase under high pressure. These results offer further understanding on the high-pressure structural evolution and physical properties in GeTe and other IV–VI semiconductors. [ABSTRACT FROM AUTHOR]

The most accepted origin for the water anomalous behavior is the phase transition between two liquids (LLPT) in the supercooled regime connected to the glassy first order phase transition at lower temperatures. Two length scale potentials are an effective approach that has long been employed to understand the properties of fluids with waterlike anomalies and, more recently, the behavior of colloids and nanoparticles. These potentials can be parameterized to have distinct shapes, as a pure repulsive ramp, such as the model proposed by de Oliveira et al. [J. Chem. Phys. 124, 64901 (2006)]. This model has waterlike anomalies despite the absence of LLPT. To unravel how the waterlike anomalies are connected to the solid phases, we employ molecular dynamics simulations. We have analyzed the fluid–solid transition under cooling, with two solid crystalline phases, BCC and HCP, and two amorphous regions being observed. We show how the competition between the scales creates an amorphous cluster in the BCC crystal that leads to amorphization at low temperatures. A similar mechanism is found in the fluid phase, with the system changing from a BCC-like to an amorphous-like structure in the point where a maxima in kT is observed. With this, we can relate the competition between two fluid structures with the amorphous clusterization in the BCC phase. These findings help to understand the origins of waterlike behavior in systems without the liquid–liquid critical point. [ABSTRACT FROM AUTHOR]

Halide perovskites are promising light-absorbing materials for high-efficiency solar cells, while the crystalline phase of halide perovskites may influence the device's efficiency and stability. In this work, we investigated the thermally driven phase transition of perovskite (CsPbIxBr3−x), which was confirmed by electron diffraction and high-resolution transmission electron microscopy results. CsPbIxBr3−x transitioned from δ phase to α phase when heated, and the γ phase was obtained when the sample was cooled down. The γ phase was stable as long as it was isolated from humidity and air. A template matching-based data analysis method enabled visualization of the thermally driven phase evolution of perovskite during heating. We also proposed a possible atomic movement in the process of phase transition based on our in situ heating experimental data. The results presented here may improve our understanding of the thermally driven phase transition of perovskite as well as provide a protocol for big-data analysis of in situ experiments. [ABSTRACT FROM AUTHOR]

Allostery in proteins involves, broadly speaking, ligand-induced conformational transitions that modulate function at active sites distal to where the ligand binds. In contrast, the concept of cooperativity (in the sense used in phase transition theory) is often invoked to understand protein folding and, therefore, function. The modern view on allostery is one based on dynamics and hinges on the time-dependent interactions between key residues in a complex network, interactions that determine the free-energy profile for the reaction at the distal site. Here, we merge allostery and cooperativity, and we discuss a joint model with features of both. In our model, the active-site reaction is replaced by the reaction pathway that leads to protein folding, and the presence or absence of the effector is replaced by mutant-vs-wild type changes in key residues. To this end, we employ our recently introduced time-lagged independent component analysis (tICA) correlation approach [Ray et al. Proc. Natl. Acad. Sci. 118(43) (2021), e2100943118] to identify the allosteric role of distant residues in the folded-state dynamics of a large protein. In this work, we apply the technique to identify key residues that have a significant role in the folding of a small, fast folding-protein, chignolin. Using extensive enhanced sampling simulations, we critically evaluate the accuracy of the predictions by mutating each residue one at a time and studying how the mutations change the underlying free energy landscape of the folding process. We observe that mutations in those residues whose associated backbone torsion angles have a high correlation score can indeed lead to loss of stability of the folded configuration. We also provide a rationale based on interaction energies between individual residues with the rest of the protein to explain this effect. From these observations, we conclude that the tICA correlation score metric is a useful tool for predicting the role of individual residues in the correlated dynamics of proteins and can find application to the problem of identifying regions of protein that are either most vulnerable to mutations or—mutatis mutandis—to binding events that affect their functionality. [ABSTRACT FROM AUTHOR]

We study low-salt, binary aqueous suspensions of charged colloidal spheres of size ratio Γ = 0.57, number densities below the eutectic number density nE, and number fractions of p = 1.00–0.40. The typical phase obtained by solidification from a homogeneous shear-melt is a substitutional alloy with a body centered cubic structure. In strictly gas-tight vials, the polycrystalline solid is stable against melting and further phase transformation for extended times. For comparison, we also prepare the same samples by slow, mechanically undisturbed deionization in commercial slit cells. These cells feature a complex but well reproducible sequence of global and local gradients in salt concentration, number density, and composition as induced by successive deionization, phoretic transport, and differential settling of the components, respectively. Moreover, they provide an extended bottom surface suitable for heterogeneous nucleation of the β-phase. We give a detailed qualitative characterization of the crystallization processes using imaging and optical microscopy. By contrast to the bulk samples, the initial alloy formation is not volume-filling, and we now observe also α- and β-phases with low solubility of the odd component. In addition to the initial homogeneous nucleation route, the interplay of gradients opens various further crystallization and transformation pathways leading to a great diversity of microstructures. Upon a subsequent increase in salt concentration, the crystals melt again. Wall-based, pebble-shaped β-phase crystals and facetted α-crystals melt last. Our observations suggest that the substitutional alloys formed in bulk experiments by homogeneous nucleation and subsequent growth are mechanically stable in the absence of solid–fluid interfaces but thermodynamically metastable. [ABSTRACT FROM AUTHOR]

The slow relaxation of interface water (IW) across three primary phases of membranes is relevant to understand the influence of IW on membrane functions at supercooled conditions. To this objective, a total of ∼ 16.26 μ s all-atom molecular dynamics simulations of 1,2-dimyristoyl-sn-glycerol-3-phosphocholine lipid membranes are carried out. A supercooling-driven drastic slow-down in heterogeneity time scales of the IW is found at the fluid to the ripple to the gel phase transitions of the membranes. At both fluid-to-ripple-to-gel phase transitions, the IW undergoes two dynamic crossovers in Arrhenius behavior with the highest activation energy at the gel phase due to the highest number of hydrogen bonds. Interestingly, the Stokes–Einstein (SE) relation is conserved for the IW near all three phases of the membranes for the time scales derived from the diffusion exponents and the non-Gaussian parameters. However, the SE relation breaks for the time scale obtained from the self-intermediate scattering functions. The behavioral difference in different time scales is universal and found to be an intrinsic property of glass. The first dynamical transition in the α relaxation time of the IW is associated with an increase in the Gibbs energy of activation of hydrogen bond breaking with locally distorted tetrahedral structures, unlike the bulk water. Thus, our analyses unveil the nature of the relaxation time scales of the IW across membrane phase transitions in comparison with the bulk water. The results will be useful to understand the activities and survival of complex biomembranes under supercooled conditions in the future. [ABSTRACT FROM AUTHOR]

We report a molecular dynamics study of the heterogeneous crystallization of high-pressure glassy water using (plastic) ice VII as a substrate. We focus on the thermodynamic conditions P ∈ [6–8] GPa and T ∈ [100–500] K, at which (plastic) ice VII and glassy water are supposed to coexist in several (exo)planets and icy moons. We find that (plastic) ice VII undergoes a martensitic phase transition to a (plastic) fcc crystal. Depending on the molecular rotational lifetime τ, we identify three rotational regimes: for τ > 20 ps, crystallization does not occur; for τ ∼ 15 ps, we observe a very sluggish crystallization and the formation of a considerable amount of icosahedral environments trapped in a highly defective crystal or in the residual glassy matrix; and for τ < 10 ps, crystallization takes place smoothly, resulting in an almost defect-free plastic fcc solid. The presence of icosahedral environments at intermediate τ is of particular interest as it shows that such a geometry, otherwise ephemeral at lower pressures, is, indeed, present in water. We justify the presence of icosahedral structures based on geometrical arguments. Our results represent the first study of heterogeneous crystallization occurring at thermodynamic conditions of relevance for planetary science and unveil the role of molecular rotations in achieving it. Our findings (i) show that the stability of plastic ice VII, widely reported in the literature, should be reconsidered in favor of plastic fcc, (ii) provide a rationale for the role of molecular rotations in achieving heterogeneous crystallization, and (iii) represent the first evidence of long-living icosahedral structures in water. Therefore, our work pushes forward our understanding of the properties of water. [ABSTRACT FROM AUTHOR]

An ionic liquid (IL) is a salt in the liquid state that consists of a cation and an anion, one of which possesses an organic component. Because of their non-volatile property, these solvents have a high recovery rate, and, hence, they are considered as environment-friendly green solvents. It is necessary to study the detailed physicochemical properties of these liquids for designing and processing techniques and find suitable operating conditions for IL-based systems. In the present work, the flow behavior of aqueous solutions of an imidazolium-based IL, 1-methyl-3-octylimidazolium chloride, is investigated, where the dynamic viscosity measurements indicate non-Newtonian shear thickening behavior in the solutions. Polarizing optical microscopy shows that the pristine samples are isotropic and transform into anisotropic after shear. These shear thickened liquid crystalline samples change into an isotropic phase upon heating, which is quantified by the differential scanning calorimetry. The small angle x-ray scattering study revealed that the pristine isotropic cubic phase of spherical micelles distort into non-spherical micelles. This has provided the detailed structural evolution of mesoscopic aggregates of the IL in an aqueous solution and the corresponding viscoelastic property of the solution. [ABSTRACT FROM AUTHOR]

The interaction of water with different substances in the earth's atmosphere lies at the heart of many processes that influence our climate. However, it is still unclear how different species interact with water on the molecular level and in which ways this interaction contributes to the water vapor phase transition. Here, we report the first measurements of water–nonane binary nucleation in the 50–110 K temperature range, along with unary nucleation data of both. The time-dependent cluster size distribution in a uniform post-nozzle flow was measured by time-of-flight mass spectrometry coupled with single-photon ionization. From these data, we extract experimental rates and rate constants for both nucleation and cluster growth. The observed mass spectra of water/nonane clusters are not or only slightly affected by the introduction of the other vapor, and the formation of mixed clusters was not observed during nucleation of the mixed vapor. Additionally, the nucleation rate of either substance is not much affected by the presence (or absence) of the other species, i.e., the nucleation of water and nonane proceeds independently, indicating that hetero-molecular clusters do not play a role during nucleation. Only at the lowest temperature of our experiment (i.e., 51 K) do the measurements suggest that interspecies interaction slows water cluster growth. The findings here are in contrast to our earlier work in which we showed that vapor components in other mixtures, e.g., CO2 and toluene/H2O, can interact to promote nucleation and cluster growth in a similar temperature range. [ABSTRACT FROM AUTHOR]

We study the effect of particle mobility on phase transitions in a spin fluid in two dimensions. The presence of a phase transition of the BKT universality class is shown in an off-lattice model of particles with purely repulsive interaction employing computer simulations. A critical spin wave region 0 < T < TBKT is found with a nonuniversal exponent η(T) that follows the shape suggested by BKT theory, including a critical value consistent with ηBKT = 1/4. One can observe a transition from power-law decay to exponential decay in the static correlation functions at the transition temperature TBKT, which is supported by finite-size scaling analysis. A critical temperature TBKT = 0.17(1) is suggested. Investigations into the dynamic aspects of the phase transition are carried out. The short-time behavior of the incoherent spin autocorrelation function agrees with the Nelson–Fisher prediction, whereas the long-time behavior differs from the finite-size scaling known for the static XY model. Analysis of coherent spin wave dynamics shows that the spin wave peak is a propagating mode that can be reasonably well fitted by hydrodynamic theory. The mobility of the particles strongly enhances damping of the spin waves, but the model still lies within the dynamic universality class of the standard XY model. [ABSTRACT FROM AUTHOR]

Monte Carlo simulations were used to study the influence of particle aspect ratio on the kinetics and phase behavior of hard gyrobifastigia (GBF). First, the formation of a highly anisotropic nucleus shape in the isotropic-to-crystal transition in regular GBF is explained by the differences in interfacial free energies of various crystal planes and the nucleus geometry predicted by the Wulff construction. GBF-related shapes with various aspect ratios were then studied, mapping their equations of state, determining phase coexistence conditions via interfacial pinning, and computing nucleation free-energy barriers via umbrella sampling using suitable order parameters. Our simulations reveal a reduction of the kinetic barrier for isotropic–crystal transition upon an increase in aspect ratio, and that for highly oblate and prolate aspect ratios, an intermediate nematic phase is stabilized. Our results and observations also support two conjectures for the formation of the crystalline state from the isotropic phase: that low phase free energies at the ordering phase transition correlate with low transition barriers and that the emergence of a mesophase provides a steppingstone that expedites crystallization. [ABSTRACT FROM AUTHOR]

Simulations of water near extended hydrophobic spherical solutes have revealed the presence of a region of depleted density and accompanying enhanced density fluctuations. The physical origin of both phenomena has remained somewhat obscure. We investigate these effects employing a mesoscopic binding potential analysis, classical density functional theory (DFT) calculations for a simple Lennard-Jones solvent, and Grand Canonical Monte Carlo (GCMC) simulations of a monatomic water (mw) model. We argue that the density depletion and enhanced fluctuations are near-critical phenomena. Specifically, we show that they can be viewed as remnants of the critical drying surface phase transition that occurs at bulk liquid–vapor coexistence in the macroscopic planar limit, i.e., as the solute radius Rs → ∞. Focusing on the radial density profile ρ(r) and a sensitive spatial measure of fluctuations, the local compressibility profile χ(r), our binding potential analysis provides explicit predictions for the manner in which the key features of ρ(r) and χ(r) scale with Rs, the strength of solute–water attraction ɛsf, and the deviation from liquid–vapor coexistence of the chemical potential, δμ. These scaling predictions are confirmed by our DFT calculations and GCMC simulations. As such, our theory provides a firm basis for understanding the physics of hydrophobic solvation. [ABSTRACT FROM AUTHOR]

In systems with frustration, the critical slowing down of the dynamics severely impedes the numerical study of phase transitions for even the simplest of lattice models. In order to help sidestep the gelation-like sluggishness, a clearer understanding of the underlying physics is needed. Here, we first obtain generic insight into that phenomenon by studying one-dimensional and Bethe lattice versions of a schematic frustrated model, the axial next-nearest neighbor Ising (ANNNI) model. Based on these findings, we formulate two cluster algorithms that speed up the simulations of the ANNNI model on a 2D square lattice. Although these schemes do not eliminate the critical slowing own, speed-ups of factors up to 40 are achieved in some regimes. [ABSTRACT FROM AUTHOR]