Utility perceived by individuals is believed to be different from the utility experienced by that individual. System dynamicists implicitly categorize this phenomenon as a form of bounded rationality and traditionally employ a simple smoothing function to capture it. We challenge this generalization by testing it against an alternative formulation of utility perception that is suggested by modern theories of behavioral economics. In particular, the traditional smoothing formulation is compared with the peak-end rule in a simple theoretical model as well as in a medium-size model of electronic health record implementation. Experimentation with the models reveals that the way utility perception is formulated is important and might affect behavior and policy implications of system dynamics models.

Polarization in online social networks has gathered a significant amount of attention in the research community and in the public sphere due to stark disagreements with millions of participants in topics surrounding politics, climate, the economy and other areas where an agreement is required. There are multiple approaches to investigating the scenarios in which polarization occurs and given that polarization is not a new phenomenon but that its virality may be supported by the low cost and latency messaging offered by online social media platforms; an investigation into the intrinsic dynamics of online opinion evolution is presented for complete networks. Extending a model which utilizes the Binary Voter Model (BVM) to examine the effect of the degree of freedom for selecting contacts based upon homophily, simulations show that different opinions are reinforced for a period of time when users have a greater range of choice for association. The facility of discussion threads and groups formed upon common views further delays the rate in which a consensus can form between all members of the network. This can temporarily incubate members from interacting with those who can present an alternative opinion where a voter model would then proceed to produce a homogeneous opinion based upon pairwise interactions.

After UV excitation, gas phase thymine returns to ground state in 5 to 7 ps, showing multiple time constants. There is no consensus on the assignment of these processes, with a dispute between models claiming that thymine is trapped either in the first (S1) or in the second (S2) excited states. In the present study, nonadiabatic dynamics simulation of thymine is performed on the basis of ADC(2) surfaces, to understand the role of dynamic electron correlation on the deactivation pathways. The results show that trapping in S2 is strongly reduced in comparison to previous simulations considering only non-dynamic electron correlation on CASSCF surfaces. The reason for the difference is traced back to the energetic cost for formation of a CO p bond in S2.

Complex systems can be described as the network of interactions between their components. In directed networks, information, energy and matter flow in dominant directions. Undirected networks, on the other hand, cannot easily capture these asymmetries. The local flow of information depends on the node spin, which is a function of the number of incoming and outgoing links attached to it. And yet, how the heterogeneous association between in-degree and out-degree constraints global dynamics remains poorly understood. Here, we develop a new theoretical model to study the relationship between node directionality and global flow dynamics. Our computational analysis of random and scale-free directed networks indicates that structural diversity constrains the dimensionality of reachable dynamical space. These results suggest how local directionality might be an universal driver of global dynamics in many systems, from social networks, to technological systems to the connectome.

Cancers remain the lead cause of disease-related, pediatric death in North America. The emerging field of complex systems has redefined cancer networks as a computational system with intractable algorithmic complexity. Herein, a tumor and its heterogeneous phenotypes are discussed as dynamical systems having multiple, strange attractors. Machine learning, network science and algorithmic information dynamics are discussed as current tools for cancer network reconstruction. Deep Learning architectures and computational fluid models are proposed for better forecasting gene expression patterns in cancer ecosystems. Cancer cell decision-making is investigated within the framework of complex systems and complexity theory.

Complex biological systems consist of large numbers of interconnected units, characterized by emergent properties such as collective computation. In spite of all the progress in the last decade, we still lack a deep understanding of how these properties arise from the coupling between the structure and dynamics. Here, we introduce the multiscale emergent functional state, which can be represented as a network where links encode the flow exchange between the nodes, calculated using diffusion processes on top of the network. We analyze the emergent functional state to study the distribution of the flow among components of 92 fungal networks, identifying their functional modules at different scales and, more importantly, demonstrating the importance of functional modules for information content of networks, quantified in terms of network spectral entropy. Our results suggest that the topological complexity of fungal networks guarantees the existence of functional modules at different scales keeping the information entropy, and functional diversity, high.

Cooperation is widely recognized to be challenging for the well-balanced development of human societies. The emergence of cooperation in populations has been largely studied in the context of the Prisoner's Dilemma game, where temptation to defect and fear to be betrayed by others often activate defective strategies. In this paper we analyze the decision making mechanisms fostering cooperation in the two-strategy Stag-Hunt and Chicken games, which include the mixed strategy Nash equilibrium, describing partially cooperative behavior. We find the conditions for which cooperation is asymptotically stable in both full and partial cases, and we show that the partially cooperative steady state is also globally stable in the simplex. Furthermore, we show that the last can be more rewarding than the first, thus making the mixed strategy effective, although people cooperate at a lower level with respect to the maximum allowed, as it is reasonably expected in real situations. Our findings highlight the importance of Stag-Hunt and Chicken games in understanding the emergence of cooperation in social networks.

In this essay, I revisit Neural Networks (NNs) through the lens of evolutionary dynamics. From the two most important features of NNs, I recover the two most general equations of evolutionary dynamics. These two equations may thus serve as a connection between NNs and a body of knowledge that has been built over multiple centuries and has found application across a vast range of disciplines. This connection may therefore also help explain why NNs with the two features in question are applicable across a much broader range of domains than initially envisioned.

This paper investigates the dynamic properties of artificial neural networks using differential equations and explores the influence of parameters on stability and neural oscillations. By analyzing the equilibrium point of the differential equations, we identify conditions for asymptotic stability and criteria for oscillation in artificial neural networks. Furthermore, we demonstrate how adjusting synaptic weights between neurons can effectively control stability and oscillation. The proposed model offers potential insights into the malfunctioning mechanisms of biological neural networks implicated in neurological disorders like Parkinson's disease tremors and epilepsy seizures, which are characterized by abnormal oscillations.

Given the high number of vehicle-crash victims, it has been established as a priority to reduce this figure in the transportation sector. For this reason, many of the recent researches are focused on including control systems in existing vehicles, to improve their stability, comfort and handling. These systems need to know in every moment the behavior of the vehicle (state variables), among others, when the different maneuvers are performed, to actuate by means of the systems in the vehicle (brakes, steering, suspension) and, in this way, to achieve a good behavior. The main problem arises from the lack of ability to directly capture several required dynamic vehicle variables, such as roll angle, from low-cost sensors. Previous studies demonstrate that low-cost sensors can provide data in real-time with the required precision and reliability. Even more, other research works indicate that neural networks are efficient mechanisms to estimate roll angle. Nevertheless, it is necessary to assess that the fusion of data coming from low-cost devices and estimations provided by neural networks can fulfill the reliability and appropriateness requirements for using these technologies to improve overall safety in production vehicles. Because of the increasing of computing power, the reduction of consumption and electric devices size, along with the high variety of communication technologies and networking protocols using Internet have yield to Internet of Things (IoT) development. In order to address this issue, this study has two main goals: 1) Determine the appropriateness and performance of neural networks embedded in low-cost sensors kits to estimate roll angle required to evaluate rollover risk situations. 2) Compare the low-cost control unit devices (Intel Edison and Raspberry Pi 3 Model B), to provide the roll angle estimation with this artificial neural network-based approach. To fulfil these objectives an experimental environment has been set up composed of a van with two set of low-cost kits, one including a Raspberry Pi 3 Model B, low cost Inertial Measurement Unit (BNO055 - 37€) and GPS (Mtk3339 - 53€) and the other having an Intel Edison System on Chip linked to a SparkFun 9 Degrees of Freedom module. This experimental environment will be tested in different maneuvers for comparison purposes. Neural networks embedded in low-cost sensor kits provide roll angle estimations very approximated to real values. Even more, Intel Edison and Raspberry Pi 3 Model B have enough computing capabilities to successfully run roll angle estimation based on neural networks to determine rollover risks situation fulfilling real-time operation restrictions stated for this problem.

Drawing inspiration from a recent construction of a polyhedral structure associated with an icosahedrally-symmetric map on the Riemann sphere, the article shows how to build such "dynamical polyhedra" for other icosahedral maps. First, icosahedral algebra is used to determine a special family of maps with 60 periodic critical points. The topological behavior of each map is worked out and results in a geometric algorithm that constructs a system of edges---the dynamical polyhedron---in natural correspondence to a map's topology. It turns out that the maps' descriptions fall into classes the presentation of which concludes the paper.

Evolution of cooperation has traditionally been studied by assuming that individuals adopt either of two pure strategies, to cooperate or defect. Recent work have considered continuous cooperative investments, turning full cooperation and full defection into two opposing ends of a spectrum and sometimes allowing for the emergence of the traditionally-studied pure strategies through evolutionary diversification. These studies have typically assumed a well-mixed population in which individuals are encountered with equal probability, Here, we allow for the possibility of assortative interactions by assuming that, with specified probabilities, an individual interacts with one or more other individuals of the same strategy. A closely related assumption has previously been made in evolutionary game theory and has been interpreted in terms of relatedness. We systematically study the effect of relatedness and find, among other conclusions, that the scope for evolutionary branching is reduced by either higher average degree of, or higher uncertainty in, relatedness with interaction partners. We also determine how different types of non-linear dependencies of benefits and costs constrain the types of evolutionary outcomes that can occur. While our results overall corroborate the conclusions of earlier studies, that higher relatedness promotes the evolution of cooperation, our investigation gives a comprehensive picture of how relatedness affects the evolution of cooperation with continuous investments.

We research the interaction between six representative carbon-based nanoparticles (CBNs) and 20 standard amino acids through all-atom molecular dynamics simulations. The six carbon-based nanoparticles are fullerene(C60), CNT55L3, CNT1010L3, CNT1515L3, CNT2020L3, and two-dimensional-graphene(Graphene33). Their curvatures decrease sequentially, and all of CNT are single-walled carbon nanotubes. We have observed that as the curvature of CBNs decreases, the adsorption effect of 20 amino acids with them has an increasing trend. In addition, we also used multi-dimensional clustering to analyze the adsorption effects of 20 amino acids on six carbon-based nanoparticles. We observed that the π-π interaction still plays an extremely important role in the adsorption of amino acids on carbon-based nanoparticles. Individual long-chain amino acids and “Benzene-like” Pro also have a strong adsorption effect with carbon-based nanoparticles.

Li₂MnO₃ is critical component in the well-studied Li-excess cathode materials xLi₂MnO₃•(1- x)LiMO₂ for achieving high lithium storage capacity. In this article, the diffusion of Li ions in Li₂MnO₃ is studied using molecular dynamics (MD) simulation with well-behaved empirical force fields obtained by fitting against the crystal structure from experiment and phonons calculated using Density Function Theory (DFT). We have found two possible tetrahedral hopping channels, 0-TM and 1-TM(Mn⁴⁺) channel, which are differentiated by the face sharing octahedral cations. Simulation results show that the 0-TM channel is active for Li hopping, while 1-TM(Mn4+) channel is inactive. During the delithiation process, the Li ions in the TM layered are firstly removed, then those in the Li layer. However, the Li ions will be trapped in the tetrahedral 0-TM channels as long as the four face sharing octahedral sites are cleared. Up to x=1.0 for Li₂-xMnO₃, almost all the Li ions are located at the tetrahedral sites, forming a regular array along a axis. The de-intercalation of tetrahedral Li ions requires a high voltage (>5.2 V vs. Li/Li+), limiting the practical capacities measured in lab. The diffusion of Mn ions into the Li layers is observed in a deeper delithiated structure (x=1.2 for Li₂-xMnO₃), indicating an initial phase transformation to a spinel-like structure. However, the Mn ions are mainly trapped in the tetrahedral sites in the Li layer, instead of the octahedral sites fin spinel-like structure. A few of Mn ions diffusing into the octahedral sites in Li layers have no face sharing tetrahedral Li ions, revealing a further Li de-intercalation is imperative for the complete phase transformation. Our model is not stable for x≥1.4 in Li₂-xMnO₃. Other charge compensation mechanism should be considered in this high delithiation stage, eg. oxygen release.

The selectivity in the simultaneous detection of ascorbic acid (AA), dopamine (DA), and uric acid (UA) has been an open problem in the biosensing field. Many surface modification methods were carried out for glassy carbon electrodes (GCE), including the use of graphene oxide and amino acids as a selective layer. In this work, molecular dynamics (MD) simulations were performed to investigate the role of serine oligomers on the selectivity of the AA, DA, UA analytes. Our models consisted of a graphene oxide (GO) sheet under a solvent environment. Serine tetramers were added into the simulation box and were adsorbed on the GO surface. Then, the adsorption of each analyte on the mixed surface was monitored from MD trajectories. It was found that the adsorption of AA was preferred by serine oligomers due to the largest number of hydrogen-bond forming functional groups of AA, while UA was the least preferred due to its highest aromaticity. Finally, the role of hydrogen bonds on the electron transfer selectivity of biosensors was discussed with some previous studies.

Quantum tunneling in a two-dimensional integrable map is studied. The orbits of the map are all confined to the curves specified by the one-dimensional Hamiltonian. It is found that the behavior of tunneling splitting for the integrable map and the associated Hamiltonian system is qualitatively the same, with only a slight difference in magnitude. However, the tunneling tails of the wave functions, obtained by superposing the eigenfunctions that form the doublet, exhibit significant difference. To explore the origin of the difference, we observe the classical dynamics in the complex plane and find that the existence of branch points appearing in the potential function of the integrable map could play the role for yielding non-trivial behavior in the tunneling tail. The result highlights the subtlety of quantum tunneling, which cannot be captured in nature only by the dynamics in the real plane.

The article presents the stages of modeling and simulation of a new design of a wheelchair with the option of moving up and down stairs. These analyzes were aimed at the synthesis of the de-sign parameters and parameters of the sensor and control systems. The simulation results were verified by experimentally testing the prototype.

The article is devoted to the issues of assessing the impact of changes in the delay value in the feed direction, what determines the regenerative nature of self-excitation of the cutting system, on the dynamics of metal processing on metal-cutting machines. It is this interconnected model of the cutting system that allows to identify the main dynamic effects that occur in the cutting system, including the effect of vibrations regeneration when cutting along the track. The article also has a section dedicated to the experiment on a real metalworking machine, on which vibrations of the top of the cutting wedge of the tool were measured and based on these data, a real calculation of the delay time and the magnitude of the change in the area of the cut layer was carried out, due to variations in the real feed. The conducted studies have shown that in addition to the vibrations of the cutting tool in the direction of feeding, the vibration activity of the tool in the direction of cutting plays an essential role in ensuring the regenerative effect. The delay time constant of the operator is formed, which determines the real value of the feed during cutting.

We present a new study of Jupiter atmosphere’s dynamics using for the first time the extremely high resolution capabilities of VLT/ESPRESSO to retrieve wind velocities in Jupiter’s troposphere, with a dedicated ground-based Doppler velocimetry method. These results are complemented by a deeper analysis of Cassini data during its flyby of Jupiter in December 2000, performing cloud tracking at visible wavelengths, obtaining a more comprehensive dynamical interpretation. We explore the effectiveness of this new method to measure winds in Jupiter, using high resolution spectroscopy data obtained from ESPRESSO observations performed in July 2019, with a Doppler velocimetry method based on back scattered solar radiation in the visible range. Coupled with our ground based results, we retrieved a latitudinal and longitudinal profile of Jupiter’s winds along select bands of the atmosphere. Comparing the results between cloud tracking methods, based on previous reference observations, and our new Doppler velocimetry approach we found a good agreement between them, demonstrating the eectiveness of this technique. The winds obtained in this exploratory study have a two-fold relevance: they contribute for the temporal and spatial variability study of Jupiter troposphere’s dynamics, and also the results presented here validate this Doppler technique to study the dynamics of Jupiter’s atmosphere and pave the way for further exploration of a broader region of Jupiter’s disk for a more comprehensive retrieval of winds and to evaluate their spatial and temporal variability.

Many species of plants are found in regions to which they are alien. Their global distributions are characterised by a family of exponential functions of the kind that arise in elementary statistical mechanics (an example in ecology is MacArthur's broken stick). We show here that all these functions are quantitatively reproduced by a model containing a single parameter – some global resource partitioned at random on the two axes of species number and site number. A dynamical model generating this equilibrium is a two fold stochastic process and suggests a curious and interesting biological interpretation in terms of niche structures fluctuating with time and productivity; with sites and species highly idiosyncratic. Idiosyncrasy implies that attempts to identify a priori those species likely to become naturalized are unlikely to be successful. Although this paper is primarily concerned with a particular problem in population biology, the two fold stochastic process may be of more general interest.

In the present work, we provide an electronic structure based method for the “on-the-fly” deter- mination of vibrational sum frequency generation (v-SFG) spectra. The predictive power of this scheme is demonstrated at the air-water interface. While the instantaneous fluctuations in dipole moment are obtained using the maximally localized Wannier functions, the fluctuations in polar- izability are approximated to be proportional to the second moment of Wannier functions. The spectrum henceforth obtained captures the signatures of hydrogen bond stretching, bending, as well as low-frequency librational modes.

Modelling errors, robust stabilization/tracking problems under parameter and model uncertainties complicate the control of the flexible underactuated systems. Chattering-free sliding-mode based input-output control law realizes robustness against the structured and unstructured uncertainties in the system dynamics and avoids excitation of unmodeled dynamics. The main purpose is to propose a robust adaptive solution for stabilizing and tracking direct-drive (DD) flexible robot arms under parameter and model uncertainties, as well as external disturbances. A lightweight robot arm subject to external and internal dynamic effects was taken into consideration. The challenges are compensating actuator dynamics with the inverter switching effects and torque ripples, stabilizing the zero dynamics under parameter/model uncertainties and disturbances while precisely track the predefined reference position. The precise control of this kind of system demands an accurate system model and knowledge of all sources that excite unmodeled dynamics. For this purpose, equations of motion for a flexible robot arm were derived and formulated for the large motion via Lagrange’s method. The goals were determined to achieve high-speed, precise position control, and satisfied accuracy by compensating the unwanted torque ripple and friction that degrades performance through an adaptive robust control approach. The actuator dynamics and their effect on the torque output were investigated due to the transmitted torque to the load side. The high-performance goals, precision&robustness issues, and stability concerns were satisfied by using robust-adaptive input-output linearization-based control law combining chattering-free sliding mode control (SMC) while avoiding the excitation of unmodeled dynamics.

The phenomenon of droplet impact on moving surfaces is widely observed in fields such as transportation, rotating machinery, and inkjet printing. Droplets exhibit non-axisymmetric behavior due to the motion of solid surfaces which significantly determines core parameters such as contact time, maximum spreading radius, and bounding velocity, thereby affecting the efficiency of related applications. In this study, we focus on kinetics and morphology of the non-axisymmetric bouncing behaviors for droplets impacting on a moving superhydrophobic surface (SHPS) within the normal (Wen) and tangential (Wet) Weber numbers. Considering the influences of moving surface on the contact area and contact time, the previous scaling formula for the horizontal velocity of droplets has been improved. Based on the velocity superposition hypothesis, we establish a theoretical model for the ratio of the maximum spreading radius at both ends depending on Wen and Wet. This research provides both experimental and theoretical evidence for understanding and controlling the non-axisymmetric behavior of droplets impacting on moving surfaces.

Single nucleotide polymorphisms (SNP) are associated with diseases and drug response variabilities in humans. Elucidating the damaging and disease-associated SNPs using wet-laboratory approaches can be challenging and resource-demanding due to the large number of SNPs in the human genome. Due to the growth in the field of computational biology and bioinformatics, algorithms have been developed to help screen and filter out the most deleterious SNPs that are worth considering for wet-laboratory studies. Here we review the existing in-silico based methods used to predict and characterize the effects of SNPs on protein structure and function. This cutting-edge approach will facilitate the search for novel therapeutics, help understand the etiology of diseases and fast-track the personalized medicine agenda.

Due to failure of the continuum hypothesis for higher Knudsen numbers, rarefied gases and microflows of gases are particularly difficult to model. Macroscopic transport equations compete with particle methods, such as DSMC to find accurate solutions in the rarefied gas regime. Due to growing interest in micro flow applications, such as micro fuel cells, it is important to model and understand evaporation in this flow regime. Here, evaporation boundary conditions for the R13 equations, which are macroscopic transport equations with applicability in the rarefied gas regime, are derived. The new equations utilize Onsager relations, linear relations between thermodynamic fluxes and forces, with constant coefficients, that need to be determined. For this, the boundary conditions are fitted to DSMC data and compared to other R13 boundary conditions from kinetic theory and Navier-Stokes-Fourier (NSF) solutions for two one-dimensional steady-state problems. Overall, the suggested fittings of the new phenomenological boundary conditions show better agreement to DSMC than the alternative kinetic theory evaporation boundary conditions for R13. Furthermore, the new evaporation boundary conditions for R13 are implemented in a code for the numerical solution of complex, two-dimensional geometries and compared to NSF solutions. Different flow patterns between R13 and NSF for higher Knudsen numbers are observed.

Copper demand for a country's copper industry has a greater pull effect. China's copper consumption in 2015 has accounted for 50% of the world. The scientific forecast of China's copper demands trend is also an important basis for analyzing its future environmental impact. This paper assumes that China's economy will be developing high, medium and low scenarios, and forecasts economic and social indicators such as total GDP, population and per capita GDP in China from 2016 to 2030. Then, predicted the demand of copper resources in China from 2016 to 2030 by the combination of system dynamics model, vector autoregressive moving average model and inverted U-type empirical model. The results show that: (1) in 2020, 2025 and 2030, China's refined copper demand will be 13 Mt, 15 Mt and 15.5 Mt. (2) China's copper demand growth slowed down significantly from 2016-2030. (3) 2025-2030, China's copper resource demand is stable, into the platform of demand growth, the highest peak value in 2027 will be 15.5 Mt. (4) 2030 years later, China's copper resource demand will enter a slow decline.

An inducer is one of the most important components of centrifugal pumps, whose presence will result in a significant increase in hydraulic performance and pump efficiency. The primary function of the inducer, however, is to delay the destructive phenomenon of cavitation, which has presented a considerable design challenge. Nevertheless, the amount of improvement and increase in inducer performance in both cavitation and non-cavitation modes depends critically on the radial laxity of the blade tip. As part of this study, the performance of the inducer in the cavitation state has been simulated and compared with the experimental data, which are in good agreement. This paper examines the effect of blade tip lagging on cavitation, and the results indicate that this will improve cavitation ­and delay this destructive phenomenon ­but will negatively affect the non-cavitation performance. With the increase in clearance, the range of the return current will increase at the tip of the inducer blade as well.

Recently, most countries suffering from global warming due to the waste gases coming out of the combustion engines used to generate electric power by traditional methods. Therefore, many countries are currently trying to find alternative solutions to this problem by using renewable energies such as solar energy, wind energy, and energy generated from the waves of the ocean and sea to overcome pollution problems. The waves of the seas and oceans have enormous and largely untapped energy, so this work presents an efficient way to use the energy of sea waves to generate electricity. And also, trying to be a way to produce electricity from wave energy in the future. The best suitable places along the shores of Jizan city were inspected to install the buoy system to revenue advantage of the wave supreme height to achieve the height amount of electrical energy. A mathematical model was made to analyze the wave energy and convert it into energy extracted by mechanical force. The mathematical analyzes used the data collected from satellite maps of the numerous severe waves in the Red sea of the Jazan area, and the data published in previous research along that area. It was found that the beach of the Al Shuqaiq area is the greatest for installing the buoy and obtaining the highest electrical energy due to the presence of the highest wave intensity, followed by the beach of the Baysh and Al Morgan area. Also, one of the objectives of this research is to study the design of a device powered by a buoy to use the waves of the Red Sea to generate electric power. The influence of the buoy system design parameters such as the buoy diameter, length of the cylinder, and the length of the connecting rod on the electrical energy generation from wave energy was investigated. The current device is designed with a gearbox to produce continuous power with a single electric generator. A floating mooring device uses the rise and fall of bulges to convert sea wave energy into electrical energy. The device consists of a float, arm, two wheels of different diameters, a gear set, and an electric generator. The effect of the design of several factors on the performance of the device for converting the sea waves energy into electrical energy, including the length of the buoy arm, the wave height by changing the cam diameter, and the conversion ratios between the gear set, to optimize the output power of the wave energy.

To study the coupling mechanism and dynamic response of the end-face movable gear trans-mission system under complex excitation, a specific configuration of end-meshing movable gear reduction mechanism was used to achieve predetermined rigid thrust transmission and mis-matched gear meshing functions, which solved the inherent defects of traditional harmonic gear mechanism thin-wall flexible wheel easy to be damaged by fatigue. Considering the phenomenon of elastic deformation of live teeth accompanied by significant changes in meshing characteristics in the transmission process of end-meshing harmonic reducer, the influence of dynamic meshing parameters, live tooth deformation, time-varying stiffness of tooth meshing, and time-varying backlash on nonlinear dynamic performance was explored, as well as the mechanism of mul-ti-parameter coupling on transmission performance. The nonlinear dynamics model of the end-meshing harmonic reducer is established to solve the chattering prediction problem. Finally, a comprehensive test bed for the transmission system of harmonic reducer with meshing type with adjustable characteristic end was built to verify the correctness of the theoretical model and provide the theoretical and technical basis for exploring the optimal parameter selection of pas-sive vibration suppression problem.

The analysis of high-temperature physicochemical properties of theK₃AlF₆-Na₃AlF₆-AlF₃ molten salt system is crucial for practical production of electrolytic aluminum. In this paper, the physicochemical properties of the K₃AlF₆-Na₃AlF₆-AlF₃ molten salt electrolyte system at 1173 K and standard atmospheric pressure were simulated using molecular dynamics calculations. The effects of Zr、B、C、W、Ti、Hf、Nb、and Ta on the radial distribution function, coordination number, viscosity and conductivity of the electrolyte system are discussed in detail. The simulation results show that the coordination number between Al-F and the conductivity is decreasing as the content of B²- increases. And the coordination number between Al-F and the viscosity is decreasing as the content of Zr⁴⁺ increases. With the addition of impurity elements such as C, W, Ti, Hf, Nb, and Ta, the data show that C has no significant effect on this electrolyte system, while the presence of W will decrease the coordination number between Al-F in the system. The presence of Ti, Hf, Nb, and Ta may intensify the decomposition reaction between one of the Al-F ligands, which increases the coordination number between Al-F and the number of ions in the electrolyte. This work provides theoretical support for the subsequent study of inert anodes containing Zr, B, Ti, Hf, Nb, and Ta.

The adsorption of proteins on polymer is widely used in biosensors. Here, molecular dynamics (MD) simulation was used to study the immobilization of angiotensin converting enzyme II (ACE2) with six initial orientations proposed on polystyrene (PS) at the ambient conditions of pH (4.5, 6, 7, 8, 9.5) and NaCl (0.01, 0.05, 0.1, 0.15, 0.2, 0.25 M). ACE2 immobilization under favorable ambient conditions was characterized by minimum distance (short), settlement time (fast), interaction energy (substantial) and protein configuration (stable). ACE2 orientations proposed in 0.15M NaCl were respectively preferable to (90, 0, 0), (0, 0, 0) and (0, 270, 0), (180, 0, 0) and (0, 90, 0), (90, 0, 0) at pH 4.5, 6, 7, 9.5. ACE2 immobilization was further evaluated at pH 7 by optimizing NaCl concentration. Its proposed orientations of (i) (0, 270, 0), (ii) (0, 0, 0) and (90, 0, 0), and (iii) (0, 90, 0) and (90, 0, 0) were preferable in 0.05, 0.1 and 0.2 M NaCl, respectively. The great significance of cubic-orientation settlement mode provides tangible improvement for the microfabrication of biochips for rapid diagnosis of severe acute respiratory syndrome coronavirus II (SARS-CoV-2).

When dealing with galactic dynamics, or more specifically, with galactic rotation curves, one basic assumption is always taken: the frame of reference relative to which the rotational velocities are given is assumed to be inertial. In other words, fictitious forces are assumed to vanish relative to the observational frame of a given galaxy. It might be interesting, however, to explore the outcomes of dropping that assumption; that is, to search for signatures of non-inertial behavior in the observed data. In this work, we show that the very discrepancy in galaxy rotation curves could be attributed to non-inertial effects. We derive a model for spiral galaxies that takes into account the possible influence of fictitious forces and find that the additional terms in the new model, due to fictitious forces, closely resemble dark halo profiles. Following this result, we apply the new model to a wide sample of galaxies, spanning a large range of luminosities and radii. It turns out that the new model accurately reproduces the structures of the rotation curves and provides very good fittings to the data.

Due to poor ventilation conditions in the tunnel, if gas pipeline leaks, the consequence of the accident will be more serious. Therefore, before the emergency repair, gas in the tunnel needs to be excharged so as not to explode during the repair process. Therefore, it is necessary to study the ventilation of gas in the tunnel. Based on the computational fluid dynamics (CFD) theory and taking the Yanyingshan tunnel section of China-Myanmar pipeline as an example, this paper uses Fluent software to establish the leakage model of the gas pipeline and fan model in the tunnel and analyzes the influence of different fan locations and number of fans on gas concentration. It can be concluded that: (1) the use of press-in method makes it more efficient to discharge gas out of the tunnel. (2) In order to make ventilation efficient, the fan should be arranged in a higher position and needs to be at a distance from the top of the tunnel. (3) Parallel use of two fans has better ventilation effect than single fan.

The gravitational charge should be the energy instead of mass. This modification will lead to some different results and the experiments to test are also presented. Especially, we propose a scheme to achieve the negative energy and gravitational repulsive force in the lab.

The coordination planning between land use and transportation system is an important premise of solving urban transportation problems and realizing land use integration. This study investigates the interactive and feedback relationship between land use and transportation system from the perspective of access management. By integrating the land use and traffic data from Las Vegas Metropolitan area with the system dynamics model, the causal relationship and causal loop diagrams (CLDs) are introduced to analyze the cause-and-effect relationship and quantitative relationship between the factors of the combined system of land use and transportation, and then sub-models partition and system simulation are performed. The systems dynamics model is established by analyzing the relationship between a series of access management techniques, traffic characteristics, and land use features. The results show that system dynamics model can be used as an effective alternative to model the symbiosis relationship of land use and transportation system for urban planning and construction.

In the computational prediction of bridge vibrations due to high-speed train traffic, the most accurate results can be obtained by considering the interaction dynamics between the train, the superstructure, and the supporting structure. To achieve this, a detailed understanding of the coupling properties of all elements is crucial as they significantly influence the calculated vibrations. The studies in this article investigate the influence of different levels of modeling complexity on the computational acceleration results of single-span girder bridges with a ballasted superstructure. A numerical study on an extensive parameter field of single-span girder bridges is conducted to investigate the influence of modeling the bridge structures as coupling beams, i.e., by considering them as two vertically coupled beams representing the track (rails and sleepers) and the supporting structure. The connection between both beams reflects the stiffness and damping properties of the ballasted superstructure and can reproduce its load-distribution capacity. The excitation is applied as either a moving load or a multi-body model of the train, an Austrian Railjet, to evaluate interdependencies of interaction effects between the vehicle and track, and between track and bridge structure. The reference model is a simply-supported Bernoulli-Euler beam excited by moving axle loads. The comparison of acceleration results allows for identifying critical combinations of structural and train parameters for which the implementation of interaction dynamics has a particularly significant impact on the calculated vibrations and quantifying that impact. These findings provide the possibility of formulating structure-dependent recommendations concerning the targeted application of more complex modeling of the structure (coupling beam model) on the one hand and train (multi-body model) on the other.

Fibroin-chitosan composites, especially those containing nanohydroxyapatite, show potential for bone tissue regeneration. The physicochemical properties of these biocomposites depend on the compatibility between their components. In this study, the intermolecular interactions of fibroin and chitosan were analyzed using a molecular dynamics approach. Two types of systems were investigated: one containing acetic acid and the other containing calcium (Ca²⁺) and hydrogen phosphate (HPO₄²⁻) ions mimicking hydroxyapatite conditions. After obtaining the optimal equilibrium structures, the distributions of several types of interactions including hydrogen bonds, ionic contacts, and hydrophobic contacts along with the structural and energetical features were examined. The calculated binding energy values for the fibroin-chitosan complexes confirm their remarkable stability. The high affinity of fibroin for chitosan can be explained by the formation of a dense network of interactions between considered biopolymers. These interactions were found to be primarily hydrogen bonds and ionic contacts involving ALA, ARG, ASN, ASP, GLN, GLU, GLY, LEU, PRO, SER, THR, TYR, and VAL residues. As it was established, the complexation of fibroin with chitosan maintains the β-sheet conformation of the peptide. β-Sheet fragments in fibroin are involved in the formation of significant amount of hydrogen bonds and ionic contacts with chitosan.

Pool boiling is extensively utilized in various fields, such as electronics cooling, battery cooling, nuclear power plants, and chemical processing. In the process of temperature control through boiling, the generation and transportation of bubbles play a crucial role in calculating the heat transfer capacity of the system. Therefore, it holds immense significance to gain a profound understanding of the mechanisms underlying bubble formation, growth, and detachment. Currently, numerous mechanistic explanations and empirical correlations have been proposed to elucidate various parameters of bubbles during the boiling process. These findings were considered valuable references in selecting appropriate boiling media and designing efficient heating surfaces. To comprehensively present the progress of bubble formation and heat transfer during the boiling process, the forces exerted on the bubbles was highlighted in this article. And a meticulous review on bubble force analysis and correlation formulas pertaining to various relevant parameters (e.g. nucleation sites density, bubble growth rate, detachment diameter, detachment frequency) were conducted. This review article also were expected to provide a novel foundation for further exploration of enhanced boiling heat transfer.

Thermal shape memory polymers (SMPs) find increasing applications in biology and medicine due to their promising potential for improved biocompatibility. In such applications, the inevitable penetration of small molecules from the ambient fluid into the polymer influences the shape recovery process and often leads to a reduction of the temperature, at which shape recovery takes place. We show here via molecular dynamics simulations that the size of additive molecules plays a key role for this process. While the effect of concentration on the recovery rate is monotonic in the investigated range, a non-monotonic dependence on the size of additive molecules emerges at temperatures close to the glass transition. This work thus identifies the additives’ size to be a qualitatively novel parameter for switching the recovery process in polymer-based shape memory materials.

In this paper, the current technologies of the regenerative shock absorber systems have been categorized and evaluated. Three drive modes of the regenerative shock absorber systems, namely the direct drive mode, the indirect drive mode and hybrid drive mode are reviewed for their readiness to be implemented. The damping performances of the three different modes are listed and compared. Electrical circuit and control algorithms have also been evaluated to maximize the power output and to deliver the premium ride comfort and handling performance. Different types of parameterized road excitations have been applied to vehicle suspension systems to investigate the performance of the regenerative shock absorbers including that of the nonlinear regenerative shock absorber. The research gaps for comparison of the different drive modes and the nonlinearity analysis of the regenerative shock absorbers are identified and, the corresponding research questions have been proposed for future work.

In this review we summarize the recent development in understanding the role of PIP2 in cellular function and signaling. We first discuss the effect of PIP2 on actin binding proteins addressing the mechanism of the actin cytoskeletal dynamics such as polymerization or depolymerization of the filamentous network or the coupling to membrane to generate forces. Next, we outline the role of PIP2 in membrane dynamics. We summarized how the membrane organization depends upon PIP2 in the presence of ions or transmembrane proteins that are sensitive to membrane curvature. We discuss how clathrin coated pits interact with adaptor proteins during the endocytosis process, which is facilitated by PIP2. Finally, we discuss the role of PIP2 in cell signaling and diseases.

One of the strategies to reduce the contents of the low density lipoproteins (LDLs) in a blood is a hemoperfusion, in which case they are selectively retracted from plasma by an adsorber located outside the patient body. Recently, a photo-controllable smart surface was developed for this purpose characterised by high selectivity and reusability [C.Guo et al., ACS Applied Materials & Interfaces 2022, DOI:10.1021/acsami.2c07193]. We present a mesoscopic model for such a setup involving the azobenzene-containing polymer brush and the model LDL particles. The latter comprise an uniform spherical core covered by a shell of elongated particles representing phospholipids. The system is studied using the coarse-grained molecular dynamics simulation. We examined the dependences of the binding energy on both the length of polymer chains and the grafting density of a brush, and established optimal conditions for the adsorption. These are explained by a competition between the concentration of azobenzenes and phospholipids in the same spatial region, flexibility of polymer chains, and excluded volume effects.

The permeation of dioxin-like pollutants, namely, chlorinated dibenzodioxins and dibenzofurans, through lipid membranes has been simulated using classic molecular dynamics (CMD) combined with the umbrella sampling approach. The most toxic forms of chlorinated dibenzodioxin and dibenzofuran, 2,3,7,8-tetrachloro-p-dibenzodioxin (TCDD) and 2,3,7,8-tetrachlorodibenzofuran (TCDF), and a dioleyl-phosphatidylcholine (DOPC) lipid membrane of 50 Å wide have been chosen for our study. The free energy profile shows the penetration process is largely favoured thermodynamically (DG≈-12 kcal/mol), with a progressively decrease of the free energy until reaching the energy minima at distances of 8Å and 9.5Å from the centre of the membrane for, respectively, TCDD and TCDF. At the centre of the membrane, both molecules display subtle local maxima with free energy differences of 0.5 and 1 kcal/mol with respect to the energy minima for TCDD and TCDF, respectively. Furthermore, the intermolecular interactions between the molecules and the lipid membrane have been characterized at the minima and the local maxima using hybrid quantum mechanics/molecular mechanics energy decomposition analysis (QM/MM-EDA). Total interaction energies of -17.5 and -16.5 kcal/mol have been found at the energy minima for TCDD and TCDF, respectively. In both cases, the dispersion forces govern the molecule-membrane interactions, no significant changes have been found at the local maxima, in agreement with the classical free energy profile. The small differences found in the results obtained for TCDD and TCDF point out the adsorption and diffusion processes through the cell membrane are not related to the different toxicity shown by these pollutants.

The United Nations Organization in its agenda for 2023 approved Health and Well-being as one of the objectives of sustainable development, in this sense, Higher Education Insti-tutions contribute with interventions to improve each human process that generates a healthy life. In Latin America as well as in European countries, studies have been con-ducted to observe psychological well-being at this stage of life. The present study analyzes in a longitudinal way the psychological well-being in three times in people who attend group dynamics programs. The sample is 668 participants, with an age range between 65 and 80 years. The sample has been distributed in two analysis groups controlling for gender and age: group 1 (with group dynamics) and group 2 (without group dynamics). The Spanish adaptation of Ryff's (2006) Psychological Well-being scale α= 0.87 was ap-plied to both groups. The scale is composed of six dimensions of psychological well-being (self-acceptance, positive relationships, autonomy, mastery of the environment, purpose in life and personal growth). Statistical analyses showed significant differences between the two groups: greater psychological well-being and differences in the dimensions of positive relationships and personal growth in group 1. In addition, logistic regression analysis was carried out to test the predictive capacity of psychological well-being by the existence or not of group dynamics, with significant results in several dimensions. With the results obtained, we can conclude that group dynamics have a relevant role in con-tributing to the development of better levels of psychological well-being of the elderly with longitudinal increase in time. Future studies should be aimed at evaluating the value of group dynamics in the way of successful aging.

The main objective of this study is to demonstrate proof of principle that computational fluid dynamics (CFD) modeling is a tool for studying the contribution of covert and overt vascular architecture to the risk of cerebrovascular disease in in sickle cell disease (SCD) as well as uncover one or more mechanism of response to therapy such as chronic red blood cell (cRBC) transfusion. We analyzed baseline (screening), pre-randomization and study exit magnetic resonance angiogram (MRA) images from 10 (5 each from the transfusion and observation arms) pediatric sickle SCD participants in the silent cerebral infarct transfusion (SIT) trial, using CFD modeling. We reconstructed the intracranial portion of the internal carotid artery and branches and extracted the geometry using 3D Slicer. We cut specific potions of the large intracranial artery to include segments of the internal carotid, middle, anterior, and posterior cerebral artery such that the vessel segment analyzed extended from the intracranial beginning of the internal carotid artery up to immediately after (~0.25 inches) the middle cerebral artery branching point. Cut models were imported into Ansys 2021R2/2022R1 and laminar and time-dependent flow simulation was performed. Change in time averaged mean velocity, wall shear stress, and vessel tortuosity were compared between the observation and cRBC arm. We did not observe a correlation between time averaged mean velocity (TAMV) and mean transcranial doppler (TCD) velocity at study entry. There was also no difference in change in time average mean velocity, wall shear stress (WSS), and vessel tortuosity between the observation and cRBC transfusion arms. WSS and TAMV were abnormal for 2 (developed TIA) out of the 3 participants (one participant had SCI) that developed neurovascular outcomes. CFD approaches allows for the evaluation of vascular topology and hemodynamics in SCD using MRA images. In this proof of principle study, we show that CFD could be a useful tool and we intend to carry out future studies with a larger sample to enable more robust conclusions.

At the end of 2018, a high school student asked a question in Zhihu community, claiming that he had proved Goldbach's conjecture. The problem caused an explosive reaction and a large number of users participated in the discussion. And has caused the widespread influence. On January 1, 2019, the questioner issued his "proof". His proof was soon proved wrong. The heated discussion caused by the incident contains a lot of information of social science analysis value. Therefore, we follow up the event in the first time and build a time series dataset for the event. Taking the questioner's "proof" as the dividing line, all the answers, comments, sub comments and user information of writing these texts before and after two days were recorded. This series of temporal information can reflect the dynamic features of the interaction between user opinions, and the impact of exogenous shocks (proof release) on community opinions. The dataset can be used not only for the demonstration of various social network analysis algorithms, but also for a series of natural language processing tasks such as fine-grained sentiment analysis for long texts, as well as multimodal tasks combining natural language processing and social network analysis. This paper introduces the characteristics and structure of the dataset, shows the visualization effect of social network, and uses the dataset train the benchmark model of emotion analysis.

InGaN/GaN samples grown on c-plane sapphire substrate with different In concentrations by metal organic chemical vapor deposition are demonstrated. The subsequent capping GaN layer growth opens a possibility for dislocation reduction due to the lateral strain relaxation in growth geometry. We present the further growth optimization and innovative characterization of InGaN layers overgrown on different structures with varying In concentrations. The photoelectrical and optical properties of the InGaN layers with or without capping GaN layer were investigated by time-resolved picosecond transient grating and temperature dependence photoluminescence. We note a 10-fold increase in carrier lifetime in the InGaN layers when the sample structure changes from PIN to single InGaN layer.

Experimental test facilities are generally characterised using linear transfer functions to relate the wavemaker forcing amplitude to wave elevation at a probe located in the wavetank. Second and third order correction methods are becoming available but are limited to certain ranges of waves in their applicability. Artificial intelligence has been shown to be a suitable tool to find even highly nonlinear functional relationships. This paper reports on a numerical wavetank implemented using the OpenFOAM software package which is characterised using artificial intelligence. The aim of the research is to train neural networks to represent non-linear transfer functions mapping a desired surface-elevation time-trace at a probe to the wavemaker input required to create it. These first results already demonstrate the viability of the approach and the suitability of a single setup to find solutions over a wide range of sea states and wave characteristics.

An axial wind turbine's performance is analyzed numerically using geometrical and operational parameters. The performance of an axial wind turbine is simulated using three-dimensional computational fluid dynamics (CFD) and compared with experimental data, which confirms the validity of the assumptions and method used. A study of the effects of air in clean and sand particlesy conditions is then conducted. A study was also conducted on mechanical power and the impact of the blade tip on the extension of the vortex. It has been found that different sand particle sizes have a significant impact on wind performance numerically from 0.1 to 0.9, and the results indicate that mechanical power and thrust increase as wind velocity increases, but larger sand particle sizes reduce power, thrust, and torque as wind velocity increases. The turbulence intensity causes the pressure distribution on the walls to be irregular, thereby increasing mechanical power and propulsion force. Moreover, as the wind velocity increases, the tip speed ratio increases and the vortexes are subsequently extended.

The conformational state of the activation loop (A-loop) is pivotal for the activation of most protein kinases. Hence, the characterization of the conformational dynamics of the A-loop is important to increase our understanding of the molecular processes related to diseases and to support the discovery of small molecule kinase inhibitors. Here, we carry out a combination of Molecular Dynamics (MD) and Essential Dynamics (ED) analyses to fully map the effects of phosphorylation, ADP, and conformation disrupting (CD) inhibitors (i.e., CD532 and MLN8054) on the dynamics of the A-loop of Aurora-A. MD revealed that the stability of the A-loop in an open conformation is enhanced by single phospho-Thr-288, while paradoxically, the presence of second phosphorylation at Thr-287 decreases such stability and renders the A-loop more fluctuant in time and space. Moreover, we found that such post-translational modification has a significant effect on the direction of A-loop motions. ED analysis suggests that the presence of the phosphate moiety induces the dynamics of Aurora-A to sample two distinct energy minima, instead of a single large minimum, as in unphosphorylated Aurora-A states. This suggests that the conformational distributions of Aurora-A with both single and double phospho-threonine modifications are remarkably different from the unphosphorylated state. In the closed states, binding of CD532 and MLN8054 inhibitors has the effect of increasing the distance of the N- and C-lobes of the kinase domain of Aurora-A, and the angle analysis between those two lobes during MD simulations showed that the N- and C-lobes are kept more open in presence of CD532, compared to MLN8054. As the A-loop is a common feature of Aurora protein kinases, our studies provide a general description of the conformational dynamics of this structure upon phosphorylation and different ligands binding.

: Arboviruses transmitted by Aedes aegypti and Ae. albopictus have been introduced to Florida on many occasions. Infrequently, these introductions lead to sporadic local transmission and, more rarely, sustained local transmission. Both mosquito species are present in Florida, with spatio-temporal variation in population composition. We developed a 2-vector compartmental, deterministic model to investigate factors influencing Chikungunya virus (CHIKV) establishment. The model includes a non-linear, temperature-dependent mosquito mortality function based on minimum mortality in a central temperature region. Latin Hypercube sampling was used to generate parameter sets used to simulate transmission dynamics, following the introduction of one infected human. The analysis was repeated for 3 values of the mortality function central temperature. Mean annual temperature was consistently important in the likelihood of epidemics, and epidemics increased as the central temperature increased. Ae. albopictus recruitment was influential at the lowest central temperature while Ae. aegypti recruitment was influential at higher central temperatures. Our results indicate that the likelihood of CHIKV establishment may vary, but overall Florida is permissive for introductions. Model outcomes were sensitive to the specifics of mosquito mortality. Mosquito biology parameters are variable, and improved understanding of this variation will improve our ability to predict the outcome of introductions.

Fluctuations of protein three-dimensional structures and large-scale conformational transitions are crucial for the biological function of proteins and their complexes. Experimental studies of such phenomena remain very challenging and therefore molecular modeling can be a good alternative or a valuable supporting tool for the investigation of large molecular systems and long-time events. In this mini-review, we present two alternative approaches to the coarse-grained (CG) modeling of dynamic properties of protein systems. We discuss two CG representations of polypeptide chains used for Monte Carlo dynamics simulations of protein local dynamics and conformational transitions and, on other hand, highly simplified structure-based Elastic Network Models of protein flexibility. In contrast to classical Molecular Dynamics the modeling strategies discussed here allow quite accurate modeling of much larger systems and longer time dynamic phenomena. We briefly describe the main features of these models and outline some of their applications, including modeling of near-native structure fluctuations, sampling of large regions of the protein conformational space, or possible support for the structure prediction of large proteins and their complexes.

Turbulent flow is the first and fundamental physical phenomena to evaluate when optimising cost and reducing emissions from an Anode Baking Furnace (ABF). Gas flow patterns, velocity field, pressure drop, shear stress, and turbulent dissipation rate variables are the main operational parameters to be optimised, considering a specific geometry. Computational Fluid Dynamics (CFD) allows simulating physical phenomena using numerical methods with computer resources. In particular, the finite element method is one of the most used methods to solve the flow equations. This method requires a discretisation of the geometry of the ABF, called mesh. Hence, mesh is the main input to the finite element method. A suitable mesh for applying a discretisation method determines whether the problem can be simulated or not. Generating an appropriate mesh remains a challenge to perform accurate simulations. In this work, a comparison between meshes generated using two mesh generation tools is presented. Results of different study cases are included.

In this paper we show how to efficiently implement parallel discrete simulations on Multi-Core and GPU architectures through a real example of application: a cellular automata model of laser dynamics. We describe the techniques employed to build and optimize the implementations using OpenMP and CUDA frameworks. We have evaluated the performance on two different hardware platforms that represent different target market segments: high-end platforms for scientific computing, using an Intel Xeon Platinum 8259CL server with48cores and also an NVIDIA Tesla6V100 GPU, both running on Amazon Web Server (AWS) Cloud, and on a consumer-oriented platform, using an Intel Core i9 9900k CPU and an NVIDIA GeForce GTX 1050 TI GPU. Performance results are compared and analysed in detail. We show that excellent performance and scalability can be obtained in both platforms, and we extract some important issues that imply a performance degradation for them. We also found that current Multi-Core CPUs with large core numbers can bring a performance very near to that of GPUs, even similar in some cases.

The Hollow fiber membrane contactor (HFMC) is an intensification process that combines membrane and absorption technologies, being a technology in expansion, so the dynamics of this process is crucial for an increase in industrial scale-up. In this manuscript a rigorous dynamic model of an HFMC is developed and simulated in Python software for the absorption of post-combustion CO2 using a monoethanolamine (MEA) solution, which was validated using available data in literature in steady state. Dynamic simulations were performed to analyze the model's capability to predict the effect of disturbances on the responses of CO2 capture ratio and loading. Four distinct experiments were conducted by step changes and pulse disturbances at the following variables: gas and liquid volumetric flow rates, CO2 composition in flue gases, and CO2 lean load. The gas phase showed a faster response compared to the liquid. The CO2 capture ratio and rich loading presented opposite responses when disturbing both liquid flows, and for the composition and lean loading variations, the capture ratio presented an inverse response characteristic, while the loading showed a delayed response. As the model proved to be capable of returning to the reference state, it can potentially be used for control studies purposes.

Effect of abrasive shape on the three-body abrasion behaviors of monocrystalline silicon was investigated by Molecular dynamics modelling. The axial ratio of abrasive particle was varied from 1.00, i.e., a complete sphere, to 0.40 to mimic abrasive shape. The particle’s movement turns toward sliding from rolling when the axial ratio becomes less than a critical value 0.46. In the abrasion process, the friction force and normal force showed an approximately sinusoid-like fluctuation for the rolling ellipsoidal particles, while the front cutting of particle caused that friction force increased and became larger than normal force for sliding particles. The phase transformation process was tracked under different particle’ movement patterns. The Si-II and Bct5 phase producing in loading process can partially transform to Si-III/Si-XII phase and backtrack to original crystal silicon under pressure release, which also occurred in the abrasion process. The secondary phase transformation showed difference for particles’ rolling and sliding movements after three-body abrasion. The rolling of particle induced the periodical and inhomogeneous deformation of substrates, while the sliding benefited producing high-quality surface in CMP process. This study aiming to construct more precise model to understand the wear mechanism benefits evaluating the MEMS wear and CMP process of crystal materials.

The β-amyloid cleaving enzyme 1 (BACE1) is regarded as an important target of drug design toward treatment of Alzheimer's disease (AD). In this study, three separate molecular dynamics (MD) simulations and calculations of binding free energies were carried out to comparatively probe identification mechanism of BACE1 on three inhibitors 60W, 954 and 60X. The analyses on MD trajectories indicate that the presence of three inhibitors influences structural stability, flexibility and internal dynamics of BACE1. Binding free energies calculated by using solvated interaction energy (SIE) and molecular mechanics generalized Born surface area (MM-GBSA) methods reveal that the hydrophobic interactions provide decisive forces for inhibitor-BACE1 binding. The calculations of residue-based free energy decomposition suggest that the sidechains of residues L91, D93, S96, V130, Q134, W137, F169 and I179 play key roles in inhibitor-BACE1 binding, which provides a direction for future drug design toward treatment of AD.

This work aims to investigate the effect of Pt concentration on crystal growth mechanism of Platinum-Palladium (Pt-Pd) binary alloy system during the annealing process starting from amorphous phase until some definite temperatures. The calculations have been performed by using molecular dynamic (MD) simulations. Interatomic interactions are described by on Sutton-Chen type Embedded Atom Potential Energy function. In order to understand the main structural properties at the stable and unstable phases, changes in RDF curves versus time have been analysed for different annealing temperatures. Crystalline type bonded pairs have been determined using MD calculations which is required for the computation of Avrami coefficients and for understanding crystal growth mechanism. The results demonstrate that the increase in concentration of Pt during annealing leads to migration of atoms in the crystal lattice points, elimination of dislocations and formation of perfect crystal structure.

Newton's second law, sometimes called the fundamental principle of dynamics, although considered as an irreducible axiom of mechanics based on and validated by experience, might paradoxically depend on the energetic content of our cosmos, through the existence of a dynamical form of aether appearing, in the ancient theory of Fatio de Duillier and Lesage, as an explanation of gravitational forces. In such a framework, Newton's theory is an approximation valid for small velocities as it was the case in special relativity theory.

Motivated by experimental success in designing linear amphiphilic copolymers with selective demicellization capabilities through head-to-tail depolymerization triggered by specific stimuli, we advanced the design paradigm by increasing the number of triggering groups on the copolymer chain through architectural modifications. Utilizing molecular dynamics simulations, we explored the demicellization and cargo release dynamics of linear and miktoarm copolymers, featuring one, two and three hydrophobic block or branches, each capable for head-to-tail depolymerization. Our findings revealed that, under stoichiometric trigger molecule concentrations, miktoarms with three branches exhibited consistently faster depolymerization rates than those with two branches and linear copolymers. Conversely, at constant trigger molecule concentrations, the depolymerization rates of copolymers exhibited more complex behavior influenced by two opposing factors: the excess of trigger molecules, which increased with a decrease in the number of hydrophobic branches or blocks, and simultaneous head-to-tail depolymerization, which intensified with an increasing number of branches. Our study elucidates the intricate interplay between copolymer architecture, trigger molecule concentrations, and depolymerization dynamics, providing valuable insights for the rational design of amphiphilic copolymers with tunable demicellization and cargo release properties.

In an effort to enhance the efficiency and safety of open-pit mines, this study explores the optimization of end slope road parameters and slope structures, specifically focusing on unmanned driving lanes.A significant aspect of the study is the development of a truck trajectory offset model, which considers the different reaction times between automated sensors and human drivers in adapting to environmental changes. To test these concepts, the study uses numerical simulations to confirm the stability of the proposed end slope designs. Using Victory West Mine No. 1 as a case study, the research determines the optimized width for unmanned driving lanes and the maximum angle for safe steepening of end slopes. The findings indicate that the optimized unmanned lane width for NTE240 mining dump trucks is 1743 mm, allowing for a 2-degree increase in the slope angle at the south end slope. This optimization leads to a steep mining stripping volume of 3.2735 million m³ and a coal output of 2.49628 million tons, maintaining a stripping ratio of 1.31m³/t. These results demonstrate that the unmanned driving road width optimization not only ensures slope safety but also significantly boosts the economic benefits of steep mining in open-pit mines.

Most works related to the identification of mathematical nonlinear systems make one suppose that such approaches can always be directly applied to any nonlinear system. This misconception produces a great deal of discouragement when the obtained results are not the expected ones. Thus, the current work hypothesizes that the most information one has about the mathematical structure of the model, the most precise the identification result should be. Therefore, a variant of the Sparse Identification of Nonlinear Dynamics (SINDY) approach is presented to obtain the full mathematical nonlinear model of a high-order system with coupled dynamics, namely, a commercial quadcopter. Furthermore, due to its high sensitivity to inputs, a control system is devised using the identified model to stabilize the quadcopter. This illustrates the effectiveness of the proposed identification method.

Molecular dynamics(MD) simulations are playing an increasingly important role in structure-based drug discovery (SBDD). Here we review the use of MD for proteins in aqueous solvation, organic/aqueous mixed solvents (MDmix) and with small ligands, to the classic SBDD problems: binding mode and binding free energy predictions. The simulation of proteins in their condensed state reveals the solvent structure and preferential interaction sites (hot spots) on the protein surface. This information is largely transferable across all classes of protein ligands (from water to drugs) and can be used very effectively to understand ligand recognition and improve the predictive capability of well-established methods such as molecular docking. MD simulations for protein and drug or drug-like compounds are now being used but are still computationally expensive and can only be applied to specific cases. On the other hand, MDmix simulations can now be used in SBDD and we will describe the latest developments and implementations. We expect to see an increase in the application of these techniques to a plethora of protein targets to identify new drug candidates with the advent of new tools and faster computers.

With the diversified development of media forms, livestreaming e-commerce has become a new sales model. Unlike the traditional sales model, this paper constructs a three-party game model composed of manufacturers, social media influencers and consumers based on the livestreaming model. It explores the equilibrium strategy selections of each participant in the supply chain system. In analyzing the evolutionary game stability of each participant, this paper obtains the equilibrium strategy and stability factors by analysis of the income and expenditure matrix. It uses the simulation model to analyze reasons for strategy selections of different game participants in the livestreaming e-commerce model. The results show that the strategies of manufacturers, social media influencers and consumers' selections have different impacts on their decision-making mechanism, and man-ufacturers are more sensitive to consumers' active participation in livestreaming e-commerce.

The four-way reversing valve is a key component determining whether a heat pump air conditioning system can work. Studying the four-way reversing valve is of great significance for improving the performance of the heat pump air conditioning system. In this paper, the control equation of the sliding block movement process is constructed to simulate the reversing process of the four-way reversing valve in fluid mechanics, and a detailed study is carried out from the sliding block speed, sliding block acceleration, sliding block force, refrigerant leakage, and other aspects. Research has shown that optimizing the structure of the slider can greatly improve the working performance of the four-way reversing valve. Increasing the height of the slider not only improves the pressure relief ability of the four-way reversing valve and prevents the slider from being damaged by high-pressure gas impact, but also reduces the leakage of refrigerant during the directional process; When reducing the length of the slider, although it will improve the pressure relief capacity of the four-way reversing valve, it also increases the amount of refrigerant leakage. The research results provide a theoretical basis and improvement direction for designers to improve the performance of four-way reversing valves.

Unsmooth reversing is one of the most common faults in the four-way reversing valve of air conditioning system, and the airflow in the valve chamber and the shape of the slider of the reversing valve are the main factors in unsmooth reversing. In order to study the influence of the airflow in the chamber and the slider shape on the reversing process, the fluid flow in the valve chamber of a four-way reversing valve is obtained using computational fluid dynamics (CFD) simulations, and the influence of slider structure on slider thrust and resistance of the four-way reversing valve is analyzed. Four-way reversing valves with no-cutting, straight-cutting and arc-cutting sliders are employed to evaluate the reversing characteristics. A pressure measuring device for a four-way reversing valve chamber is designed, which is performed to test the chamber pressure of three different types of sliders of four-way reversing valves. The maximum error of experimental tests and simulations is within 5% of engineering tolerance, so as to verify the reliability of numerical simulation. The results indicate that the chamber pressure and resistance-thrust ratio with the cutting slider are reduced, which is beneficial to the reversing of the four-way reversing valve. Compared with the straight cutting structure, the arc cutting slider four-way reversing valve has a more stable reversing process and better comprehensive reversing performance. The results can provide a theoretical basis for engineers to improve the performance of four-way reversing valves.

Water scarcity is a worldwide concern, considering that water is a limited resource and essential for life. In Brazil, approximately 30% of its population lives in one of the areas that most suffer from lack of water, the semi-arid region covering about 20% of the country's territorial extension. The lack of water, mainly in the northeast of the country, has been a problem for years, where people who live in this territory have been suffering for months from the lack of water and the poor distribution of the resource, which increases the degree of inequality between the regions of the country. The research aims to show the effect of the hydrological cycle on the quality of vegetation and how such processing can end up affecting people's lives and the environment. The municipality of Floresta is within the scope of sub-basin 73, the territory is part of the municipalities benefited from the PISF- São Francisco River integration project works, and the study has raised discussion on the part of the scientific community regarding the effectiveness of the population project and the possible risks caused to the environment. This study carried out a temporal analysis, 1961-2021, to investigate water availability and how the hydrological system can affect society, raising the degree of its importance for a better quality of life and consequently reducing hydric inequalities. It was also evaluated how climate changes can interfere with these environmental changes and the quality of life at the local scale. The hydrological model system used to assess water availability was the Pernambuco Hydrological Response Units SUPer-System. In addition, this system helped to understand the changes that occurred in the studied sub-basin. Yet was used UAV-Unmanned Aerial Vehicles to view relationship between lives in house and environment conditions. The study can obtain accurate and necessary information to understand the dynamics of the hydrological system and the possible changes of social hydric inequalities. Thus, it is concluded that climate change will have different impacts on small scales and lives conditions due to different characteristic's environment. It is very important to carry out studies on a detailed scale to provide better public policies for mitigating the effects of climate change on people's lives.

Climate projections show that southern Indonesia such as West Nusa Tenggara is projected to experience a lower precipitation and higher temperatures. To date, research on climate change impact on Indonesian rice production yield is limited. As climate change is projected to decrease rainfall and to increase temperatures, this paper offers a qualitative analysis using system archetypes to understand the impacts of climate change on rice production. Two system archetypes are identified including Limits to Growth and Success to Successful. Both archetypes explain that rice production is hampered by high minimum temperature as photosynthesis output is decreased by increasing respiration. This paper shows that using a simple tool, system archetypes, we can describe the impacts of climate change on rice production. The outputs of this study such as a causal loop diagram and system archetypes can be a basis to develop a simulation model in understanding the impacts of climate change on main crops.

The buildup dynamics of diode-pumped passively mode-locked solid-state laser is thoroughly explored by the real-time measurement with temporal sampling rate up to 40 GHz. A concise cavity is developed to ensure the transient dynamics purely arising from the gain medium and saturable absorber. Experimental results reveal that the laser output in the buildup process exhibits numerous passively Q-switched pulses followed with a damped relaxation oscillation prior to the stable mode locking. Furthermore, it is confirmed that the laser output has already displayed single clean mode-locked pulses inside the first several Q-switched envelopes before stepping into the stage of relaxation oscillation. The present real-time exploration is expected to provide important information for practical applications with temporal modulation of the pump intensity.

The airborne inoculum of Fusarium circinatum, the fungal pathogen causing Pine Pitch Canker (PPC), is one of the main means of spread of the disease in forest stands and forest nurseries. Since this world-wide known pathogen was introduced in Europe, its biology in this newly infected area still remains scarcely known. To shed more light on this topic, we set an experiment on a naturally PPC infected forest of Monterey pine in Galicia (NW Spain) with the following two goals: (i) to describe the seasonal spore dispersal pattern during one year of regular sampling and (ii) to assess the spatial spore dispersal pattern around the infested plot. Portable rotating arm spore traps were used and complemented with meteorological measurements. The abundance of F. circinatum spores in the samples was evaluated by quantitative PCR (qPCR) with hydrolysis probe. The results showed almost permanent occurrence of the air inoculum throughout the whole year, being detected in 27 of the 30 samplings. No clear temporal trends were observed, but higher air inoculum was favoured by previous lower air temperatures and lower leaf wetness. Conversely, neither rainfall nor air humidity seemed to have any significant importance. The spatial spread of the inoculum was noted to be successful up to a distance of 1000 m in the wind direction, even with winds of just 5 m s^-1. Our study shows that rotating arm spore traps combined with qPCR may be an efficient tool for F. circinatum detection.

The combination of wet-lab experimental data on multi-site combinatorial mutations and machine learning is an innovative method in protein engineering. In this study, we present an improved innovative sequence–activity relationship (innov'SAR) methodology based on novel descriptors and digital signal processing (DSP) to construct a predictive model. In this improved approach, 21 experimental (R)-selective amine transaminases from Aspergillus terreus (AT-ATA) were used as an input to predict higher thermostability than that predicted using the existing data. We successfully improved the determination coefficient (R2) of the model from 0.66 to 0.92. In addition, root mean square deviation (RMSD) and root mean square fluctuation (RMSF) were estimated and conformation analysis based on molecular dynamics simulations was performed to verify the enhanced thermal stability of the screened mutants. The improved innov'SAR algorithm enhanced the predictive accuracy, suggesting a method for modifying the stability of AT-ATA, which may help in directed evolutionary screening and open up new avenues for protein engineering.

Droplet epitaxy allows the efficient fabrication of a plethora of 3D, III-V-based nanostructures on different crystalline orientations. Quantum dots grown on (311)A-oriented surface are obtained with record surface density, with or without a wetting layer. These are appealing features for quantum-dot lasing, thanks to the large density of quantum emitters and a truly 3D lateral confinement. However, the intimate photophysics of this class of nanostructures has not yet been investigated. Here we address the main optical and electronic properties of s-shell excitons in individual quantum dots grown on (311)A substrates with photoluminescence spectroscopy experiments. We show the presence of neutral exciton and biexciton as well as positive and negative charged excitons. We investigate the origins of spectral broadening, identifying them in spectral diffusion at low temperature and phonon-interaction at higher temperature, the presence of fine interactions between electron and hole spin, and a relevant heavy-hole/light-hole mixing. We interpret the level filling with a simple Poissonian model reproducing the power excitation dependence of the s-shell excitons. These results are relevant for the further improvement of this class of quantum emitters and their exploitation as single photon sources for low density samples as well as for efficient lasers for high density samples.

This paper proposes the use of non-uniform extended surfaces installed externally to the tubes of the radiation section of fired heaters, in order to obtain a better heat flux distribution to the coils. To this end, the heat transfer mechanisms present in such equipment were studied through computational fluid dynamics (CFD), using simplified geometries that represent typical sizes of fired heaters. Also, a simplified model for the combustion was considered. Although this model oversimplifies the physics of the problem, it was able to give satisfactory results for the parameters being optimized, considering the main objective of this paper, that is to minimize the non-uniformity of heat flux in the tubes of the radiant section of fired heaters. It was possible to obtain optimized geometric parameters for different types of extended surfaces evaluated, coupling the results of these models with the Particle Swarm optimization method through the use of a response surface technique,. The results indicate a significant improvement in the uniformity of the heat flux distribution to the tubes through the use of the proposed extended surfaces. Thus, this solution reveals to be an interesting alternative to reduce the risks of fluid degradation and coking formation. Future studies must investigate the non-uniformity of the heat flux due to the presence of the flame and consider the interaction between the reactive flow and the participating medium. Nevertheless, this paper presents some results that justify the optimization of such extended surfaces taking into consideration thermal radiation.

The Brazilian livestock is predominantly extensive, with approximately 90% of the production being sustained on pasture, which occupies around 20% of the territory. In the current climate change scenario and where cropland is becoming a limited resource, there is a growing need for a more efficient land use and occupation. It is estimated that more than half of the Brazilian pastures have some level of degradation; however there is still no mapping of the quality of pastures on a national scale. In this study, we mapped and evaluated the spatio-temporal dynamics of pasture quality in Brazil, between 2010 and 2018, considering three classes of degradation: Absent (D0), Intermediate (D1), and Severe (D2). There was no variation in the total area occupied by pastures in the evaluated period, in spite of the accentuated spatial dynamics, with a retraction in the center-south and expansion to the north, over areas of ​​native vegetation. The percentage of non-degraded pastures increased ~12%, due to the recovery of degraded areas and the emergence of new pasture areas as a result of the prevailing spatial dynamics. However, about 44 Mha of the pasture area is currently severely degraded. The dynamics in pasture quality were not homogeneous in property size classes. We observed that in the approximately 2.68 million properties with livestock activity, the proportion with quality gains was twice as low in small properties compared to large ones, and the proportion with losses was three times greater, showing an increase in inequality between properties with more and less resources (large and small, respectively). The areas occupied by pastures in Brazil present an unique opportunity to increase livestock production and make available areas for agriculture, without the need for new deforestation in the coming decades.

As a widely used plastic, the aging and degradation of polyethylene (PE) are inevitable problems, whether the goal is to prolong the life of PE products or address the issue of white pollution. Molecular simulation is a vital scientific tool in elucidating the mechanisms and processes of chemical reactions. To obtain the distribution and evolution process of PE's thermal oxidation products, this work employs the self-consistent charge density functional tight binding (SCC-DFTB) method to perform molecular simulations of the thermal oxidation of PE with different crystallinity and branched structures. We discovered that crystallinity does not affect the thermal oxidation mechanism of PE, but higher crystallinity makes PE more susceptible to cross-linking and carbon chain growth, reducing the degree of PE carbon chain breakage. The branched structure of PE results in differences in pore size between the carbon chains, with larger pores leading to a concentrated distribution of O2 and free radicals subsequently formed. The breakdown of PE is slowed down when free radicals are localized in low-density regions of the carbon chain. The specifics and mechanism of PE's thermal oxidation are clearly revealed in this paper, which is essential for understanding the process in-depth and for the development of anti-aging PE products.

Shut-in after fracturing is generally adopted for wells in shale oil reservoirs, and imbibition occurring in matrix nanopores has been proven as an effective way to improve recovery. In this research, the molecular dynamics (MD) simulation was used to investigate the effects of wettability and pressure on nanopores imbibition during shut-in for a typical shale reservoir Jimsar. The results indicate that the microscopic advancement mechanism of the imbibition front is the competitive adsorption between “interfacial water molecules” at the imbibition front and “adsorbed oil molecules” on the pore wall. The essence of spontaneous imbibition involves the adsorption and aggregation of water molecules onto the hydroxyl groups on the pore wall. The flow characteristics of shale oil suggest that the overall push of the injected water to the oil phase is the main reason for the displacement of adsorbed oil molecules. Thus, shale oil especially the heavy hydrocarbon component in the adsorbed layer tends to slip on the walls. However, the weak slip ability of heavy components on the wall surface is an important reason that restricts the displacement efficiency of shale oil during spontaneous imbibition. The effectiveness of spontaneous imbibition is strongly dependent on the hydrophilicity of the matrix pore's wall. The better hydrophilicity of the matrix pore's wall facilitates higher levels of adsorption and accumulation of water molecules on the pore wall and requires less time for “interfacial water molecules” to compete with adsorbed oil molecules. During the forced imbibition process, the pressure difference acts on both the bulk oil and the boundary adsorption oil, but mainly on the bulk oil, which leads to the occurrence of wetting hysteresis. Meanwhile, shale oil still existing in the pore always maintains a good stratified adsorption structure. Because of the wetting hysteresis phenomenon, as the pressure difference increases, the imbibition effect gradually increases, but the actual capillary pressure gradually decreases and there is a loss in the imbibition velocity relative to the theoretical value. Simultaneously, the decline in hydrophilicity further weakens the synergistic effect of pressure difference because of the more pronounced wetting hysteresis. Thus, selecting an appropriate well pressure enables cost savings and maximizes the utilization of the formation's natural power for EOR.

Explosive boiling is a fast phase transition from an ultra-thin liquid film to vapor under an extremely high heat flux, which typically has been studied using molecular dynamics simulation (MDS) method. The present MDS study investigated the explosive boiling of a liquid argon nanofilm over different solid copper surfaces with different nanowall patterns including parallel and cross nanowalls. For each surface, atomic motion trajectories, number of liquid and vapor argon atoms, heat flux, and, mainly, onset time of explosive boiling were investigated. The simulation results indicated that explosive boiling occurs earlier on parallel and cross nanowall surfaces than on an ideally smooth surface, regardless of the topology and configuration of nanowalls. Moreover, the results revealed that by using the cross nanowall surfaces, the onset time of explosive boiling decreased by 0.7–4% compared to the parallel nanowall surfaces. In addition, it was found that the onset time of explosive boiling strongly depends on the potential energy barrier and the movement space between nanowalls for both parallel and cross nanowall surfaces. Furthermore, simulation findings showed that even though increasing the height of cross nanowalls increases heat flux and temperature of the fluid argon domain, it does not necessarily result in a shorter onset time of explosive boiling. These findings demonstrate the capability of cross nanowall surfaces for explosive boiling, thereby could be utilized in future surface design for thermal management applications.

State transition in the multiple-levels system has the great potential applications in the quantum technology. In this article we employ a deterministic approach in complex space to analyze the dynamics of the 1s-2p electron transition in the hydrogen atom. The electron’s spin motion is embodied in the framework of quantum Hamilton mechanics that allows us to examine the transition dynamics more precisely. The transition is driven by an oscillating electric field in the z-direction. The electron’s transition process can be visualized by monitoring its motion in the complex space. The quantum potential and the total energy proposed in this paper provide new indices to observe the dynamic changes of electrons in the transition process.

Complex human behavior, including interlimb and interpersonal coordination, has been studied from a dynamical system perspective. We review the applications of a dynamical system approach to a sporting activity, which includes continuous, discrete, and switching dynamics. Continuous dynamics identified switching between in- and anti-phase synchronization, controlled by an interpersonal distance of 0.1 m during expert kendo matches, using a relative phase analysis. As discrete dynamics, return map analysis was applied to the time series of movements during kendo matches, and six coordination patterns were classified. Furthermore, state transition probabilities were calculated based on the two states, which clarified the coordination patterns and switching behavior. We introduced switching dynamics with temporal inputs to clarify the simple rules underlying the complex behavior corresponding to switching inputs in a striking action as a non-autonomous system. As a result, we determined that the time evolution of the striking action was characterized as fractal-like movement patterns generated by a simple Cantor set rule with rotation. Finally, we propose a switching hybrid dynamics to understand both court-net sports, as strongly coupled interpersonal competition, and weakly coupled sports, such as martial arts.

In this paper, the nanoscale dissipative mechanisms of a Cu pad in a Ball Grid Array (BGA) packaging structure during isothermal ageing and uniaxial tension were investigated by the molecular dynamics (MD) method and experiments. From the result of the isothermal ageing test, a nonuniform consumption of Cu and large amount of Kirkendall voids were observed at the interface of Cu and Cu₃Sn. To study the effect of pressure and orientation on this phenomenon, MD simulations were conducted on four types of Cu-Cu₃Sn interface structures with different orientations of Cu. By comparing the diffusion coefficients of atoms in those cases, it was found that the tensile stress would inhibit the diffusion of atoms, whereas compressive stress would accelerate it, and this would be more significant under a larger magnitude of stress and temperature. Note that, in the model with the (101) surface Cu at the interface, both Cu and Cu₃Sn have a higher diffusion coefficient compared with the model with (001) surface Cu. Thus, the orientation of Cu will also contribute to the uniform consumption of the pad. Uniaxial tension simulation combined with DXA and CSP analyses on those models also shows the model with (001) surface Cu has a greater mechanical reliability in our simulations and related experiments.

In this work we study the recently developed parametrized partition function formulation and show how we can infer the thermodynamic properties of fermions based on numerical simulation of bosons and distinguishable particles at various temperatures. In particular, we show that in the three dimensional space defined by energy, temperature and the parameter characterizing parametrized partition function, we can map the energies of bosons and distinguishable particles to fermionic energies through constant-energy contours. We apply this idea to both noninteracting and interacting Fermi systems and show it is possible to infer the fermionic energies at all temperatures, thus providing a practical and efficient approach to obtain thermodynamic properties of large fermion systems with numerical simulation. As an example, we present energies for up to 50 noninteracting fermions and up to 20 interacting fermions at all temperatures and show good agreement with the analytical result for noninteracting case.

The most cited evidence for (nonbaryonic) dark matter has been an apparent lack of visible mass to gravitationally support the observed orbital velocity of matters in rotating disk galaxies. Yet measurement of the mass of celestial objects cannot be straightforward, requiring theories derived from the known physical laws along with some empirically established semi-quantitative relationship. The most reliable means for determining the mass distribution in rotating disk galaxies is to solve a force balance equation according to Newton’s laws from measured rotation curves, similar to calculating the Sun’s mass from the Earth’s orbital velocity. Another common method to estimate galactic mass distribution is to convert measured brightness from surface photometry based on empirically established mass-to-light ratio. For convenience, most astronomers commonly assumed a constant mass-to-light ratio for estimation of the so-called “luminous” or “visible” mass, which should not be expected accurate. The mass determined from rotation curve typically exhibit an exponential-like decline with galactrocentric distance, qualitatively consistent with observed surface brightness. This fact scientifically suggests variable mass-to-light ratio of baryonic matter in galaxies without the need for dark matter.

The seasonal and vertical spatial variation characteristics of total N, NH4+—N and NO3-—N in the soil in Picea mongolica forest at different stand ages were studied. The results showed that: (1) NH4+ —N in the soil in Picea mongolica forest was positively correlated with sample date, soil NO3-—N and forest age (P <0.05), and negatively correlated with soil depth (P < 0.01). NO3-—N in the soil in Picea mongolica forest was negatively correlated with sample date and forest age (P < 0.05), and significantly negatively correlated with soil depth (P < 0.01); (2) the content of NH4+—N is greater than that of NO3-—N in each soil layer of Picea mongolica forest at three forest ages. The correlation analysis showed which had an impact on the content of NH4+—N and NO3-—N in the soil in Picea mongolica forest. The content was lower in February and November, and higher in May and August；(3) The contents of total N, NH4+—N and NO3-—N in the soil of Picea mongolica forest in the three forest ages increased with the depth of soil layer (0-50cm), and showed an overall downward trend. The contents of NH4+—N in the soil layer of young forest (0-10cm, 10-20cm and 20-30cm, 30-40cm and 40-50cm) were significantly different (P＜0.05), and the contents of NO3-—N were also significantly different (P＜0.05), while their contents of middle-aged forest and near mature forest increased with the depth of soil layer, There was no significant difference in the content of NH4+—N; (4) The content of NH4+—N in the soil in Picea mongolica forest generally showed the following change trend near mature forest > middle-aged forest > young forest. The change trend of soil NO3-—N content is consistent with that of the content of NH4+—N in the soil in Picea mongolica forest.

This paper deals with the characteristics of the cyclone separator from the Lagrangian perspective to design important dependent variables, develops a neural network model for predicting the separation performance parameter, and compares the predictive performance between the traditional surrogate model and the neural network model. In order to design the important parameters of the cyclone separator based on the particle separation theory, the force acting until the particles are separated was calculated using the Lagrangian-based CFD methodology. As a result, it was proved that the centrifugal force and drag acting on the critical diameter having a separation efficiency of 50% were similar, and the particle separation phenomenon in the cyclone occurred from the critical diameter, and it was set as an important dependent variable. For developing a critical diameter prediction model based on machine learning and multiple regression methods, Unsteady-RANS analyzes according to shape dimensions were performed. The input design variables for predicting the critical diameter were selected as four geometry parameters that affect the turbulent flow inside the cyclone. As a result of comparing the model prediction performances, the ML model showed the 32.5 % of improvement rate of R2 compared to the traditional MLR considering the nonlinear relationship between the cyclone design variable and the critical diameter. The proposed techniques have proven to be fast and practical tools for cyclone design.

High-efficiency rock breaking by hydraulic jetting is the key to radial horizontal drilling technology. In order to improve the drilling efficiency of hydraulic jet rock breaking in radial horizontal wells, based on LS-Dyna display dynamics, a numerical simulation model of single nozzle jet rock breaking is established to analyze the influence of different nozzle parameters on rock breaking effect. Then, the numerical simulation model of the spin multi-nozzle jet bit is established, and the influence of different rotation speeds on the rock breaking effect of the jet bit is analyzed. Finally, the rock breaking drilling characteristics of the spin multi-nozzle jet bit and the conventional multi-nozzle jet bit are compared and analyzed. The results show that when the jet impacts the rock surface, the larger the inclination angle is, the larger the rock breaking width formed by the jet is. The smaller the dip angle, the greater the rock breaking depth. When the inclination angle is greater than 60 °, it is difficult to meet the needs of reaming. The width and depth of the nozzle gradually increase with the increase of the diameter. When the nozzle diameter is greater than 1.3mm, the growth rate of rock breaking depth begins to decrease. The optimum nozzle diameter is 1.3 mm. When v = 50m / s, the damage caused by the jet to the rock surface is very small, because the condition of rock fracture is not reached at this time. This shows that there is a critical value of the water jet impact velocity, and only when the velocity is reached, the rock will break. When the velocity is v = 150m / s, v = 200m / s, v = 250m / s, v = 300m / s, the rock breaks. At the same time, the higher the speed, the higher the degree of rock fracture, the greater the fracture depth, the greater the fracture area, the better the fracture effect. The tangential and radial velocity of the jet increases the shear and tensile failure rate of the sample, and improves the rock breaking efficiency of the jet, which has certain guiding significance for improving the rock breaking drilling efficiency of radial horizontal well drilling.

We propose a Gauss model (GM), a map from time to the bell-shaped Gauss function to model the deaths per day and country, as a quick and simple model to make predictions on the coronavirus epidemic. Justified by the sigmoidal nature of a pandemic, i.e. initial exponential spread to eventual saturation, we apply the GM to existing data, as of April 2, 2020, from 25 countries during first corona pandemic wave and study the model's predictions. We find that logarithmic daily fatalities caused by Covid-19 are well described by a quadratic function in time. By fitting the data to second order polynomials from a statistical chi2-fit with 95\% confidence, we are able to obtain the characteristic parameters of the GM, i.e. a width, peak height and time of peak, for each country separately, with which we extrapolate to future times to make predictions. We provide evidence that this supposedly oversimplifying model might still have predictive power and use it to forecast the further course of the fatalities caused by Covid-19 per country, including peak number of deaths per day, date of peak, and duration within most deaths occur. While our main goal is to present the general idea of the simple modeling process using GMs, we also describe possible estimates for the number of required respiratory machines and the duration left until the number of infected will be significantly reduced.

The molecular dynamics simulation has been used to investigate the structural and transport properties of (Ba_0.5-xSr_x)La_0.5InO_3-δ (x=0, 0.1, 0.2) oxygen-ion conductor. The previous studies reported that the ionic conductivity of Ba-doped LaInO3 decreases because Ba dopant forms narrow oxygen path in the lattice, which could hinder the diffusion of oxygen ion. In this study, we reveal the mechanism to improve the ionic conductivity by Ba and Sr co-doping on La site in LaInO3 perovskite oxide. The results show that the ionic conductivity of (Ba_0.5-xSr_x)La_0.5InO_3-δ increases with increasing numbers of Sr ions, which oxygen diffusion paths including Sr ion have larger critical radius than Ba ions. The RDF calculations showed the heights of peak in composition including Sr ions is lower and broaden, so oxygen ions moved easily into other oxygen sites.

To obtain a given structure in 3D building-up modes, it is necessary to use optimal surfacing modes that determine the amount of heat input in local areas. For this purpose, aluminum bronze was surfacing onto a deformed base and the structure and its characteristics in the surfacing and heat-affected zones were studied by the EBSD method. The heterogeneity of the formation of the structure in each selected zone is established, which indicates the heterogeneity of heat input in local areas of the material in one mode of surfacing. For typical cases of crystallization, a molecular dynamics simulation of crystallization processes with different heat input to the base with characteristics specified based on experimental data was carried out. It has been established that the amount of heat input determines the degree of melting and the inherited defectiveness of growing crystals. The formation of misorientation boundaries and crystallization centers of new grains is determined by the conditions of joint growth of grains with given crystallographic parameters of the computational model. Numerical calculations agree with the experimentally observed results.

Detailed flow distributions in vascular systems are the key to identifying hemodynamic risk factors for the development and progression of vascular diseases. Although computational fluid dynamics (CFD) has been widely used in bioengineering research on hemodynamics predictions, not only are high-fidelity CFD simulations time-consuming and computing-expensive, but also not friendly to clinical applications due to the difficulty of comprehensive numerical calculations. Machine learning (ML) algorithums to estimate the flow field in vascular systems based on the angiographic images of the blood flow using existed diagnostic tools are emerging as a new pathway to facilitate the mapping of hemodynamics. In present work, the dye injection in a water flow was simulated as an analogy of the contrast perfusion in blood flow using CFD. In the simulation, the light passes through the flow field and generates projective images, as an analogy of X-ray imaging. The simulations provide both the ground truth velocity field and the projective images of the flow with dye patterns. A rough velocity field was estimated using the optical flow method (OFM) based on projective images. ML algorithums are then trained using the ground truth CFD data and the OFM velocity estimation as the input. Finally, the interpretable (logistic regression) and deep (neural networks, convolutional neural networks, long short term memory) machine learning models are validated by using parallel in vitro experiments on the same flow setup. The validation results showed that the employed ML model significantly reduced the error rate from 53.5% to 2.5% in average for the v-velocity estimation.

Drug nanosuspensions offer a promising approach to improve bioavailability for poorly soluble drug candidates. Such formulations often necessitate the inclusion of an excipient to stabilize the drug nanoparticles. However, a rationale to choose the correct excipient for a given drug candidate remains unclear. To gain molecular insight into formulation design, this work first utilizes molecular dynamics simulation to computationally investigate drug-excipient interactions for a number of combinations that have been previously studied experimentally. We find that hydrophobic interactions drive excipient adsorption to drug nanoparticles, and that the fraction of polar surface area serves as a predictor for experimental measurements of nanosuspension stability. To test these ideas prospectively, we applied our model to an uncharacterized drug compound, GDC-0810. Our simulations predicted that salt forms of GDC-0810 would lead to more stable nanosuspensions than the neutral form; therefore, we tested the stability of salt GDC-0810 nanosuspensions and found that the salt form readily formed nanosuspensions even without excipient. To avoid computationally expensive simulations in the future, we extended our model by showing that simple, two-dimensional properties of the single drug molecule can be used to rationalize nanosuspension design without simulation. In all, our work demonstrates how computational tools can provide molecular insight into drug-excipient interactions, and aid in rational formulation design.

We develop an action principle to construct the field equations for dissipative/resistive general relativistic two-temperature plasmas, including a neutrally charged component. The total action is a combination of four pieces: an action for a multi-fluid/plasma system with dissipation/resistivity and entrainment; the Maxwell action for the electromagnetic field; the Coulomb action with a minimal coupling of the four-potential to the charged fluxes; and the Einstein-Hilbert action for gravity (with the metric being minimally coupled to the other action pieces). We use a pull-back formalism from spacetime to abstract matter spaces to build unconstrained variations for the neutral, positively, and negatively charged fluid species and for three associated entropy flows. The full suite of field equations is recast in the so-called ``$3 + 1$'' form (suitable for numerical simulations), where spacetime is broken up into a foliation of spacelike hypersurfaces and a prescribed ``flow-of-time''. A previously constructed phenomenological model for the resistivity is updated to include the modified heat flow and the presence of a neutrally charged species. We impose baryon number and charge conservation as well as the Second Law of Thermodynamics in order to constrain the number of free parameters in the resistivity. Finally, we take the Newtonian limit of the ``$3 + 1$'' form of the field equations which can be compared to existing non-relativistic formulations. Applications include main sequence stars, neutron star interiors, accretion disks, and the early universe.

Freezing/thawing indices are important indicators of dynamics of frozen ground on the Qinghai-Tibet Plateau (QTP), especially for areas with limited observations. Based on the numerical outputs of Community Land Surface Model version 4.5 (CLM4.5) from 1961 to 2010, this study compared the spatial and temporal variations between air freezing/thawing indices (2 m above the ground) and ground surface freezing/thawing indices in permafrost and seasonally frozen ground (SFG) across the QTP, after presenting changes of frozen ground distribution in each decade in the context of warming and wetting. The results indicate that an area of 0.60×106 km2 permafrost of QTP degraded to SFG in 1960s-2000s, and the primary shrinkage period occurred in 2000s. The air freezing index (AFI) and ground freezing index (GFI) decreased dramatically at the rates of 71.00℃·d/decade and 34.33℃·d/decade from 1961 to 2010, respectively. In contrast, the air thawing index (ATI) and ground thawing index (GTI) increased strikingly with the values of 48.13℃·d/decade and 40.37℃·d/decade in the past five decades. Permafrost showed more pronounced changes in freezing/thawing indices since 1990s compared to SFG. The changes of thermal regimes in frozen ground shows closely relations to air warming until the late 1990s, especially in 1998 when the QTP has undergone the most progressively warming. However, a sharp increase in the annual precipitation from 1998 seems to take a stronger control role on thermal degradation in frozen ground than the air warming in 2000s. Meanwhile, the prominent vegetation expansion further promotes thermally unstable of frozen ground in this highly wetting period.

Aquaporin-3 (AQP3) is one of the aquaglyceroporins, which is expressed in the basolateral layer of the skin membrane. Studies have reported that human skin squamous cell carcinoma overexpresses AQP3 and inhibition of its function may alleviate skin tumorigenesis. In the present study, we have applied a virtual screening method that encompasses filters for physicochemical properties and molecular docking to select potential hit compounds that bind to the Aquaporin-3 protein. Based on molecular docking results, the top 20 hit compounds were analyzed for stability in the binding pocket using unconstrained molecular dynamics simulations and further evaluated for binding free energy. Furthermore, examined the ligand-unbinding pathway of the inhibitor from its bound form to explore possible routes for inhibitor approach to the ligand-binding site. With a good docking score, stability in the binding pocket, and free energy of binding, these hit compounds can be developed as Aquaporin-3 inhibitors in the near future.

Electrowetting – induced detachment of droplets from solid surfaces is important for numerous applications in the fields of heat transfer and fluid mechanics. The forced oscillations of droplets on solid surfaces and their ability to detach are studied. In this study the process is efficiently simulated by implementing a powerful methodology developed by our team. Our results agree with experiments, showing that optimal detachment, in terms of actuation energy, is achieved when the application of voltage is synchronized with the spreading time of the droplet. Under these conditions, the droplet oscillates with a period close to that of a mirrored Rayleigh droplet. The relation between the droplet’s oscillation period and its physical properties is examined. During voltage-droplet synchronization, the droplet’s ability to detach depends mostly on its contact angle, its viscosity, and the applied voltage. An energy analysis is also conducted, revealing how energy is supplied to the droplet by electrowetting-induced detachment.

In this paper we introduce an ODE system to model the interaction between nature resources and human exploitation in a rich consumeristic society. In this model the rate of change of population depends on wealth per capita, and the rate of consumption has a quadratic growth with respect to population and wealth. As in our previous work \cite{BaCra_1}, we distinguish between renewable and non renewable resources, and we introduce a replenishment term for non renewable resources. We first get some information on the asymptotic behaviour of wealth and population, then we compute all the equilibria of the system and we study their stability when the parameter of exploitation of resources is small. Some numerical simulations confirm the theoretical results, and give suggestions for future research.

The effects of anisotropic properties (wettability and roughness) on microgroove surfaces in pool boiling heat transfer using Novec-7100 as a working fluid were experimentally investigated. The idea of introducing the concept of anisotropic wettability in boiling experiments draws inspiration from biphilic surfaces. The investigation is also motivated by two-phase immersion cooling, which involves phase-change heat transfer by using dielectric liquid as a working fluid. Very few studies have focused on the effects of surfaces with anisotropic properties on boiling performance. Thus, this study aims to examine the pool boiling heat transfer performance on surfaces with microgroove-induced anisotropic properties under the saturation condition. A femtosecond laser–texturing method was employed to create microgroove surfaces with different groove spacing. The results indicated that anisotropic properties affected the heat transfer coefficient and critical heat flux. Relative to the plain surface, microgroove surfaces enhanced the heat transfer performance due to the increased number of bubble nucleation sites and higher bubble detachment frequency. An analysis of bubble dynamics under different surface conditions was conducted with the assistance of high-speed images. The microgroove surface with groove spacing of 100 μm maximally increased the BHTC by 37% compared to the plain surface. Finally, the CHF results derived from experiments were compared with related empirical correlations. Good agreement was achieved between the results and the prediction correlation.

The increasingly common occurrence of antibiotic-resistant bacteria has become an urgent public health issue. There are currently some infections without any effective treatment, which require new therapeutic strategies. An attractive alternative is the design of compounds capable of disrupting bacterial communication known as quorum sensing (QS). In gram-negative bacteria, such communication is regulated by acyl-homoserine lactones (AHLs). QS allows bacteria to proliferate before expressing virulence factors. Our group previously reported that hexyloxy phenyl imidazoline (9) demonstrated 71% inhibitory activity of QS at 100 µM (IC50=90.9 µM) in Chomobacterium violaceum, a gram-negative bacterium. The aim of the present study was to take 9 as a lead compound to design and synthesize three 2-imidazolines (13-15) and three 2-oxazolines (16-18), to be evaluated as quorum sensing inhibitors on C. violaceum CV026. We were looking for compounds with a higher affinity towards the Cvi receptor of this bacterium and the ability to inhibit QS. The binding mode of the test compounds on the Cvi receptor was explored with docking studies and molecular dynamics. It was found that 8-pentyloxyphenyl-2-imidazoline 13 reduced the production of violacein (IC50=56.38 µM) without affecting bacterial growth, suggesting inhibition of quorum sensing. Indeed, compound 13 is apparently one of the best QS inhibitors known to date. Molecular docking revealed the affinity of compound 13 for the orthosteric site of the N-hexanoyl homoserine lactone (C₆-AHL) on the CviR protein. Ten aminoacid residues in the active site of C₆-AHL interacted with 13, and 7 of these are the same as those interacting with AHL. Contrarily, 8-octyloxyphenyl-2-imidazoline 14, 8-decyloxyphenyl-2-imidazoline 15 and 9-decyloxyphenyl-2-oxazoline 18 bound only to an allosteric site and thus did not compete with C₆-AHL for the orthosteric site.

In current transportation hubs, it is a common phenomenon that passengers travel with necessary wheeled-luggage or suitcases. Due to the fact that most luggage occupies certain space within dense passenger crowds with giant mass inertia, its abnormal motions, such as drifting, frequently trig unavoidable local disturbances and turbulence in surrounding pedestrian flows, further increasing security risk. While, current researches primarily concentrate on examining the impact of luggage on crowd evacuation efficiency through scenario-based experiments and analyzing the spatial characteristics of crowd distribution in surveillance videos. There still lacks of the dynamic disturbance mechanisms and crowd stability associated with the combined motion of passengers and luggage. Therefore, this study considers the luggage-laden passenger as a deformable particle, the resulting disturbance on surrounding non-luggage passengers is analyzed and quantified into a nonlinear pressure term. Subsequently, the disturbance propagation model of passenger-luggage is developed by adapting the classical Aw-Rascle traffic flow model with incorporation of the pressure term. Simulation experiments of disturbances caused by luggage drifting and retrograding validate that the proposed propagation model aligns with real-world situations in transportation hubs. Quantitively, when luggage-laden passenger drifted left, the disturbing force of left side crowd can reach a peak of 238N in the passenger density of 3.0 p/m2, and the maximum difference between the left and right sides disturbing force can reach 153N. Furthermore, it is recommended that the proposed model can be applied to crowd flow analysis and risk assessment.

The aeroelastic response of lightweight low-speed aircraft with slender wings under extreme flow turbulence intensity is not well-understood. Experiments on a commmercial unmanned aerial vehicle (UAV) with a 3-meter wingspan and aspect ratio of 13.6 were performed in a large open-return wind tunnel with extreme flow turbulence intensity of T∞≈10%. The structural dynamics of the wing in the bending mode and the flow beneath the wing to capture the effect of aileron deflection was measured using laser displacement sensors and tomographic particle image velocimetry (PIV) techniques. The unsteady lift produced by the wing was also measured using a high-capacity load cell at an angle of attack, α of 2 degrees for three freestream velocities U∞ of 13.4 m/s, 17.9 m/s, and 26.8 m/s representing the UAV’s stall through cruise speed. It was found that high flow turbulence intensity with large integral length scales relative to the wing chord plays a dominant role in the large unsteady lift and wing displacements measured. The power spectral density (PSD) of wing structural vibration shows that flow shedding from the wing and the integral length scales have a significant impact on the overall power inherent in the bending vibration of the wing. Computation of vorticity iso-surfaces in the flow measurement volume surrounding the aileron reveal a striking observation; aileron deflection, δa of 10∘becomes less effective in producing additional spanwise vorticity, which is proportional to circulation and lift, at U∞ of 26.8 m/s because the freestream already has elevated levels of vorticity. A paradigm shift in design is suggested for light aircraft structures with slender wings operating in highly turbulent flow, which is to employ multiple control surfaces in order to respond to this flow and mitigate large bending or torsion displacements and the probability of structural failure.

Several therapeutic monoclonal antibodies are approved by FDA against the PD-1/PD-L1 (programmed death-1/programmed death ligand-1) immune checkpoint, which has been a great success in cancer treatment. However, existing therapeutics, including small molecules inhibitors against PD-L1 checkpoint have certain drawbacks such as high cost and drug resistance that challenge the current available anti-PD-L1 therapy. Thereof, this study presents the screening of 32552 compounds from Natural Product Atlas database against PD-L1, including three steps of structure-based virtual screening and binding free energy refinement for the selection of potent PD-L1 inhibitors. Subsequently, five natural compounds, i.e., Neoenactin B1, Actinofuranone I, Cosmosporin, Ganocapenoid A, and 3-[3-hydroxy-4-(3-methylbut-2-enyl)phenyl]-5-(4-hydroxybenzyl)-4-methyldihydrofuran-2(3H)-one, were collected based on the ADMET (absorption, distribution, metabolism, excretion, and toxicity) profiling and binding free energy (>-70 kcal/mol) for further computational investigation by comparison to co-crystallised ligand, i.e. JQT inhibitor. Based on interaction mapping, explicit 100 ns molecular dynamics simulation, and post binding free energy calculations, all the compounds exhibited intermolecular interactions (hydrogen and hydrophobic) with essential residues and substantial complex stability by comparison to the JQT inhibitor. Collectively, the calculated results advocate the selected natural compounds as the putative potent inhibitors of PD-L1 and, therefore, can be considered for further development of PD-L1 immune checkpoint inhibitor in cancer immunotherapy.