Oxygen and other light elements comprise up to 5 wt% of the Earth's outer-core, and may significantly influence its physical properties and the operation of the geodynamo. Here we report in situ x-ray diffraction measurements of Fe, Fe + 4.5 FeO (atomic proportion), and Fe2O3 melts at 177-438 GPa, achieved using laser-driven shock compression at an x-ray free-electron laser. The melts exhibit Fe-O coordination numbers between 4.0(0.4) and 4.5(0.4), indicating predominantly four-fold coordination environments. These coordination states are significantly smaller than those of Fe-bearing lower-mantle phases such as bridgmanite and ferropericlase. Shorter Fe-Fe interatomic distances in compressed iron oxide melts drive the denser packing relative to ambient melts, while the structural differences between Fe + 4.5 FeO and Fe2O3 melts under shock indicate that the oxidation state modulates oxygen solubility in liquid Fe. At around 177 GPa (380 km below the core-mantle boundary), Fe2O3 melts exhibit higher Fe-O coordination, suggesting that local variations in oxygen content could contribute to the stratification in the uppermost outer-core inferred from seismological and geomagnetic obse
Block copolymer (BCP) melts play a critical role in the design of thermoplastics, owing in large part to the creation of alternating nano-scale domains of soft and stiff components. Considerable attention has been given to the short-to-intermediate time response of BCP melts, when the storage modulus is expected to dominate the viscoelastic properties. However, less attention has been paid to the long-time relaxation and rigidity of microphase separated BCP melts or the role that domain morphology plays in modulating near-equilibrium response. We take advantage of the ability of self-consistent field theory (SCFT) to calculate equilibrium properties of BCP melts to explore the anisotropic elastic response of ordered ABA and AB copolymer melts as quasistatic deformation processes. This allows us to determine the anisotropic stiffness of the liquid crystal-like lamellar and columnar phases due to modulations in domain spacing, as well as the full stiffness tensor of the cubic BCC sphere and double gyroid phases. We explore elastic modulus landscapes for both single grain materials and random polycrystals over architectural parameters and segregation strengths, using AB diblock and AB
We estimate the plateau moduli of highly entangled end-pinned bead-spring polymer melts with Z = 100 and Z = 200 from the time-dependent elastic response to a step strain, which we first extrapolate to infinite time and then interpolate to zero strain. We present data for systems deformed in the melt state as well as for systems deformed at the primitive path level following the recent iPPA protocol. We observe excellent agreement between the plateau moduli obtained via the two deformation protocols and good agreement with the available experimental data for commodity polymer melts using a common mapping on the Kuhn scale.
In liquid pnictogens, quasi-stable structures can be formed near melting temperature. The nature of their stability does not have the unified point of view. In the present work, the task of determining the degree of atomic bonding in these structures is solved using the Crystal Orbital Hamilton Population (COHP) method. The original results of ab-initio simulation of arsenic, antimony and bismuth melts near their melting temperatures are used. It is shown that the features of the electron interaction at the level of $p$-orbitals determine the characteristic bond lengths and angles between atoms. It has been established that the stability of structures decreases according to a power law with an increase in the atomic mass of a chemical element and the number of atoms in the structure. The obtained results clarify the understanding the mechanisms of formation of quasi-stable structures in pnictogen melts from first principles.
We investigate with numerical simulations the influence of topology and stiffness on macroscopic rheological properties of polymer melts consisting of unknotted, knotted or concatenated rings. While melts of flexible, knotted oligomer rings tend to be significantly more viscous than their unknotted counterparts, differences vanish in a low shear rate scenario with increasing degree of polymerization. Melts of catenanes consisting of two rings on the other hand are consistently more viscous than their unconcatenated counterparts. These topology-based differences in rheological properties can be exploited to segregate mixtures of otherwise chemically similar polymers, e.g., in microfluidic devices, which is demonstrated by exposing a blend of flexible knotted and unknotted oligomer rings to channel flow.
The effect of shear flows on the thermal conductivity of polymer melts is investigated using a reversed nonequilibrium molecular dynamics (RNEMD) method. We extended the original RNEMD method to simultaneously produce spatial gradients of temperature and flow velocity in a single direction. This method enables accurate measurement of the thermal conductivity in the direction transverse to shear flow. The Weissenberg number defined with the shear rate and the relaxation time of the polymer conformation can uniformly differentiate the occurrence of shear rate dependence of the thermal conductivity across different chain lengths. The stress-thermal rule (STR) (i.e., the linear relationship between anisotropic parts of the stress tensor and the thermal conductivity tensor) holds for entangled polymer melts even under shear flows but not for unentangled polymer melts. Furthermore, once entanglements form in polymer chains, the stress-thermal coefficient in the STR remains independent of the polymer chain length. These observations align with the theoretical foundation of the STR, which focuses on energy transmission along the network structure of entangled polymer chains [Brule, Rheol.
The viscous flow of polymer chains in dense melts is dominated by topological constraints whenever the single chain contour length, N, becomes larger than the characteristic scale Ne, defining comprehensively the macroscopic rheological properties of the highly entangled polymer systems. Even though the latter are naturally connected to the presence of hard constraints like knots and links within the polymer chains, the difficulty of integrating the rigorous language of mathematical topology with the physics of polymer melts has limited somehow a genuine topological approach to the problem of classifying these constraints and to how they are related to the rheological entanglements. In this work, we tackle this problem by studying the occurrence of knots and links in lattice melts of randomly knotted and randomly concatenated ring polymers of various bending stiffness. Specifically, by introducing an algorithm which shrinks the chains to their minimal shapes which do not violate topological constraints and by analyzing those in terms of suitable topological invariants, we provide a detailed characterization of the topological properties at the intra-chain level (knots) and of links
Polymer melts undergoing large deformation by uniaxial elongation are studied by molecular dynamics simulations of bead-spring chains in melts. Applying a primitive path analysis to strongly deformed polymer melts, the role of topological constrains in highly entangled polymer melts is investigated and quantified. We show that the over-all, large scale conformations of the primitive paths (PPs) of stretched chains follow affine deformation while the number and the distribution of entanglement points along the PPs do not. Right after deformation, PPs of chains retract in both directions parallel and perpendicular to the elongation. Upon further relaxation we observe a long-lived clustering of entanglement points. Together with the delayed relaxation time this leads to a metastable inhomogeneous distribution of topological constraints in the melts.
In polymer physics it is typically assumed that excluded volume interactions are effectively screened in polymer melts. Hence, chains could be described by an effective random walk without excluded volume interactions. In this letter, we show that this mapping is problematic by analyzing the occurrence of knots, their spectrum and sizes in polymer melts, corresponding random walks and chains in dilute solution. The effective random walk severely overrates the occurrence of knots and their complexity, particularly when compared to melts of flexible chains, indicating that non-trivial effects due to remnants of self-avoidance still play a significant role for the chain lengths considered in this numerical study. For melts of semiflexible chains, the effect is less pronounced. In addition, we find that chains in a melt are very similar in structure and topology to dilute single chains close to the collapse transition, which indicates that the latter are also not well-represented by random walks. We finally show that typical equilibration procedures are well-suited to relax the topology in melts.
Equilibration of polymer melts containing highly entangled long polymer chains in confinement or with free surfaces is a challenge for computer simulations. We approach this problem by first studying polymer melts based on the soft-sphere coarse-grained model confined between two walls with periodic boundary conditions in two directions parallel to the walls. Then we apply backmapping to reinsert the microscopic details of the underlying bead-spring model. Tuning the strength of the wall potential, the monomer density of confined polymer melts in equilibrium is kept at the bulk density even near the walls. In a weak confining regime, we observe the same conformational properties of chains as in the bulk melt showing that our confined polymer melts have reached their equilibrated state. Our methodology provides an efficient way of equilibrating large polymer films with different thicknesses and is not confined to a specific underlying microscopic model. Switching off the wall potential in the direction perpendicular to the walls, enables to study free-standing highly entangled polymer films or polymer films with one supporting substrate.
The widely used double-bridging hybrid (DBH) method for equilibrating simulated entangled polymer melts [R. Auhl et al., J. Chem. Phys. v. 119, p. 12718, 2003] loses its effectiveness as chain stiffness increases into the semiflexible regime because the energy barriers associated with double-bridging Monte Carlo moves become prohibitively high. Here we overcome this issue by combining DBH with the use of core-softened pair potentials. This reduces the energy barriers substantially, allowing us to equilibrate melts with $N \simeq 40 N_e$ and chain stiffnesses all the way up to the isotropic-nematic transition using simulations of no more than 100 million timesteps. For semiflexible chains, our method is several times faster than standard DBH; we exploit this speedup to develop improved expressions for Kremer-Grest melts' chain-stiffness-dependent Kuhn length $\ell_K$ and entanglement length $N_e$.
The lattice cluster theory (LCT) for semiflexible linear telechelic melts, developed in paper I, is applied to examine the influence of chain stiffness on the average degree of self-assembly and the basic thermodynamic properties of linear telechelic polymer melts. Our calculations imply that chain stiffness promotes self-assembly of linear telechelic polymer melts that assemble on cooling when either polymer volume fraction $φ$ or temperature $T$ is high, but opposes self-assembly when both $φ$ and $T$ are sufficiently low. This allows us to identify a boundary line in the $φ$-$T$ plane that separates two regions of qualitatively different influence of chain stiffness on self-assembly. The enthalpy and entropy of self-assembly are usually treated as adjustable parameters in classical Flory-Huggins type theories for the equilibrium self-assembly of polymers, but they are demonstrated here to strongly depend on chain stiffness. Moreover, illustrative calculations for the dependence of the entropy density of linear telechelic polymer melts on chain stiffness demonstrate the importance of including semiflexibility within the LCT when exploring the nature of glass formation in models o
The paper presents the results of molecular dynamics study of the viscosity of nickel-containing binary metal melts for a wide range of temperatures, including the region of the equilibrium liquid phase and supercooled melt. It is shown that the temperature dependencies of the viscosity of binary metal melts are described by the Kelton's quasi-universal model. Based on the analysis of the viscosity coefficient of the binary melt composition within the framework of the Rosenfeld's scale transfor\-mations, it has been established that to correctly describe the viscosity of binary/multi\-component metal melts within the framework of entropy models, it is necessary to use a more complex representation of the excess entropy $S_{ex}$ than in the approximation of pair correlation entropy $S_2$.
The generalized entropy theory (GET) offers many insights into how molecular parameters influence polymer glass-formation. Given the fact that chain rigidity often plays a critical role in understanding the glass-formation of polymer materials, the GET was originally developed based on models of semiflexible chains. Consequently, all previous calculations within the GET considered polymers with some degree of chain rigidity. Motivated by unexpected results from computer simulations of fully flexible polymer melts concerning the dependence of thermodynamic and dynamic properties on the cohesive interaction strength ($ε$), the present paper employs the GET to explore the influence of $ε$ on glass-formation in models of polymer melts with a vanishing bending rigidity, i.e., fully flexible polymer melts. In accord with simulations, the GET for fully flexible polymer melts predicts that basic dimensionless thermodynamic properties (such as the thermal expansion coefficient and isothermal compressibility) are universal functions of the temperature scaled by $ε$ in the regime of low pressures. Similar scaling behavior is also found for the configurational entropy density in the GET for fu
Chromite is a key magmatic mineral frequently used as petrogenetic indicator of physico-chemical conditions of mafic magma crystallization. In this work, magnesiochromite and chromite solubility in a natural basalt and an iron-free haplobasalt was investigated at 1440°C and atmospheric pressure under controlled CO-CO2 gas mixtures corresponding to the range two log units below to two log units above the fayalite-magnetite-quartz buffer. The source of chromium was either natural chromite or synthetic Cr2O3, the latter reacting with the basaltic liquids to form a magnesiochromite. The highest concentrations of Cr in haplobasaltic melts are found in equilibrium with magnesiochromite, depending on redox conditions. In detail, at low fO2, liquids have high Cr contents, but the variation of log [Cr, ppm] is not a linear function of log fO2. Using our new data and data from the literature a model for Cr concentrations at chromite/magnesiochromite saturation in silicate melts has been developed based upon a thermodynamic formalism. Our model may be used to assess the effect of melt composition on chromite and magnesiochromite solubility in silicate melts during peridotite melting and assim
We present an effective and simple multiscale method for equilibrating Kremer Grest model polymer melts of varying stiffness. In our approach, we progressively equilibrate the melt structure above the tube scale, inside the tube and finally at the monomeric scale. We make use of models designed to be computationally effective at each scale. Density fluctuations in the melt structure above the tube scale are minimized through a Monte Carlo simulated annealing of a lattice polymer model. Subsequently the melt structure below the tube scale is equilibrated via the Rouse dynamics of a force-capped Kremer-Grest model that allows chains to partially interpenetrate. Finally the Kremer-Grest force field is introduced to freeze the topological state and enforce correct monomer packing. We generate $15$ melts of $500$ chains of $10.000$ beads for varying chain stiffness as well as a number of melts with $1.000$ chains of $15.000$ monomers. To validate the equilibration process we study the time evolution of bulk, collective and single-chain observables at the monomeric, mesoscopic and macroscopic length scales. Extension of the present method to longer, branched or polydisperse chains and/or
The electrical conductivity of melts of boron and its carbide (B4C), nitride (BN), and phosphide (BP) has been studied at pressures to 7.7 GPa and temperatures to 3500 K. It has been shown that these melts are good conductors with specific electrical conductivity values comparable with that of iron melt at ambient pressure.
Using in-situ synchrotron tomography, we investigate the coarsening dynamics of barium borosilicate melts during phase separation. The 3-D geometry of the two interconnected phases is determined thanks to image processing. We observe a linear growth of the size of domains with time, at odds with the sublinear diffusive growth usually observed in phase-separating glasses or alloys. Such linear coarsening is attributed to viscous flow inside the bicontinuous phases, and quantitative measurements show that the growth rate is well explained by the ratio of surface tension over viscosity. The geometry of the domains is shown to be statistically similar at different times, provided that the microstructure is rescaled by the average domain size. Complementary experiments on melts with a droplet morphology demonstrate that viscous flow prevails over diffusion in the large range of domain sizes measured in our experiments (1 - 80 microns).
Despite the long-recognized fact that chemical structure and specific interactions greatly influence the thermodynamic properties of polymer systems, a predictive molecular theory that enables systematically addressing the role of chemical structure and specific interactions has been slow to develop even for polymer melts. While the lattice cluster theory (LCT) provides a powerful vehicle for understanding the influence of various molecular factors, such as monomer structure, on the thermodynamic properties of polymer melts and blends, the application of the LCT has heretofore been limited to the use of the simplest polymer model in which all united atom groups within the monomers of a species interact with a common monomer averaged van der Waals energy. Thus, the description of a compressible polymer melt involves a single van der Waals energy. As a first step towards developing more realistic descriptions to aid in the analysis of experimental data and the design of new materials, the LCT is extended here to treat models of polymer melts in which the backbone and side groups have different interaction strengths, so three energy parameters are present, namely, backbone-backbone, s
Melting is omnipresent in nature and technology, with applications ranging from metallurgy, biology, food science, and latent thermal energy storage to oceanography, geophysics, and climate science, and occurring on all scales from sub-millimeter to global scales. The key objective is to understand the rate at which an object melts as a function of its size and of the ambient conditions. To achieve this it is important to be able to extrapolate from small scale experiments and observations to large or even global scales. This is done by scaling laws. However, these are only meaningful if there is no transition from one scaling relation to another one. Here we show, however, that for both fixed and freely-advected melting objects immersed in a turbulent flow a melting transition does exist, namely from slow melting at the small scales to fast melting at the large scales. We do so by controlled melting experiments and corresponding direct numerical simulations, covering four orders of magnitude in scale. The transition corresponds to the transition from a laminar-type boundary layer around the melting object to a turbulent-type boundary layer, i.e., from so-called classical turbulenc