Fluid State Research Articles

Hole statistics of equilibrium 2D and 3D hard-sphere crystals.

The probability of finding a spherical "hole" of a given radius r contains crucial structural information about many-body systems. Such hole statistics, including the void conditional nearest-neighbor probability functions GV(r), have been well studied for hard-sphere fluids in d-dimensional Euclidean space Rd. However, little is known about these functions for hard-sphere crystals for values of r beyond the hard-sphere diameter, as large holes are extremely rare in crystal phases. To overcome these computational challenges, we introduce a biased-sampling scheme that accurately determines hole statistics for equilibrium hard spheres on ranges of r that far extend those that could be previously explored. We discover that GV(r) in crystal and hexatic states exhibits oscillations whose amplitudes increase rapidly with the packing fraction, which stands in contrast to GV(r) that monotonically increases with r for fluid states. The oscillations in GV(r) for 2D crystals are strongly correlated with the local orientational order metric in the vicinity of the holes, and variations in GV(r) for 3D states indicate a transition between tetrahedral and octahedral holes, demonstrating the power of GV(r) as a probe of local coordination geometry. To further study the statistics of interparticle spacing in hard-sphere systems, we compute the local packing fraction distribution f(ϕl) of Delaunay cells and find that, for d ≤ 3, the excess kurtosis of f(ϕl) switches signat a certain transitional global packing fraction. Our accurate methods to access hole statistics in hard-sphere crystals at the challenging intermediate length scales reported here can be applied to understand the important problem of solvation and hydrophobicity in water at such length scales.

Holographic fluids from 5D dilaton gravity

We study a solvable class of five-dimensional dilaton gravity models that continuously interpolate between anti-de Sitter (AdS5), linear dilaton (LD5) and positively curved spacetimes as a function of a continuous parameter ν. The dilaton vacuum expectation value is set by a potential localized on a flat brane. We chart the elementary properties of these backgrounds for any admissible ν, and determine stability conditions of the brane-dilaton system. We find that the spectrum of metric fluctuations can be either continuous or discrete. It features a massless graviton mode confined between the brane and the curvature singularity, and a massive radion mode tied to brane-dilaton stability. We show that, in the presence of a bulk black hole, the holographic theory living on the brane features a perfect fluid. The equation of state of the holographic fluid interpolates between radiation, pressureless matter and vacuum energy as a function of ν. This extends earlier findings on holographic fluids. Our results suggest that the thermodynamics of the fluid mirrors precisely the thermodynamics of the bulk black hole.

All-optical method to directly measure the pressure–volume–temperature equation of state of fluids in the diamond anvil cell

We have developed a new all-optical method to directly measure the pressure–volume–temperature (PVT) equation of state (EOS) of fluids and transparent solids in the diamond anvil high pressure cell by measuring the volume of the sample chamber. Our method combines confocal microscopy and white light interference with a new analysis method, which exploits the mutual dependence of sample density and refractive index: Experimentally, the refractive index determines the measured sample chamber thickness (and therefore the measured sample volume/density), yet the sample density is by far the dominant factor in determining the variation in the refractive index with pressure. Our analysis method allows us to obtain a set of values for the density and refractive index, which are mutually consistent and agree with the experimental data within error. We have conducted proof-of-concept experiments on a variety of samples (H2O, CH4, C2H6, C3H8, KCl, and NaCl) at ambient temperature and at high temperatures up to just above 500 K. Our proof-of-concept data demonstrate that our method is able to reproduce known fluid and solid EOS within error. Furthermore, we demonstrate that our method allows us to directly and routinely measure the PVT EOS of simple fluids at GPa pressures up to, at least, 514 K (the highest temperature reached in our study). A reasonable estimation of the known sources of error in our volume determinations indicates that the error is currently ±2.7% at high temperature and that it is feasible to reduce it to ca. ±1% in future work.

Deciphering ore-forming fluid evolution processes of the Yuerya gold deposit, eastern Hebei, China: Insights from pyrite texture, trace element and sulfur isotope compositions

The Yuerya gold deposit, a representative gold deposit associated with Mesozoic granites in the eastern Hebei region, is located in the eastern Hebei-western Liaoning gold ore-concentrated area, the northern margin of the North China Plate. The orebody is mainly preserved in the Yanshanian granite and its contact zone with the carbonate rocks of the Gaoyuzhuang Formation, and the main ore mineral is pyrite. Based on the pyrite texture, geochemical characteristics of the in-situ trace element analysis, and sulfur isotope, the pyrite in the mining area is classified into five generations (Py1, Py2, Py3, Py4a, and Py4b). Five generations of pyrite all have an average Te/Se ratio > 1, and their average Co/Ni ratios are 1087.27, 43.97, 4.17, 12.62, and 4.99, respectively, indicating that the ore-forming fluid predominantly originates from magmatic-hydrothermal. In-situ sulfur isotopic analysis using LA–MC–ICP–MS shows that the δ34S values of the five generations of pyrites ranged from 0.3 ‰ to 9.6 ‰; the mean values are 3.8 ‰ (Py1, n = 4), 4.2 ‰ (Py2, n = 10), 2.6 ‰ (Py3, n = 18), 4.0 ‰ (Py4a, n = 20) and 4.6 ‰ (Py4b, n = 29), respectively. Approximately 32 % of the results exceed 5.0 ‰, surpassing the values indicative of a magmatic-hydrothermal source. Thermodynamic simulations reveal that the fluctuations in the oxidation state of the ore-forming fluid is not the predominant cause of δ34S enrichment in the pyrite. Our studies suggest that sulfur is primarily derived from mantle-crust mixtures, with some contributions from the wallrock of the Gaoyuzhuang Formation. The results from LA–ICP–MS trace element analysis of pyrites reveal a significant enrichment of As-Bi-Te-Tl alongside Au in the ore-forming fluid. Additionally, the results of ore-forming elements in three typical cross sections of different wallrocks in the mining area also show that Au, Ag, As, and Pb are closely related. Therefore, it is considered that the over-enrichment of As could be an important factor leading to the precipitation of Au. Pyrite undergoes expansion of the mineral lattice parameters or lattice dislocation due to the substitution of S with As, creating space for the growth of solid-solution Au. Subsequently, some of the solid-solution Au within pyrite is liberated through internal oscillations, forming the visible gold particles.

Photocontrolled Reversible Solid-Fluid Transitions of Azopolymer Nanocomposites for Intelligent Nanomaterials.

Intelligent polymer nanocomposites are multicomponent and multifunctional materials that show immense potential across diverse applications. However, to exhibit intelligent traits such as adaptability, reconfigurability and dynamic properties, these materials often require a solvent or heating environment to facilitate the mobility of polymer chains and nanoparticles, rendering their applications in everyday settings impractical. Here intelligent azopolymer nanocomposites that function effectively in a solvent-free, room-temperature environment based on photocontrolled reversible solid-fluid transitions via switching flow temperatures (Tfs) are shown. A range of nanocomposites is synthesized through the grafting of Au nanoparticles, Au nanorods, quantum dots, or superparamagnetic nanoparticles with photoresponsive azopolymers. Leveraging the reversible cis-trans photoisomerization of azo groups, the azopolymer nanocomposites transition between solid (Tf above room temperature) and fluid (Tf below room temperature) states. Such photocontrolled reversible solid-fluid transitions empower the rewriting of nanopatterns, correction of nanoscale defects, reconfiguration of complex multiscale structures, and design of intelligent optical devices. These findings highlight Tf-switchable polymer nanocomposites as promising candidates for the development of intelligent nanomaterials operative in solvent-free, room-temperature conditions.

An improved artificial spin ice structure for restoring ice degeneracy

It is known that rare-earth-based pyrochlore oxide may accommodate the well-defined two-in-two-out spin ice state with its tetrahedral unit. Low-energy excitation is argued to favor the highly concerned monopole state which attracts essential attention. However, such an excitation cannot be directly tracked and imaged using advanced characterizations, raising challenges to our understanding of the physics of monopoles. In this work, we propose an improved two-dimensional artificial spin ice structure on the Shastry–Sutherland lattice to restore the degeneracy of realistic pyrochlore systems. Such a structure avoids the deficiency of inequivalent nearest and next-nearest exchanges in the planar quadrate unit, which, however, is equivalent to the tetrahedral unit of realistic pyrochlore oxides. Therefore, this spin ice model restores state degeneracy that is lost in conventional planar artificial spin ice structures, representing an improved simulator of real spin ice systems. Our careful investigations of such improved structures reveal the rich physics of spin ice excitations, including the phase diagram, which allows different ordered phases and interesting critical phase transitions between spin ice phase I and phase II. Energy spectrum analysis suggests that restoration of state degeneracy substantially reduces monopole excitation energy, resulting in a striking monopole emergency at the critical point. Furthermore, the emergent spin dimer phase in this improved model allows high-density monopole excitations and exhibits high-correlated monopole fluid states.

Surface Chemistry of Cytosporone-B Incorporated in Models for Microbial Biomembranes as Langmuir Monolayers.

Cytosporone-B, a polyketide renowned for its antimicrobial properties, was integrated into Langmuir monolayers composed of dipalmitoylphosphoethanolamine (DPPE) and dioleoylphosphoethanolamine (DOPE) lipids, effectively emulating microbial cytoplasmic membranes. This compound exhibited an expansive influence on DPPE monolayers while inducing condensation in DOPE monolayers. This led to a notable reduction in the compressibility modulus for both lipids, with a more pronounced effect observed for DPPE. The heightened destabilization observed in DOPE monolayers subjected to biologically relevant pressures was particularly noteworthy, as evidenced by surface pressure-time curves at constant area. In-depth analysis using infrared spectroscopy at the air-water interface unveiled alterations in the alkyl chains of the lipids induced by cytosporone-B. This was further corroborated by surface potential measurements, indicating a heightened tilt in the acyl chains upon drug incorporation. Notably, these observed effects did not indicate an aggregating process induced by the drug. Overall, the distinctive impact of cytosporone-B on each lipid underscores the importance of understanding the nuanced effects of microbial drugs on membranes, whether in condensed or fluid states.

Development of performance analysis model with numerical and dynamic energy modeling for an integrated PVT-HP system

This study delivers a comprehensive numerical analysis of an integrated photovoltaic/thermal-heat pump (PVT-HP) system, underscored by an innovative three-dimensional thermal energy storage (TES) model. Addressing the limitations of conventional one-dimensional TES models, which struggled with accurately simulating horizontal flows and the fluid state at various facility inlets, this research marks a significant advancement. Our development of a 3D TES model optimizes finite elements to reduce computation time and enhance accuracy, vital for precise energy analysis in buildings and air conditioning systems. The integration of this model with a network model allows for detailed analysis of flow rate changes and their impacts on heat transfer, crucial for the integrated PVT-HP system. Key findings include a 30 % increase in the temperature difference between the top and bottom layers of the TES when the flow rate is reduced from 50 LPM to 20 LPM, and a 1 % decrease in the heat pump's coefficient of performance (COP) with a 10 LPM reduction in flow rate. The study emphasizes that managing the flow rate is key to optimizing temperature gradients and system efficiency. Strategic adjustments in flow rate significantly affect thermal storage and transfer capabilities, proving essential in system performance. This research offers vital insights into optimizing PVT-HP systems, focusing on the importance of flow rate control in improving temperature management and overall system performance.

Ionic strength effect on the structure and dynamics of colloidal dispersions with weak attractive interactions

We present a macro- and microrheological study of the effect of ionic strength on the phase behavior, structure and flow properties of colloidal dispersions stabilized by short-range repulsive interactions inferred via co-polymerization of weak acid groups. We investigated the impact of ionic strength, by varying the concentration of natrium chloride (NaCl), on dispersions with only repulsive interactions, but also when additional attractive depletion forces are present due to added non-absorbent polymer, namely polyethyleneoxide (PEO). For dispersions with only repulsive electrosteric interactions, at low particle volume fractions (ϕ < 0.4), increasing ionic strength hardly affects the relative viscosity whereas at higher particle volume fractions, a decrease in viscosity is observed due to a reduced range of electrosteric repulsion between particles, corresponding to a reduction of the effective particle volume fraction ϕeff. For dispersions including attractive interactions, at volume fractions ϕ = 0.45 below the hard sphere freezing point (ϕc=0.5), independently of the ionic strength, the bulk viscosity increases monotonically with increasing PEO concentration due to the transition from a fluid to a fluid/crystalline and finally to a gel state. However, the increase in ionic strength shifts the concentrations of both phase transitions to lower polymer concentrations, indicating that in presence of salt a weaker attraction is required to induce these transitions. At ϕ=0.52, just above ϕc, for dispersions without added salt, broadening of the fluid-crystalline coexistence regime, due to added PEO, results in a viscosity minimum corresponding to different size and packing density of the crystalline regions as observed earlier in Weis et al. [1]. When salt is added, the initial state for dispersions without PEO is a fluid state due to the reduction of ϕeff and the addition of a small amount of PEO leads to an increase in bulk viscosity due to the formation of large crystals for dispersions containing 10 mM NaCl and a network of small, dense crystals of size ≈10 µm when 100 mM NaCl is added. At ϕ=0.54, a result similar to ϕ=0.52 is obtained for dispersions without added salt, whereas in presence of salt, we find that the fluid/crystalline coexistence regime narrows down as the range of electrosteric repulsion decreases and the addition of PEO results in gel formation and gels become more uniform with increasing attraction strength. These investigations demonstrate how superimposed short range electrosteric repulsion and weak depletion attraction affect microstructure and flow behavior of colloidal dispersions.

Breaking the limitation of terahertz resonances in ferrites through 4D printing of metamaterials

4D printing allows the creation of 3D-printed structures with a variety of new functionalities that respond to external stimuli, including electromagnetic properties, for which research is limited, such as terahertz resonances in ferrite materials; and therefore holding great promise for a wide range of applications. However, the limited physical properties of ferrite materials, including their anisotropic properties and limited material manipulation space, pose significant challenges to the development of ferrite materials designed for progressively increasing terahertz frequencies. Therefore, this paper proposes a platform for overcoming the limitations of terahertz resonance in ferrites through the 4D printing of metamaterials. Both the elastic and viscous moduli of the slurry are sensitive to the applied shear stress; as this stress increases, these moduli sharply decrease by several orders of magnitude, resulting in the transformation of the slurry into a fluid state similar to the liquefaction of ink. Electromagnetic coupling modes between electric dipoles and magnetic dipoles endow these ferrite metamaterials with induction and switching of terahertz resonances due to orientation stimulation. Our work sheds light on the creation of electromagnetic properties and their control by 4D printing.

Quantitative response of the spatial distribution of diesel oil to freezing and thawing temperatures in groundwater

The mobilization and redistribution of organic contaminants in groundwater is the basis and key to explore its dynamic evolution and appropriate remediation. The naturally occurring diametrical temperature gradient during freezing and thawing cycle leads to distinct behaviors of organic contaminants in groundwater. In this study, the pore-scale distribution of diesel oil in the porous media was quantitatively divided into capillary fluid state (CFS) and free fluid state (FFS) based on multiphase flow dynamics, employing low-field nuclear magnetic resonance (LF-NMR) technology. The pore-scale distribution of diesel oil depends not only on the freezing and thawing cycle but also on the temperature gradient according to LF-NMR results. The content of diesel oil in the CFS generally increases with a positive temperature gradient (e.g. freezing) compared to a negative temperature gradient (e.g. thawing), while the content of diesel oil in the FFS generally decreases. This dependence of the temperature gradient on pore-scale distribution of the diesel oil is positively correlated with the particle size of the porous medium. Furthermore, the pore-scale distribution of the diesel oil during the freezing and thawing cycle is influenced by the kinematic viscosity of the diesel oil. There is an exponential relationship between the diesel oil content and the kinematic viscosity, independent of the freezing or thawing process. During the freezing process, the diesel oil migrates from FFS to CFS, while this migration is reversed during the thawing process. The reverse migration of the diesel oil between the freezing and thawing processes leads to a spatial redistribution of the diesel oil, which is controlled by both the fluid energy and the capillary force. The present work provide meaningful guidance for the remediation of groundwater contamination in cold regions.

Molecular dynamics simulations of wettabilities on shale in three- and four-phase systems

The interfacial properties between three-phase fluid and shale are crucial in carbon dioxide flooding and sequestration. Here, we present an investigation on fluid-solid interfacial properties in the H2O/CO2/decane/shale four-phase system under geological conditions using molecular simulations. The four-phase system exhibits distinct three-phase contacts due to the varying occurrence state of fluids. The method (J. Colloid Interface Sci. 585 (2021) 320-327) for calculating the wettabilities of three-phase systems is not directly applicable to four-phase systems due to the lack of interfaces. The method is extended to four-phase systems by conducting simultaneous simulations of two such systems with different fluid phase positions, ensuring all interfaces are accessible. We found that the relative interfacial free energies at the fluid-solid interfaces and the wettabilities of the H2O/CO2/shale and CO2/decane/shale contacts in the four-phase system are similar to those observed in the corresponding three-phase systems. However, the wettabilities of the H2O/decane/shale contact in the four-phase system significantly differ from those of the H2O/decane/shale three-phase system. The introduction of a CO2 phase not only enhances the wettability but also amplifies this effect, especially under conditions of lower temperatures and higher pressures. For instance, the wettability increases by 35° at 298 K and 6 MPa. This variation in wettabilities is primarily attributed to the substantial influence of CO2 on interfacial energies in the interfaces of decane-solid and H2O-decane. The elevated impact of CO2 at lower temperatures and higher pressures is explained by the enhanced CO2 adsorption in both fluid-solid and fluid-fluid interfaces.

Nonequilibrium hyperuniform states in active turbulence

We demonstrate that the complex spatiotemporal structure in active fluids can feature characteristics of hyperuniformity. Using a hydrodynamic model, we show that the transition from hyperuniformity to nonhyperuniformity and antihyperuniformity depends on the strength of active forcing and can be related to features of active turbulence without and with scaling characteristics of inertial turbulence. Combined with identified signatures of Levy walks and nonuniversal diffusion in these systems, this allows for a biological interpretation and the speculation of nonequilibrium hyperuniform states in active fluids as optimal states with respect to robustness and strategies of evasion and foraging.

Arrested states in colloidal fluids with competing interactions: A static replica study.

We present the first systematic application of the integral equation implementation of the replica method to the study of arrested states in fluids with microscopic competing interactions (short-range attractive and long-range repulsive, SALR), as exemplified by the prototype Lennard-Jones-Yukawa model. Using a wide set of potential parameters, we provide as many as 11 different phase diagrams on the density (ρ)-temperature (T) plane, embodying both the cluster-phase boundary, TC(ρ), and the locus below which arrest takes place, TD(ρ). We describe how the interplay between TC and TD-with the former falling on top of the other, or the other way around, depending on thermodynamic conditions and potential parameters-gives rise to a rich variety of non-ergodic states interspersed with ergodic ones, of which both the building blocks are clusters or single particles. In a few cases, we find that the TD locus does not extend all over the density range subtended by the TC envelope; under these conditions, the λ-line is within reach of the cluster fluid, with the ensuing possibility to develop ordered microphases. Whenever a comparison is possible, our predictions favorably agree with previous numerical results. Thereby, we demonstrate the reliability and effectiveness of our scheme to provide a unified theoretical framework for the study of arrested states in SALR fluids, irrespective of their nature.

Rotational inertia-induced glassy transition in chiral particle systems.

The dense active matter exhibits characteristics reminiscent of traditional glassy phenomena, yet the role of rotational inertia in glass dynamics remains elusive. In this study, we investigate the glass dynamics of chiral active particles influenced by rotational inertia. Rotational inertia endows exponential memory to particle orientation, restricting its alteration and amplifying the effective persistence time. At lower spinning frequencies, the diffusion coefficient exhibits a peak function relative to rotational inertia for shorter persistence times, while it steadily increases with rotational inertia for longer persistence times. In the realm of high-frequency spinning, the impact of rotational inertia on diffusion behavior becomes more pronounced, resulting in a nonmonotonic and intricate relationship between the diffusion coefficient and rotational inertia. Consequently, the introduction of rotational inertia significantly alters the glassy dynamics of chiral active particles, allowing for the control over transitions between fluid and glassy states by modulating rotational inertia. Moreover, our findings indicate that at a specific spinning temperature, there exists an optimal spinning frequency at which the diffusion coefficient attains its maximum value.

Investigation of the kinematics behavior of spherical detector with different diameter ratios in oil pipelines

The spherical detector is an important technical equipment for water supply pipeline leakage detection. Its stable operation inside the pipeline will be affected by the fluid state in the pipeline, and the motion law of detectors with different diameters is different. In this paper, the dynamic analysis is carried out using the numerical simulation method, and the motion of spherical detectors with different diameters in horizontal oil pipelines is simulated. It is found that the flow state of the detector changes from the acceleration stage to the stable stage during the motion process, and with the decrease of the diameter ratio, the pressure difference between the front and back of the detector gradually increases, the thrust shows an upward trend, and the speed shows a greater speed in the stable stage. The motion behavior of the detector under different pipe diameter ratios is obtained, which provides help for the subsequent optimization of the shape design of the micro-leakage internal spherical detector.

Analysis of the number of topological defects in active nematic fluids under applied shear flow.

The number of topological defects in the shear flow of active nematic fluids is numerically investigated in this study. The evolution of the flow state of extensile active nematic fluids is explored by increasing the activity of active nematic fluids. Evidently, medium-activity active nematic fluids exhibit a highly ordered vortex lattice fluid state. However, high-activity active nematic fluids exhibit a meso-scale turbulent flow accompanied by topological defects. The number of topological defects (Ndef) increases with increasing shear Reynolds number (Res). Fluid viscosity strongly influences Ndef, while the influence of fluid density is relatively weak. Ndef decreases with increasing activity length scale (lζ) value. A small Res value strongly influences Ndef, whereas a large lζ value only weakly influences Ndef. As the activity increases, Ndef in contractile active nematic fluids becomes larger than that of extensile active nematic fluids.

The Chain of Processes Forming Porphyry Copper Deposits—An Invited Paper

Abstract Porphyry-related mineral deposits are giant geochemical anomalies in the Earth’s crust with orders-of-magnitude differences in the content and proportion of the three main ore metals Cu, Au, and Mo. Deposit formation a few kilometers below surface is the product of a chain of geologic processes operating at different scales in space and time. This paper explores each process in this chain with regard to optimizing the chances of forming these rare anomalies. On the lithosphere scale, deposits with distinct metal ratios occur in provinces that formed during brief times of change in plate motions. Similar metal ratios of several deposits in such provinces compared with global rock reservoirs suggest preceding enrichment of Au or Mo in lithospheric regions giving rise to distinct ore provinces. The largest Cu-dominated deposits and provinces are traditionally explained by selective removal of Au during generation or subsequent evolution of mantle magmas, but the possibility of selective Cu pre-enrichment of lithosphere regions by long-term subduction cannot be dismissed, even though its mechanism remains speculative. Evolution of hydrous basaltic melts to fertile magmas forming porphyry Cu deposits requires fractionation toward more H2O-rich magmas in the lower crust, as shown by their adakite-like trace element composition. The prevailing interpretation that this fractionation leads to significant loss of chalcophile ore metals by saturation and removal of magmatic sulfide might be inverted to a metal enrichment step, if the saturating sulfides are physically entrained with the melt fraction of rapidly ascending magmas. Ascent of fertile magma delivers a large mass of H2O-rich ore fluid to the upper crust, along points of weakness in an overall compressive stress regime, within a limited duration as required by mass and heat balance constraints. Two mechanisms of rapid magma ascent are in debate: (1) wholesale emplacement of highly fractionated and volatile-rich granitic melt into a massive transcrustal channelway, from which fluids are exsolved by decompression starting in the lower crust, or (2) partly fractionated magmas filling a large upper crustal magma chamber, from which fluids are expelled by cooling and crystallization. Transfer of ore-forming components to a hydrothermal ore fluid is optimized if the first saturating fluid is dense and Cl rich. This can be achieved by fluid saturation at high pressure, or after a moderately H2O rich intermediate-composition melt further crystallizes in an upper crustal reservoir before reaching fluid saturation. In either case, metals and S (needed for later hydrothermal sulfide precipitation) are transferred to the fluid together, no matter whether ore components are extracted from the silicate melt or liberated to the ore fluid by decomposition of magmatic sulfides. Production and physical focusing of fluids in a crystallizing upper crustal magma chamber are controlled by the rate of heat loss to surrounding rocks. Fluid focusing, requiring large-scale lateral flow, spontaneously occurs in mushy magma because high water content and intermediate melt/crystal ratio support a network of interconnected tubes at the scale of mineral grains. Calculated cooling times of such fluid-producing magma reservoirs agree with the duration of hydrothermal ore formation measured by high-precision zircon geochronology, and both relate to the size of ore deposits. Ore mineral precipitation requires controlled flow of S- and metal-rich fluids through a vein network, as shown by fluid inclusion studies. The degree of hydrothermal metal enrichment is optimized by the balance between fluid advection and the efficiency of cooling of the magmatic fluid plume by heat loss to convecting meteoric water. The depth of fluid production below surface controls the pressure-temperature (P-T) evolution along the upflow path of magmatic fluids. Different evolution paths controlling density, salinity, and phase state of fluids contribute to selective metal precipitation: porphyry Au deposits can form at shallow subvolcanic levels from extremely saline brine or salt melt; high-grade Au-Cu coprecipitation from coexisting and possibly rehomogenizing brine and vapor is most efficient at a depth of a few kilometers; whereas fluids cooling at greater depth tend to precipitate Cu ± Mo but transport Au selectively to shallower epithermal levels. Exhumation and secondary oxidation and enrichment by groundwater finally determine the economics of a deposit, as well as the global potential of undiscovered metal resources available for future mining.

Research on the reconstruction of porous bronze structures based on powder metallurgy simulation

The microscopic morphology of pores within the porous throttle determines the gas flow state, which directly affects the macroscopic performance of the air bearing. Geometric modeling of porous materials is a prerequisite for pore characterization and fluid state analysis. Therefore, in this study, a porous model was established by simulating each step of the powder metallurgy process of porous bronze, in which the sintering fusion simulation was innovatively accomplished by median filtering of the material cross-sectional images. The feasibility of the reconstruction method was verified by comparing pore morphology parameters and permeability results with the actual porous materials. On this basis, the flow simulation was performed to reveal the fluid flow state inside the pore channel structures. This method avoids the complex algorithms of mathematical methods and the dependence on high-precision instruments of experimental methods.

Thermalization and Hydrodynamics in an Interacting Integrable System: The Case of Hard Rods

We consider the relaxation of an initial non-equilibrium state in a one-dimensional fluid of hard rods. Since it is an interacting integrable system, we expect it to reach the Generalized Gibbs Ensemble (GGE) at long times for generic initial conditions. Here we show that there exist initial conditions for which the system does not reach GGE even at very long times and in the thermodynamic limit. In particular, we consider an initial condition of uniformly distributed hard-rods in a box with the left half having particles with a singular velocity distribution (all moving with unit velocity) and the right half particles in thermal equilibrium. We find that the density profile for the singular component does not spread to the full extent of the box and keeps moving with a fixed effective speed at long times. We show that such density profiles can be well described by the solution of the Euler equations almost everywhere except at the location of the shocks, where we observe slight discrepancies due to dissipation arising from the initial fluctuations of the thermal background. To demonstrate this effect of dissipation analytically, we consider a second initial condition with a single particle at the origin with unit velocity in a thermal background. We find that the probability distribution of the position of the unit velocity quasi-particle has diffusive spreading which can be understood from the solution of the Navier–Stokes (NS) equation of the hard rods. Finally, we consider an initial condition with a spread in velocity distribution for which we show convergence to GGE. Our conclusions are based on molecular dynamics simulations supported by analytical arguments.