Comparative study on phase transition behaviors of fractional molecular field theory and random-site Ising model

Fractional molecular field theory(FMFT)is a phenomenological theory that describes phase transitions in crystals with randomly distributed components,such as the relaxor-ferroelectrics and spin glasses.In order to verify the feasibility of this theory,this paper fits it to the Monte Carlo simulations of specific heat and susceptibility versus temperature of two-dimensional(2D)random-site Ising model(2D-RSIM).The results indicate that the FMFT deviates from the 2D-RSIM significantly.The main reason for the deviation is that the 2D-RSIM is a typical system of component random distribution,where the real order parameter is spatially heterogeneous and has no symmetry of space translation,but the basic assumption of FMFT means that the parameter is spatially uniform and has symmetry of space translation.

Emergent topological ordered phase for the Ising-XY model revealed by cluster-updating Monte Carlo method

The two-component cold atom systems with anisotropic hopping amplitudes can be phenomenologically described by a two-dimensional Ising-XY coupled model with spatial anisotropy.At low temperatures,theoretical predictions[Phys.Rev.A 72053604(2005)]and[arXiv:0706.1609]indicate the existence of a topological ordered phase characterized by Ising and XY disorder but with 2XY ordering.However,due to ergodic difficulties faced by Monte Carlo methods at low temperatures,this topological phase has not been numerically explored.We propose a linear cluster updating Monte Carlo method,which flips spins without rejection in the anisotropy limit but does not change the energy.Using this scheme and conventional Monte Carlo methods,we succeed in revealing the nature of topological phases with half-vortices and domain walls.In the constructed global phase diagram,Ising and XY-type transitions are very close to each other and differ significantly from the schematic phase diagram reported earlier.We also propose and explore a wide range of quantities,including magnetism,superfluidity,specific heat,susceptibility,and even percolation susceptibility,and obtain consistent and reliable results.Furthermore,we observed first-order transitions characterized by common intersection points in magnetizations for different system sizes,as opposed to the conventional phase transition where Binder cumulants of various sizes share common intersections.The critical exponents of different types of phase transitions are reasonably fitted.The results are useful to help cold atom experiments explore the half-vortex topological phase.

Bifurcation Analysis Reveals Solution Structures of Phase Field Models

The phase field method is playing an increasingly important role in understanding and predicting morphological evolution in materials and biological systems.Here,we develop a new analytical approach based on the bifurcation analysis to explore the mathematical solution structure of phase field models.Revealing such solution structures not only is of great mathematical interest but also may provide guidance to experimentally or computationally uncover new morphological evolution phenomena in materials undergoing electronic and structural phase transitions.To elucidate the idea,we apply this analytical approach to three representative phase field equations:the Allen-Cahn equation,the Cahn-Hilliard equation,and the Allen-Cahn-Ohta-Kawasaki system.The solution structures of these three phase field equations are also verified numerically by the homotopy continuation method.

Properties of the phase diagram from the Nambu-Jona-Lasino model with a scalar-vector interaction

We investigated the properties of the phase diagram of high-order susceptibilities,speed of sound,and polytropic index based on an extended Nambu-Jona-Lasinio model with an eight-quark scalar-vector interaction.Non-monotonic behavior was observed in all these quantities around the phase transition boundary,which also revealed the properties of the critical point.Further,this study indicated that the chiral phase transition boundary and critical point could vary depending on the scalarvector coupling constant G_(SV).At finite densities and temperatures,the negative G_(SV)term exhibited attractive interactions,which enhanced the critical point temperature and reduced the chemical potential.The G_(SV)term also affected the properties of the high-order susceptibilities,speed of sound,and polytropic index near the critical point.The non-monotonic(peak or dip)structures of these quantities shifted to a low baryon chemical potential(and high temperature)with a negative G_(SV).G_(SV)also changed the amplitude and range of the nonmonotonic regions.Therefore,the scalar-vector interaction was useful for locating the phase boundary and critical point in QCD phase diagram by comparing the experimental data.The study of the non-monotonic behavior of high-order susceptibilities,speed of sound,and polytropic index is of great interest,and further observations related to high-order susceptibilities,speed of sound,and polytropic index being found and applied to the search for critical points in heavy-ion collisions and the study of compact stars are eagerly awaited.

Modeling of strain-inducedγ→α′phase transformation for 201Cu metastable austenitic stainless steel and its verification in rolling processes

To model the strain-inducedγ→α′phase transformation for the Cr-Mn metastable austenitic stainless steel,the 201Cu steel was chosen as the analytical material and the cylindrical samples of this steel with size ofϕ5 mm×10 mm were compressed at strains of 0.2–0.6 in the temperature range of 25–150°C and in the strain rate range of 0.1–5.0 s^(−1).The flaky samples were prepared by wire cutting from the cylindrical samples and the volume fraction of the strain-inducedα′phase was detected in the test point of the flaky samples.The volume fraction changing with the process parameters was modeled,and the critical temperatures and the critical strains to preventγ→α′phase transformation were calculated as other different process parameters changed.The linear fitting goodness of the model between the calculated volume fraction values and the tested ones is 0.986 and the validity of the model was verified by application in cold and warm rolling experiments.

Phase field model for electric-thermal coupled discharge breakdown of polyimide nanocomposites under high frequency electrical stress

In contrast to conventional transformers, power electronic transformers, as an integral component of new energy power system, are often subjected to high-frequency and transient electrical stresses, leading to heightened concerns regarding insulation failures. Meanwhile, the underlying mechanism behind discharge breakdown failure and nanofiller enhancement under high-frequency electrical stress remains unclear. An electric-thermal coupled discharge breakdown phase field model was constructed to study the evolution of the breakdown path in polyimide nanocomposite insulation subjected to high-frequency stress. The investigation focused on analyzing the effect of various factors, including frequency, temperature, and nanofiller shape, on the breakdown path of Polyimide(PI) composites. Additionally, it elucidated the enhancement mechanism of nano-modified composite insulation at the mesoscopic scale. The results indicated that with increasing frequency and temperature, the discharge breakdown path demonstrates accelerated development, accompanied by a gradual dominance of Joule heat energy. This enhancement is attributed to the dispersed electric field distribution and the hindering effect of the nanosheets. The research findings offer a theoretical foundation and methodological framework to inform the optimal design and performance management of new insulating materials utilized in high-frequency power equipment.

Fracture propagation laws of staged hydraulic fracture in fractured geothermal reservoir based on phase field model

Hydraulic fracturing is widely used in geothermal resource exploitation, and many natural fractures exist in hot dry rock reservoirs due to in-situ stress and faults. However, the infuence of natural fractures on hydraulic fracture propagation is not considered in the current study. In this paper, based on the phase feld model, a thermo-hydro-mechanical coupled hydraulic fracture propagation model was established to reveal the infuence of injection time, fracturing method, injection fow rate, and natural fracture distribution on the fracture propagation mechanism. The results show that fracture complexity increases with an increase in injection time. The stress disturbance causes the fracture initiation pressure of the second cluster signifcantly higher than that of the frst and third clusters. The zipper-type fracturing method can reduce the degree of stress disturbance and increase fracture complexity by 7.2% compared to simultaneous hydraulic fracturing. Both low and high injection fow rate lead to a decrease in fracture propagation time, which is not conducive to an increase in fracture complexity. An increase in the natural fracture angle leads to hydraulic fracture crossing natural fracture, but has a lesser efect on fracture complexity. In this paper, we analyzed the infuence of diferent factors on initiation pressure and fracture complexity, providing valuable guidance for the exploitation of geothermal resources.

First-order quantum phase transition and entanglement in the Jaynes-Cummings model with a squeezed light

We study the quantum phase transition and entanglement in the Jaynes-Cummings model with squeezed light,utilize a special transformation method to obtain the analytical ground state of the model within the near-resonance regime,and numerically verify the validity of the analytical ground state.It is found that the ground state exhibits a first-order quantum phase transition at the critical point linearly induced by squeezed light,and the ground state entanglement reaches its maximum when the qubit-field coupling strength is large enough at the critical point.

Multi-dimensional Simulation of Phase Change by a 0D-2D Model Coupling via Stefan Condition

Considering phase changes associated with a high-temperature molten material cooled down from the outside,this work presents an improvement of the modelling and the numerical simulation of such processes for an application pertaining to the safety of light water nuclear reactors.Postulating a core meltdown accident,the behaviour of the core melt(aka corium)into a steel vessel is of tremendous importance when evaluating the vessel integrity.Evaluating correctly the heat fluxes requires the numerical simulation of the interaction between the liquid material and its solid counterpart which forms during the solidification process,but also may melt back.To simulate this configuration,encoun-tered in various industrial applications,one considers a bi-phase model constituted by a liquid phase in contact and interaction with its solid phase.The liquid phase may solidify in presence of low energetic source,while the solid phase may melt due to an intense heat flux from the high-temperature liquid.In the frame of the in-house legacy code,several simplifying assumptions(0D multi-layer discretization,instantaneous heat transfer via a quadratic temperature profile in solids)are made for the modelling of such phase changes.In the present work,these shortcomings are illustrated and further overcome by solving a 2D heat conduction model in the solid by a mixed Raviart-Thomas finite element method coupled to the liquid phase due to heat and mass exchanges through Stefan condition.The liquid phase is modeled with a 0D multi-layer approach.The 0D-liquid and 2D-solid mod-els are coupled by a Stefan like phase change interface model.Several sanity checks are performed to assess the validity of the approach on 1D and 2D academical configurations for which exact or reference solutions are available.Then more advanced situations(genu-ine multi-dimensional phase changes and an"industrial-like scenario")are simulated to verify the appropriate behavior of the obtained coupled simulation scheme.

Ground-state phase diagram,symmetries,excitation spectra and finite-frequency scaling of the two-mode quantum Rabi model

We investigate the two-mode quantum Rabi model(QRM)describing the interaction between a two-level atom and a two-mode cavity field.The quantum phase transitions are found when the ratioηof transition frequency of atom to frequency of cavity field approaches infinity.We apply the Schrieffer–Wolff(SW)transformation to derive the low-energy effective Hamiltonian of the two-mode QRM,thus yielding the critical point and rich phase diagram of quantum phase transitions.The phase diagram consists of four regions:a normal phase,an electric superradiant phase,a magnetic superradiant phase and an electromagnetic superradiant phase.The quantum phase transition between the normal phase and the electric(magnetic)superradiant phase is of second order and associates with the breaking of the discrete Z_(2) symmetry.On the other hand,the phase transition between the electric superradiant phase and the magnetic superradiant phase is of first order and relates to the breaking of the continuous U(1)symmetry.Several important physical quantities,for example the excitation energy and average photon number in the four phases,are derived.We find that the excitation spectra exhibit the Nambu–Goldstone mode.We calculate analytically the higher-order correction and finite-frequency exponents of relevant quantities.To confirm the validity of the low-energy effective Hamiltonians analytically derived by us,the finite-frequency scaling relation of the averaged photon numbers is calculated by numerically diagonalizing the two-mode quantum Rabi Hamiltonian.

Numerical simulation of microwave-induced cracking and melting of granite based on mineral microscopic models

This study introduces a coupled electromagnetic–thermal–mechanical model to reveal the mechanisms of microcracking and mineral melting of polymineralic rocks under microwave radiation.Experimental tests validate the rationality of the proposed model.Embedding microscopic mineral sections into the granite model for simulation shows that uneven temperature gradients create distinct molten,porous,and nonmolten zones on the fracture surface.Moreover,the varying thermal expansion coefficients and Young's moduli among the minerals induce significant thermal stress at the mineral boundaries.Quartz and biotite with higher thermal expansion coefficients are subjected to compression,whereas plagioclase with smaller coefficients experiences tensile stress.In the molten zone,quartz undergoes transgranular cracking due to theα–βphase transition.The local high temperatures also induce melting phase transitions in biotite and feldspar.This numerical study provides new insights into the distribution of thermal stress and mineral phase changes in rocks under microwave irradiation.

Topological phases and edge modes of an uneven ladder

We investigate the topological properties of a two-chain quantum ladder with uneven legs,i.e.,the two chains differ in their periods by a factor of 2.Such an uneven ladder presents rich band structures classified by the closure of either direct or indirect bandgaps.It also provides opportunities to explore fundamental concepts concerning band topology and edge modes,including the difference of intracellular and intercellular Zak phases,and the role of the inversion symmetry(IS).We calculate the Zak phases of the two kinds and find excellent agreement with the dipole moment and extra charge accumulation.We also find that configurations with IS feature a pair of degenerate two-side edge modes emerging as the closure of the direct bandgap,while configurations without IS feature one-side edge modes emerging as not only the closure of both direct and indirect bandgaps but also within the band continuum.Furthermore,by projecting to the two sublattices,we find that the effective Bloch Hamiltonian corresponds to that of a generalized Su–Schrieffer–Heeger model or the Rice–Mele model whose hopping amplitudes depend on the quasimomentum.In this way,the topological phases can be efficiently extracted through winding numbers.We propose that uneven ladders can be realized by spin-dependent optical lattices and their rich topological characteristics can be examined by near future experiments.

Phase equilibrium data prediction and process optimizationin butadiene extraction process

In response to the lack of reliable physical parameters in the process simulation of the butadiene extraction,a large amount of phase equilibrium data were collected in the context of the actual process of butadiene production by acetonitrile.The accuracy of five prediction methods,UNIFAC(UNIQUAC Functional-group Activity Coefficients),UNIFAC-LL,UNIFAC-LBY,UNIFAC-DMD and COSMO-RS,applied to the butadiene extraction process was verified using partial phase equilibrium data.The results showed that the UNIFAC-DMD method had the highest accuracy in predicting phase equilibrium data for the missing system.COSMO-RS-predicted multiple systems showed good accuracy,and a large number of missing phase equilibrium data were estimated using the UNIFAC-DMD method and COSMO-RS method.The predicted phase equilibrium data were checked for consistency.The NRTL-RK(non-Random Two Liquid-Redlich-Kwong Equation of State)and UNIQUAC thermodynamic models were used to correlate the phase equilibrium data.Industrial device simulations were used to verify the accuracy of the thermodynamic model applied to the butadiene extraction process.The simulation results showed that the average deviations of the simulated results using the correlated thermodynamic model from the actual values were less than 2%compared to that using the commercial simulation software,Aspen Plus and its database.The average deviation was much smaller than that of the simulations using the Aspen Plus database(>10%),indicating that the obtained phase equilibrium data are highly accurate and reliable.The best phase equilibrium data and thermodynamic model parameters for butadiene extraction are provided.This improves the accuracy and reliability of the design,optimization and control of the process,and provides a basis and guarantee for developing a more environmentally friendly and economical butadiene extraction process.

Nonlinear wave dispersion in monoatomic chains with lumped and distributed masses:discrete and continuum models

The main objective of this paper is to investigate the influence of inertia of nonlinear springs on the dispersion behavior of discrete monoatomic chains with lumped and distributed masses.The developed model can represent the wave propagation problem in a non-homogeneous material consisting of heavy inclusions embedded in a matrix.The inclusions are idealized by lumped masses,and the matrix between adjacent inclusions is modeled by a nonlinear spring with distributed masses.Additionally,the model is capable of depicting the wave propagation in bi-material bars,wherein the first material is represented by a rigid particle and the second one is represented by a nonlinear spring with distributed masses.The discrete model of the nonlinear monoatomic chain with lumped and distributed masses is first considered,and a closed-form expression of the dispersion relation is obtained by the second-order Lindstedt-Poincare method(LPM).Next,a continuum model for the nonlinear monoatomic chain is derived directly from its discrete lattice model by a suitable continualization technique.The subsequent use of the second-order method of multiple scales(MMS)facilitates the derivation of the corresponding nonlinear dispersion relation in a closed form.The novelties of the present study consist of(i)considering the inertia of nonlinear springs on the dispersion behavior of the discrete mass-spring chains;(ii)developing the second-order LPM for the wave propagation in the discrete chains;and(iii)deriving a continuum model for the nonlinear monoatomic chains with lumped and distributed masses.Finally,a parametric study is conducted to examine the effects of the design parameters and the distributed spring mass on the nonlinear dispersion relations and phase velocities obtained from both the discrete and continuum models.These parameters include the ratio of the spring mass to the lumped mass,the nonlinear stiffness coefficient of the spring,and the wave amplitude.

Study on theoretical model for electrical explosion resistivity of Al/Ni reactive multilayer foil

Al/Ni reactive multilayer foil(RMF)possesses excellent comprehensive properties as a promising substitute for traditional Cu bridge.A theoretical resistivity model of Al/Ni RMF was developed to guide the optimization of EFIs.Al/Ni RMF with different bilayer thicknesses and bridge dimensions were prepared by MEMS technology and electrical explosion tests were carried out.According to physical and chemical reactions in bridge,the electrical explosion process was divided into 5 stages:heating of condensed bridge,vaporization and diffusion of Al layers,intermetallic combination reaction,intrinsic explosion,ionization of metal gases,which are obviously shown in measured voltage curve.Effects of interface and grain boundary scattering on the resistivity of film metal were considered.Focusing on variations of substance and state,the resistivity was developed as a function of temperature at each stage.Electrical explosion curves were calculated by this model at different bilayer thicknesses,bridge dimensions and capacitor voltages,which showed an excellent agreement with experimental ones.

CMIP 6 models simulation of the connection between North/South Pacific Meridional Mode and ENSO

The subtropical North and South Pacific Meridional Modes(NPMM and SPMM)are well known precursors of El Niño-Southern Oscillation(ENSO).However,relationship between them is not constant.In the early 1980,the relationship experienced an interdecadal transition.Changes in this connection can be attributed mainly to the phase change of the Pacific decadal oscillation(PDO).During the positive phase of PDO,a shallower thermocline in the central Pacific is responsible for the stronger trade wind charging(TWC)mechanism,which leads to a stronger equatorial subsurface temperature evolution.This dynamic process strengthens the connection between NPMM and ENSO.Associated with the negative phase of PDO,a shallower thermocline over southeastern Pacific allows an enhanced wind-evaporation-SST(WES)feedback,strengthening the connection between SPMM and ENSO.Using 35 Coupled Model Intercomparison Project Phase 6(CMIP6)models,we examined the NPMM/SPMM performance and its connection with ENSO in the historical runs.The great majority of CMIP6 models can reproduce the pattern of NPMM and SPMM well,but they reveal discrepant ENSO and NPMM/SPMM relationship.The intermodal uncertainty for the connection of NPMM-ENSO is due to different TWC mechanism.A stronger TWC mechanism will enhance NPMM forcing.For SPMM,few models can simulate a good relationship with ENSO.The intermodel spread in the relationship of SPMM and ENSO owing to SST bias in the southeastern Pacific,as WES feedback is stronger when the southeastern Pacific is warmer.

Mathematical Wave Functions and 3D Finite Element Modelling of the Electron and Positron

The wave/particle duality of particles in Physics is well known. Particles have properties that uniquely characterize them from one another, such as mass, charge and spin. Charged particles have associated Electric and Magnetic fields. Also, every moving particle has a De Broglie wavelength determined by its mass and velocity. This paper shows that all of these properties of a particle can be derived from a single wave function equation for that particle. Wave functions for the Electron and the Positron are presented and principles are provided that can be used to calculate the wave functions of all the fundamental particles in Physics. Fundamental particles such as electrons and positrons are considered to be point particles in the Standard Model of Physics and are not considered to have a structure. This paper demonstrates that they do indeed have structure and that this structure extends into the space around the particle’s center (in fact, they have infinite extent), but with rapidly diminishing energy density with the distance from that center. The particles are formed from Electromagnetic standing waves, which are stable solutions to the Schrödinger and Classical wave equations. This stable structure therefore accounts for both the wave and particle nature of these particles. In fact, all of their properties such as mass, spin and electric charge, can be accounted for from this structure. These particle properties appear to originate from a single point at the center of the wave function structure, in the same sort of way that the Shell theorem of gravity causes the gravity of a body to appear to all originate from a central point. This paper represents the first two fully characterized fundamental particles, with a complete description of their structure and properties, built up from the underlying Electromagnetic waves that comprise these and all fundamental particles.

Thermodynamic consistent phase field model for sintering process with multiphase powders

A thermodynamic consistent phase field model is developed to describe the sintering process with multiphase powders. In this model, the interface region is assumed to be a mixture of different phases with the same chemical potential, but with different compositions. The interface diffusion and boundary diffusion are also considered in the model. As an example, the model is applied to the sintering process with Fe-Cu powders. The free energy of each phase is described by the well-developed thermodynamic models, together with the published optimized parameters. The microstructure and solute distribution during the sintering process can both be obtained quantitively.

Simulation of pre-precipitation in Ni_(75)Al_(14)Mo_(11) alloy by microscopic phase-field model

The early precipitation process of Ni（75）Al（14）Mo（11） alloy was simulated by microscopic phase-field model at different temperatures.The microstructure of the alloy,the precipitation time of Llo structure and occupation probability of the three kinds of atoms were investigated.It is indicated that the non-stoichiometric Ll0（Ⅰ/Ⅱ） phases are found in the precipitation process.With the temperature increasing,the appearance time of Ll0 is brought forward.The Ll0（Ⅱ） structure always precipitates earlier than the Ll0（Ⅰ） structure.Compared with lower temperature,higher temperature brings the formation time of Ll0 phase forward and makes Ll0 phase have a higher order degree.But lower temperature shortens the process time of the Ll0 phase to the Ll2 phase.Al and Mo atoms tend to occupy γ site,Ni atom tends to occupy a and β sites.At the same temperature,Al atom has stronger occupation ability than Mo atom in the same site.Ni,Al and Mo collectively form the composited Ll2 structure.

Three-dimensional modeling of effect of surface intermetallic phase on surface defects of Al-Fe-Si aluminum foils during twin-roll casting

Scanning electron microscopy and X-ray energy dispersive spectrum analysis show that the clusters of intermetallic AlFeSi particle are distributed on or near the aluminum foil stock surfaces heterogeneously. 3D finite element modeling shows that these clusters of hard particles induce the fracture of the nano-scale lubricant oil film at first and further lead to severe deformation in the nearby aluminum foil substrate along the rolling direction. Consequently, the optical property in this region differs from that in the surroundings, resulting in surface defects.