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]

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]

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]

Amorphous solids, such as glasses and gels, arise as the asymptotic limit of non-equilibrium and irreversible relaxation aging processes. These amorphous solids form when the system is suddenly and deeply quenched in the dynamic arrest region. We use the non-equilibrium self-consistent generalized Langevin equation (NE-SCGLE) theory to investigate the formation of such structures via arrested spinodal decomposition in the screened symmetric restricted primitive model. We propose a direct correlation function that allows us to derive an expression for the functional derivative of the chemical potential, which serves as the necessary input in the NE-SCGLE theory. By analyzing the asymptotic localization length and the asymptotic characteristic length, we identify different classes of dynamic arrest states as a function of the density and the final quench temperature. The system features simultaneously attractive and repulsive interactions, resulting in different arrested regions in the non-equilibrium phase diagram for a given screened parameter: (i) ionic glasses, (ii) electrostatic gels, and (iii) attractive glasses. Finally, by collapsing the asymptotic effective structure factors at a typical fractal dimension of df = 3, we confirm the formation of gels slightly above the glass-gel transition in the arrested phase diagram. [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]

Mesogens, which are typically stiff rodlike or disklike molecules, are able to self-organize into liquid crystal (LC) phases in a certain temperature range. Such mesogens, or LC groups, can be attached to polymer chains in various configurations including within the backbone (main-chain LC polymers) or at the ends of side-chains attached to the backbone in an end-on or side-on configuration (side-chain LC polymers or SCLCPs), which can display synergistic properties arising from both their LC and polymeric character. At lower temperatures, chain conformations may be significantly altered due to the mesoscale LC ordering; thus, when heated from the LC ordered state through the LC to isotropic phase transition, the chains return from a more stretched to a more random coil conformation. This can cause macroscopic shape changes, which depend significantly on the type of LC attachment and other architectural properties of the polymer. Here, to study the structure–property relationships for SCLCPs with a range of different architectures, we develop a coarse-grained model that includes torsional potentials along with LC interactions of a Gay–Berne form. We create systems of different side-chain lengths, chain stiffnesses, and LC attachment types and track their structural properties as a function of temperature. Our modeled systems indeed form a variety of well-organized mesophase structures at low temperatures, and we predict higher LC-to-isotropic transition temperatures for the end-on side-chain systems than for analogous side-on side-chain systems. Understanding these phase transitions and their dependence on polymer architecture can be useful in designing materials with reversible and controllable deformations. [ABSTRACT FROM AUTHOR]

Many molecular systems and physical phenomena are controlled by local fluctuations and microscopic dynamical rearrangements of the constitutive interacting units that are often difficult to detect. This is the case, for example, of phase transitions, phase equilibria, nucleation events, and defect propagation, to mention a few. A detailed comprehension of local atomic environments and of their dynamic rearrangements is essential to understand such phenomena and also to draw structure–property relationships useful to unveil how to control complex molecular systems. Considerable progress in the development of advanced structural descriptors [e.g., Smooth Overlap of Atomic Position (SOAP), etc.] has certainly enhanced the representation of atomic-scale simulations data. However, despite such efforts, local dynamic environment rearrangements still remain difficult to elucidate. Here, exploiting the structurally rich description of atomic environments of SOAP and building on the concept of time-dependent local variations, we developed a SOAP-based descriptor, TimeSOAP (τSOAP), which essentially tracks time variations in local SOAP environments surrounding each molecule (i.e., each SOAP center) along ensemble trajectories. We demonstrate how analysis of the time-series τSOAP data and of their time derivatives allows us to detect dynamic domains and track instantaneous changes of local atomic arrangements (i.e., local fluctuations) in a variety of molecular systems. The approach is simple and general, and we expect that it will help shed light on a variety of complex dynamical phenomena. [ABSTRACT FROM AUTHOR]

The objective of this study was to investigate the prediction of the wetting characteristics obtained from the equilibrium adsorption analysis using the Zeta adsorption isotherm approach with an experimental study. Water vapor's adsorption and wetting characteristics on a hydroxylated and nano-polished silica substrate were studied in near-equilibrium conditions at temperatures near 298 K. Using a UV–visible interferometer, water vapor adsorbate film thicknesses were measured and converted into amount adsorbed per unit area. The current results show that the wetting transition occurred at an average subcooling value of 0.39 K, less than the predicted value of 0.49 K. All the different experimental observations showed growth of film thickness as a function of subcooling value with a maximum film thickness of 12.6 nm. The analysis of the results further showed that the maximum stable film was in a metastable state that then condensed in a dropwise manner, if perturbed by increasing the subcooling. The study further revealed that the adsorbate is unstable after transitioning. The solid surface energy calculated by including the near-equilibrium observations was comparable and close to that of the equilibrium studies, thus supporting solid surface energy as a material property. [ABSTRACT FROM AUTHOR]

Photoluminescence (PL) spectroscopy of individual semiconductor nanocrystals (NCs) is a powerful method for understanding the intrinsic optical properties of these materials. Here, we report the temperature dependence of the PL spectra of single perovskite FAPbBr3 and CsPbBr3 NCs [FA = HC(NH2)2]. The temperature dependences of the PL linewidths were mainly determined by the Fröhlich interaction between excitons and longitudinal optical phonons. For FAPbBr3 NCs, a redshift in the PL peak energy appeared between 100 and 150 K, which was because of the orthorhombic-to-tetragonal phase transition. We found that the phase transition temperature of FAPbBr3 NCs decreases with decreasing NC size. [ABSTRACT FROM AUTHOR]

Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology—a machine learning approach—to study this ubiquitous phase transition, starting from the phase diagram in the liquid–vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. We find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid–vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid–vapor coexistence and water cavitation. [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]

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]

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]

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]

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]

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]

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]

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]

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]

Organic ionic plastic crystals (OIPCs), which consist of organic molecular ions, are considered excellent candidates for solid electrolytes due to their high ionic conductivity in solid phases. Molecular ions undergo either rotational or conformational relaxation at certain temperatures in OIPCs. There have been molecular simulations to understand the rotational motion. The polarizability of ions was, however, often ignored in simulations due to the high computational cost. Since the polarizability may affect the translational diffusion, the ionic conductivity, and the phase transition of ionic liquids, it should be of interest to investigate how the polarizability would affect the rotational diffusion of ions in solid phases. In this work, we perform extensive atomistic molecular dynamics simulations for two different kinds of OIPCs, 1-methyl-3-methylimidazolium hexafluorophosphate ([MMIM][PF6]) and 1-butyl-3-methylimidazolium hexafluorophosphate ([BMIM][PF6]). We employ various simulation models for ions by turning on and off the polarization in their interaction potentials. We find that the polarizability hardly affects the density, the crystalline structure, and the phase transition of both OIPCs. However, a certain rotational motion, especially the rotational diffusion of PF 6 − in [MMIM][PF6] OIPCs, is enhanced by a factor of up to four when the polarizability is turned on. The PF 6 − in [MMIM][PF6] OIPCs undergoes rotational hopping motions more significantly due to polarizability. We find that the rotational diffusion of a certain ion can be enhanced only when the polarization results in a significant change in the dipole moment of the neighboring ions around the ion. [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 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]

The thermodynamics of hard spheres tethered to a Face-Centered Cubic (FCC) lattice is investigated using event-driven molecular-dynamics. The particle–particle and the particle–tether collision rates are related to the phase space geometry and are used to study the FCC and fluid states. In tethered systems, the entropy can be determined by at least two routes: (i) through integration of the tether collision rates with the tether length rT or (ii) through integration of the particle–particle collision rates with the hard-sphere diameter σ (or, equivalently, the density). If the entropy were an entirely analytic function of rT and σ, these two methods for calculating the entropy should lead to the same results; however, a non-analytic region exists as an extension of the solid–fluid phase transition of the untethered hard-sphere system, and integration paths that cross this region will lead to values for the entropy that depend on the particular path chosen. The difference between the calculated entropies appears to be related to the communal entropy, and the location of the non-analytic region appears to be related to conditions where the regions of phase space associated with the FCC configuration become separated from those associated with the disordered fluid. The non-analytic region is finite in extent, vanishing below rT/a ≈ 0.55, where a is the lattice spacing, and there are many continuous paths that connect the fluid and solid phases that can be used to determine the crystal free energy with respect to the fluid. [ABSTRACT FROM AUTHOR]

Freezing of water is the most common liquid-to-crystal phase transition on Earth; however, despite its critical implications on climate change and cryopreservation among other disciplines, its characterization through experimental and computational techniques remains elusive. In this work, we make use of computer simulations to measure the nucleation rate (J) of water at normal pressure under different supercooling conditions, ranging from 215 to 240 K. We employ two different water models: mW, a coarse-grained potential for water, and TIP4P/ICE, an atomistic nonpolarizable water model that provides one of the most accurate representations of the different ice phases. To evaluate J, we apply the Lattice Mold technique, a computational method based on the use of molds to induce the nucleus formation from the metastable liquid under conditions at which observing spontaneous nucleation would be unfeasible. With this method, we obtain estimates of the nucleation rate for ice Ih and Ic and a stacking mixture of ice Ih/Ic, reaching consensus with most of the previously reported rates, although differing with some others. Furthermore, we confirm that the predicted nucleation rates obtained by the TIP4P/ICE model are in better agreement with experimental data than those obtained through the mW potential. Taken together, our study provides a reliable methodology to measure nucleation rates in a simple and computationally efficient manner that contributes to benchmarking the freezing behavior of two popular water models. [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 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]

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]

Modulating the interaction potential between colloids suspended in a fluid can trigger equilibrium phase transitions as well as the formation of non-equilibrium "arrested states," such as gels and glasses. Faithful representation of such interactions is essential for using simulation to interrogate the microscopic details of non-equilibrium behavior and for extrapolating observations to new regions of phase space that are difficult to explore in experiments. Although the extended law of corresponding states predicts equilibrium phases for systems with short-ranged interactions, it proves inadequate for equilibrium predictions of systems with longer-ranged interactions and for predicting non-equilibrium phenomena in systems with either short- or long-ranged interactions. These shortcomings highlight the need for new approaches to represent and disambiguate interaction potentials that replicate both equilibrium and non-equilibrium phase behavior. In this work, we use experiments and simulations to study a system with long-ranged thermoresponsive colloidal interactions and explore whether a resolution to this challenge can be found in regions of the phase diagram where temporal effects influence material state. We demonstrate that the conditions for non-equilibrium arrest by colloidal gelation are sensitive to both the shape of the interaction potential and the thermal quench rate. We exploit this sensitivity to propose a kinetics-based algorithm to extract distinct arrest conditions for candidate potentials that accurately selects between potentials that differ in shape but share the same predicted equilibrium structure. The algorithm selects the candidate that best matches the non-equilibrium behavior between simulation and experiments. Because non-equilibrium behavior in simulation is encoded entirely by the interparticle potential, the results are agnostic to the particular mechanism(s) by which arrest occurs, and so we expect our method to apply to a range of arrested states, including gels and glasses. Beyond its utility in constructing models, the method reveals that each potential has a quantitatively distinct arrest line, providing insight into how the shape of longer-ranged potentials influences the conditions for colloidal gelation. [ABSTRACT FROM AUTHOR]

Despite decades of intense research, whether the transformation of supercooled liquids into glass is a kinetic phenomenon or a thermodynamic phase transition remains unknown. Here, we analyzed optical microscopy experiments on 2D binary colloidal glass-forming liquids and investigated the structural links of a prominent kinetic theory of glass transition. We examined a possible structural origin for localized excitations, which are building blocks of the dynamical facilitation theory—a purely kinetic approach for the glass transition. To accomplish this, we utilize machine learning methods to identify a structural order parameter termed "softness" that has been found to be correlated with reorganization events in supercooled liquids. Both excitations and softness qualitatively capture the dynamical slowdown on approaching the glass transition and motivated us to explore spatial and temporal correlations between them. Our results show that excitations predominantly occur in regions with high softness and the appearance of these high softness regions precedes excitations, thus suggesting a causal connection between them. Thus, unifying dynamical and thermodynamical theories into a single structure-based framework may provide a route to understand the glass transition. [ABSTRACT FROM AUTHOR]

A deep understanding for collective behavior in an active matter system with complex interactions has far-reaching impact in biology. In the present work, we adopt Langevin dynamics simulations to investigate diffusion dynamics and phase separation in an anisotropic active particle system with a tunable biased angle α defined as the deviation between the active force direction and anisotropic orientation. Our results demonstrate that the biased angle can induce super-rotational diffusion dynamics characterized by a power-law relationship between the mean square angle displacement (MSAD) and the time interval Δt in the form of MSAD ∼ Δtβ with β > 1 and also result in non-trivial phase separation kinetics. As activity is dominant, nucleation time shows a non-monotonic dependence on the biased angle. Moreover, there arises a distinct transition of phase separation, from spinodal decomposition without apparent nucleation time to binodal decomposition with prominent nucleation delay. A significant inhibition effect occurs at right and obtuse angles, where the remarkable super-rotational diffusion prevents particle aggregation, leading to a slow nucleation process. As active force is competitive to anisotropic interactions, the system is almost homogeneous, while, intriguingly, we observe a re-entrant phase separation as a small acute angle is introduced. The prominent super-rotational diffusion under small angles provides an optimum condition for particle adsorption and cluster growth and, thus, accounts for the re-entrance of phase separation. A consistent scenario for the physical mechanism of our observations is achieved by properly considering the modulation of the biased angle on the interplay between activity and anisotropic interactions. [ABSTRACT FROM AUTHOR]

Room temperature adsorption of molecular iodine on Ag(111) has been studied by scanning tunneling microscopy (STM), low energy electron diffraction, Auger electron spectroscopy with factor analysis, and density functional theory (DFT). At the chemisorption stage, iodine first forms a ( 3 × 3 ) R30° structure. Further iodine dosing leads to continuous commensurate–incommensurate phase transition, taking place via the formation of striped superheavy domain walls. As a result, the uniaxially compressed (13 × 3 -R30°) phase is formed at an iodine coverage (θ) of 0.38 ML. At θ > 0.38 ML, first-order phase transition begins, leading to the formation of hexagonal moiré-like phases, which exhibit an anomalously large corrugation in STM (0.8–2.3 Å). In the range of 0.40–0.43 ML, the compression of hexagonal phases occurs, which ends at the formation of the (7 × 7)R21.8° structure at saturation. The DFT calculations allow us to explain the anomalous atomic corrugation of the hexagonal phases by the strong violation of the atomic structure of the substrate including up to ten layers of silver. Iodine dosing above 0.43 ML leads to the growth of 2D islands of silver iodide. The STM images of the silver iodide surface demonstrate a clear visible hexagonal superstructure with a periodicity of 25 Å superimposed with a quasi-hexagonal atomic modulation. DFT calculations of the atomic structure of AgI islands point to the formation of a sandwich-like double layer honeycomb structure similar to the case of I/Ag(100). [ABSTRACT FROM AUTHOR]

We present here the results of high-pressure broadband dielectric spectroscopy (BDS) measurements for a mixture of liquid-crystalline drug itraconazole (ITZ) and glycerol (GLY) at a critical concentration of 5% w/w in which the nematic order is eliminated. In the investigated system, smectic-A to isotropic phase transition leaves a clear fingerprint on the dielectric response, allowing for a phase diagram creation using BDS data. By following the α-relaxation dynamics under different thermodynamic conditions, we provide insights into the effect of pressure on temperature and the phenomenology of smectic-A to the isotropic phase transition. Additional measurements of specific volume as a function of pressure and temperature provide us with a deeper insight into material properties that could be analyzed comprehensively via the equation of state. We proved the validity of the density scaling concept, showing that the mixture's complexity does not exclude thermodynamic scaling of dynamic properties related to the α-process in the smectic-A phase. The low value of scaling exponent γ = 2.00 ± 0.02 and a high ratio of the activation energy at constant volume, EV, to the activation enthalpy at constant pressure, HP, indicate that temperature is a dominant variable controlling α-relaxation dynamics in the ordered smectic-A phase of the ITZ-GLY mixture. [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]

The phase behavior of a Biroli–Mézard model on the two dimensional square lattice in which hard-core particles can have at most three nearest neighboring occupied sites is investigated by means of grand-canonical Monte Carlo simulations. Finite-size scaling analysis of relevant thermodynamic quantities obtained via the histogram reweighting technique reveals that at high-density, the model undergoes a first-order phase transition with preferential sublattice occupation to a crystal phase with enantiomorph ground state configurations, in close analogy to the hard-core lattice gas with the exclusion range extended up to the third shell of nearest neighbors. [ABSTRACT FROM AUTHOR]

Monolayer and two-dimensional (2D) systems exhibit rich phase behavior, compared with 3D systems, in particular, due to the hexatic phase playing a central role in melting scenarios. The attraction range is known to affect critical gas–liquid behavior (liquid–liquid in protein and colloidal systems), but the effect of attraction on melting in 2D systems remains unstudied systematically. Here, we have revealed how the attraction range affects the phase diagrams and melting scenarios in a 2D system. Using molecular dynamics simulations, we have considered the generalized Lennard-Jones system with a fixed repulsion branch and different power indices of attraction from long-range dipolar to short-range sticky-sphere-like. A drop in the attraction range has been found to reduce the temperature of the gas–liquid critical point, bringing it closer to the gas–liquid–solid triple point. At high temperatures, attraction does not affect the melting scenario that proceeds through the cascade of solid–hexatic (Berezinskii–Kosterlitz–Thouless) and hexatic–liquid (first-order) phase transitions. In the case of dipolar attraction, we have observed two triple points inherent in a 2D system: hexatic–liquid–gas and crystal–hexatic–gas, the temperature of the crystal–hexatic–gas triple point is below the hexatic–liquid–gas triple point. This observation may have far-reaching consequences for future studies, since phase diagrams determine possible routes of self-assembly in molecular, protein, and colloidal systems, whereas the attraction range can be adjusted with complex solvents and external electric or magnetic fields. The results obtained may be widely used in condensed matter, chemical physics, materials science, and soft matter. [ABSTRACT FROM AUTHOR]

If a binary liquid mixture, composed of two alternative species with equal amounts, is quenched from a high temperature to a low temperature, below the critical point of demixing, then the mixture will phase separate through a process known as spinodal decomposition. However, if the two alternative species are allowed to interconvert, either naturally (e.g., the equilibrium interconversion of enantiomers) or forcefully (e.g., via an external source of energy or matter), then the process of phase separation may drastically change. In this case, depending on the nature of interconversion, two phenomena could be observed: either phase amplification, the growth of one phase at the expense of another stable phase, or microphase separation, the formation of nongrowing (steady-state) microphase domains. In this work, we phenomenologically generalize the Cahn–Hilliard theory of spinodal decomposition to include the molecular interconversion of species and describe the physical properties of systems undergoing either phase amplification or microphase separation. We apply the developed phenomenology to accurately describe the simulation results of three atomistic models that demonstrate phase amplification and/or microphase separation. We also discuss the application of our approach to phase transitions in polyamorphic liquids. Finally, we describe the effects of fluctuations of the order parameter in the critical region on phase amplification and microphase separation. [ABSTRACT FROM AUTHOR]

Titanium dioxide has been extensively studied in the rutile or anatase phase, while its high-pressure phases are less well-understood, despite that many are thought to have interesting optical, mechanical, and electrochemical properties. First-principles methods, such as density functional theory (DFT), are often used to compute the enthalpies of TiO2 phases at 0 K, but they are expensive and, thus, impractical for long time scale and large system-size simulations at finite temperatures. On the other hand, cheap empirical potentials fail to capture the relative stabilities of various polymorphs. To model the thermodynamic behaviors of ambient and high-pressure phases of TiO2, we design an empirical model as a baseline and then train a machine learning potential based on the difference between the DFT data and the empirical model. This so-called Δ-learning potential contains long-range electrostatic interactions and predicts the 0 K enthalpies of stable TiO2 phases that are in good agreement with DFT. We construct a pressure–temperature phase diagram of TiO2 in the range 0 < P < 70 GPa and 100 < T < 1500 K. We then simulate dynamic phase transition processes by compressing anatase at different temperatures. At 300 K, we predominantly observe an anatase-to-baddeleyite transformation at about 20 GPa via a martensitic two-step mechanism with a highly ordered and collective atomic motion. At 2000 K, anatase can transform into cotunnite around 45–55 GPa in a thermally activated and probabilistic manner, accompanied by diffusive movement of oxygen atoms. The pressures computed for these transitions show good agreement with experiments. Our results shed light on how to synthesize and stabilize high-pressure TiO2 phases, and our method is generally applicable to other functional materials with multiple polymorphs. [ABSTRACT FROM AUTHOR]

Raman spectroscopy investigations on L-methionine (L-Met) performed in a large temperature range (170–420 K) and in a wide spectral window (5–3600 cm−1) have revealed an extended disordering mechanism triggered by thermally activated motions of the terminal side-chain atoms, from 250 up to 390 K. This very progressive disordering process is characterized by two thermodynamic features, the first corresponding to a broad endotherm (250 → 310 K) marking the beginning of the process, while the second ending the disordering transformation is a sharp endothermic peak at 390 K. These thermodynamic events are correlated with the softening of lattice vibrations and the increase of the quasielastic scattering, considered as the signatures of displacive phase transitions. The amorphous-like band-shape of the low-frequency Raman spectrum collected above 390 K, resulting from the strong anharmonicity of local motions, is contrasting with the detection of additional Bragg peaks above 390 K by x-ray diffraction, consistent with the Cp jump accompanying the endothermic peak. These observations suggest that L-Met is progressively dynamically disordered adopting additional configurations in the crystalline lattice through rotations of CH3 and the side-chain flexibility not clearly detected by x-ray diffraction. These results should be crucial for considering the stability of dried proteins composed of methionine residues. [ABSTRACT FROM AUTHOR]