Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is pandemic. Prevention and control strategies require an improved understanding of SARS-CoV-2 dynamics. We did a rapid review of the literature on SARS-CoV-2 viral dynamics with a focus on infective dose. We sought comparisons of SARS-CoV-2 with other respiratory viruses including SARS-CoV-1 and MERS-CoV. We examined laboratory animal, and human studies. The literature on infective dose, transmission, and routes of exposure was limited specially in humans, and varying endpoints were used for measurement of infection. We propose the minimum infective dose of COVID-19 in humans, is higher than 100 particles, possibly slightly lower than the 700 particles estimated for H1N1 influenza. Despite variability in animal studies, there was some evidence that increased dose at exposure correlated with higher viral load clinically, and severer symptoms. Higher viral load measures did not reflect COVID-19 severity. Aerosol transmission seemed to raise the risk of more severe respiratory complications in animals. An accurate quantitative estimate of the infective dose of SARS-CoV-2 in humans is not currently feasible and needs further research. Further work is also required on the relationship between routes of transmission, infective dose, co-infection, and outcomes.

Starting from December 2019, coronavirus disease 2019 (COVID-19) has emerged as a once-in-a-century pandemic with deadly consequences, which urgently calls for new treatments, cures and supporting apparatuses. Remdesivir was reported by World Health Organization (WHO) as the most promising drug currently available for the treatment of COVID-19. Here, we use molecular dynamics simulations and free energy perturbation methods to study the inhibition mechanism of remdesivir to its target SARS-CoV-2 virus RNA-dependent RNA polymerase (RdRp). In the absence of a crystal structure of the SARS-CoV-2 RdRp, we first construct the homology model of this polymerase based on a previously available structure of SARS-CoV NSP12 RdRp (with a sequence identify of 95.8%). We then build the putative binding mode by aligning the remdesivir + RdRp complex to the ATP bound poliovirus RdRp. The putative binding structure is further optimized with molecular dynamics simulations and demonstrated to be stable, indicating a reasonable binding mode for remdesivir. The relative binding free energy of remdesivir is calculated to be -8.28 ± 0.65 kcal/mol, much stronger than the natural substrate ATP (-4.14 ± 0.89 kcal/mol) which is needed for the polymerization. The ~800-fold improvement in the K_d from remdesivir over ATP indicates an effective replacement of APT in blocking of the RdRp binding pocket. Key residues D618, S549 and R555 are found to be the contributors to the binding affinity of remdesivir. These findings demonstrate that remdesivir can potentially act as a SARS-CoV-2 RNA-chain terminator, effectively stopping its RNA reproduction, with key residues also identified for future lead optimization and/or drug resistance studies.

Few data are available on the health impact of tobacco heating products (THPs) at the population level. We used systems dynamics modelling to estimate effects in the established THP market in Japan. We projected effects of THP use in overall mortality up to 2100 and compare those projections against a baseline scenario based on smoking rates pre-THP launch, i.e., smoking only. The model was informed using data from publicly available sources and the literature, including population size, yearly deaths and smoking prevalence with initialisation year (2004) and, births and migration from 2004 to 2065. Transitions between products were estimated from cross-sectional population surveys in Japan. Potential life-years saved with the introduction of THPs was 13 million by 2100 compared with smoking only. In worst-case scenario, population health gains would be seen with THPs risk 10–50% lower risk than smoking. Assuming equal risk for dual use and smoking, THP risk would need to be at least 10% lower than smoking to achieve a population health benefit by 2100. In credible scenarios, substantial population health gains will follow the introduction of THPs in Japan in a relatively short time frame.

Biomolecular structure drives function, and computational capabilities have progressed such that the prediction and computational design of biomolecular structures is increasingly feasible. Because computational biophysics attracts students from many different backgrounds and with different levels of resources, teaching the subject can be challenging. One strategy to teach diverse learners is with interactive multimedia material that promotes self-paced, active learning. We have created a hands-on education strategy with a set of fifteen modules that teach topics in biomolecular structure and design, from fundamentals of conformational sampling and energy evaluation to applications like protein docking, antibody design, and RNA structure prediction. Our modules are based on PyRosetta, a Python library that encapsulates all computational modules and methods in the Rosetta software package. The workshop-style modules are implemented as Jupyter Notebooks that can be executed in the Google Colaboratory, allowing learners access with just a web browser. The digital format of Jupyter Notebooks allows us to embed images, molecular visualization movies, and interactive coding exercises. This multimodal approach may better reach students from different disciplines and experience levels as well as attract more researchers from smaller labs and cognate backgrounds to leverage PyRosetta in their science and engineering research. All materials are freely available at https://github.com/RosettaCommons/PyRosetta.notebooks.

The SARS-CoV-2 was confirmed to cause the regional outbreak of coronavirus disease 2019 (COVID-19) in Wuhan, China. The 3C-like protease (3CL^pro), an essential enzyme for viral replication, is a valid target to compacts SARS-CoV and MERS-CoV. In this research, an integrated library consisting of 1000 compounds from Asinex Focused Covalent (AFCL) library and 16 FDA-approved protease inhibitors were screened against SARS-CoV-2 3CL^pro. Top compounds with significant docking scores and making stable interactions with catalytic dyad residues were obtained. The screening results in identification of compound 621 from AFCL library as well as Paritaprevir and Simeprevir from FDA-approved protease inhibitors as potential inhibitors of SARS-CoV-2 3CL^pro. The mechanism and dynamic stability of binding between the identified compounds and SARS-CoV-2 3CL^pro were characterized using 50 nanoseconds (ns) molecular dynamic (MD) simulation approach. The identified compounds are potential inhibitors worthy of further development as SARS-CoV-2 3CLpro inhibitors/drugs. Importantly, the identified FDA-approved therapeutics could be ready for clinical trials to treat infected patients and help to curb the COVID-19.

The global pandemic caused by the new SARS-COV-2 coronavirus makes it necessary to search for drugs for its control. Within of this research it has been known that the ivermectin drug, a FDA-approved drugs which is formulated as an 80:20 mixture of ivermectin B1a and B1b and used commonly for parasitic infections, has an inhibitory effect on viruses, includes SARS-COV-2 at in vitro level. In the particular case of SARS-COV-2 its mechanism of action remains elusive and controversial. Interestingly, the energy of interaction of ivermectin with any of the proteins the SARS-CoV-2 and the possible structural alterations at the protein level that this drug can cause have not been reported. In this sense, we carried out a bioinformatics study with docking strategies and molecular dynamics to predict the binding and disturbance induced by ivermectin in proteins associated with SARS-CoV-2. We use DockThor and Molegro docking scores. The DockThor server and myPresto software were used to build complexes and dynamics studies, respectively. The results obtained suggested that ivermectin is capable of docking with the 3CL protease and the HR2 domain, and may promote structural changes in these proteins by inducing unfolding/folding. Specifically, ivermectin brings protease to a significantly more deployed conformational state and the HR2 domain to a more compact state compared to the native state. Finally, it is shown that B1a and B1b macrocyclic lactones have a behavior different from to each target protein. These results suggest a possible inhibitory effect against SARS-CoV-2 due to a synergistic role of this drug to spontaneously bind to two important proteins involve in the proliferation of this virus. However, more studies are required on this possible mechanism of action.

The computing power of smartphones has not received considerable attention in the mainstream education system. Most of the education-oriented smartphone applications (apps) are limited to general purpose services like Massive Open Online Courses (MOOCs), language learning, and calculators (performing basic mathematical calculations). Greater potential of smartphones lies in educators and researchers developing their customized apps for learners in highly specific domains. In line with this, we present Fluid Dynamics, a highly accurate Android application for measuring flow properties in compressible flows. This app can determine properties across the stationary normal and oblique shock, moving normal shock and across Prandtl $-$ Meyer expansion fan. This app can also measure isentropic flows, Fanno flows, and Rayleigh flows. The functionality of this app is also extended to calculate properties in the atmosphere by assuming the International Standard Atmosphere (ISA) relations and also flows across the Pitot tube. Such measurements are complicated and time-consuming since the relations are implicit and hence require the use of numerical methods, which give rise to repetitive calculations. The app is an efficient semi-implicit solver for gas dynamics formulations and uses underlying numerical methods for the computations in a graphical user interface (GUI), thereby easing and quickening the learning of concerned users. The app is designed for the Android operating system, the most ubiquitous and capable surveillance platform, and its calculations are based on JAVA based code methodology. In order to check its accuracy, the app's results are validated against the existing data given in the literature.

Cancers are complex, adaptive ecosystems. They remain the leading cause of disease-related death among children in North America. As we approach computational oncology and Deep Learning Healthcare, our mathematical models of cancer dynamics must be revised. Recent findings support the perspective that cancer-microenvironment interactions consist of turbulent flows. As such, cancer pattern formation, protein-folding and metastatic invasion are discussed herein as processes driven by chemical turbulence within the framework of complex systems theory. Current state-of-the-art quantitative approaches used in reconstructing cancer stem cell networks are reviewed. To conclude, cancer stem cells are presented as strange attractors of the Waddington landscape.

Coronavirus disease 2019 (COVID-19) is an emerging infectious disease that was first reported in Wuhan, China and has subsequently spread worldwide. In the absence of any antiviral or immunomodulatory therapies, the disease is spreading at an alarming rate. 5 to 10% of recovered patients in Wuhan test positive again; this suggest that for controlling COVID-19, vaccines may be better option than drugs. A clinical trial to evaluate an anti-COVID-19 vaccine has started recently. However, its efficacy and potency have to be evaluated and validated. As an alternative, we are presenting a first-of-its-kind, designed multi-peptide subunit based epitope vaccine against COVID-19. The vaccine construct comprise an adjuvant, CTL, HTL, and B-cell epitopes joined by linkers. The vaccine is non-toxic, non-allergenic, thermostable and immunogenic with the capability to elicit a humoral and cell-mediated immune response. The findings are validated with high-end computation-based methods. This unique vaccine is made up of 33 highly antigenic epitopes from three proteins that have a prominent role in host receptor recognition, viral entry, and pathogenicity. We advocate this vaccine must be synthesized and tested urgently as public health priority.

A new form of coronavirus called severe acute respiratory disease coronavirus type 2 (SARS-CoV-2) is currently causing a pandemic. A six-month evolutionary history of SARS-CoV-2 is witnessed by characterising the total genome of 821 samples using comparative phylogenomic approaches. Our analyses produced striking inclusive results that may guide scientists/professionals for the past/future of pandemic. Phylogenetic and time estimation analyses suggest the proximate origin of pandemic strain as Guangdong and the origin time as first half of September 2019, not Wuhan and December 2019, respectively. The viral genome experienced a substitution rate similar to other RNA viruses, but it is particularly high in some of the peptides encoding sequences such as leader protein, E gene, orf8, orf10, nsp10, N gene, S gene and M gene and nsp4, while low in nsp11, orf7a, 3C-like proteinase, nsp9, nsp8 and endoRNase. Most strikingly, the divergence rate of amino acid sequences is high proportional to nucleotide divergence. Additionally, specific non-synonymous mutations in nsp3 and nsp6 evolved under positive selection. The exponential growth rate (r), doubling time (Td) and R₀ were estimated to be 47.43 per year, 5.39 days and 2.72, respectively. Comparison of synapomorphies distinguishing the SARS-CoV-2 and the candidate ancestor bat coronavirus indicates that mutation pattern in nsp3 and S gene enabled the new strain to invade human and become a pandemic strain. We arrive at the following main conclusions: (i) six months evolution of viral genome is nearly neutral, (ii) origin of pandemic is not Wuhan and predates formal reports, (iii) although viral population is ongoing an exponential growth, the doubling time is evolving towards shortening, and (iv) divergence rate of total genome is similar to other RNA viruses, but it is prominently high in some genes while low in some others and evolution in these genes should be closely monitored as their protein products intervening to pathogenicity, virulence and immune response.

Terms related to gradients of scalar fields are introduced as scalar products into the formulation of entropy. A Lagrange density is then formulated by adding constraints based on known conservation laws. Applying the Lagrange formalism to the resulting Lagrange density leads to the Poisson equation of gravitation and also includes terms being related to curvature of space. The formalism further leads to terms possibly explaining nonlinear extensions known from modified Newtonian dynamics approaches. The article concludes with a short discussion of the presented methodology and provides an outlook on other phenomena, which might be tackled using this new approach.

This paper numerically investigates a 3D mesh independency study of a straight blade horizontal axis tidal turbine modelled using Computational Fluid Dynamics (CFD). The solution was produced by employing two turbulence models, the standard k-ε model and Shear Stress Transport (SST) in ANSYS CFX. Three parameters were investigated: mesh resolution, turbulence model, and power coefficient in the initial CFD, analysis. It was found that the mesh resolution and the turbulence model affect the power coefficient results. The power coefficients obtained from the standard k-ε model are 15% to 20% lower than the accuracy of the SST model. It can also be demonstrated that the torque coefficient increases with the increasing Tip Speed Ratio (TSR), but drops drastically after TSR = 5 and k-ε model failing to capture the non-linearity in the torque coefficient with the increasing TSR.

The evolution of integrable classical systems leads to conserved quantities and vanishing Poisson brackets. In contrast, such invariants do not exist in the dynamics of non-integrable systems, which include (but are not limited to) deterministic models with long-term chaotic behavior. The object of this review is to briefly survey the mathematical background of nonintegrability and its role in the physics unfolding well-above the Standard Model (SM) scale.

Leigh Van Valen was an American evolutionary biologist who made major contributions to evolutionary theory and is particularly remembered by his groundbreaking paper "A New Evolutionary Law" (1973) where he provided evidence from fossil record data that this law maintains that the probability of extinction within any group remains es­sentially constant through time. In order to explain such unexpected result, Van Valen formulated a very influential idea that he dubbed the "Red Queen hypothesis". It states that the constant decay must be a consequence of evolutionary interactions among connected species within ecological networks. In Van Valen's picture, species do not merely evolve: they also coevolve with other species. As a consequence, when thinking in adaptation to an external environment, the other species must be considered as part (may be a major part) of such external world. Van Valen's law provided the first complex systems theory of coevolutionary dynamics and inspired a whole range of theoretical and experimental developments and scholars from very diverse fields, from economics to physics. In that respect, Leigh Van Valen's contribution percolated far beyond its original formulation. Red Queen arms races are nowadays considered a widespread feature of complex adaptive systems.

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.

Caveolin-1 (Cav-1) is 22 kDa caveolae protein, acts as a scaffold within caveolar membranes. It interacts with alpha subunits of G-protein and thereby regulates their activity. Earlier studies reported elevated or up-regulated levels of caveolin-1 in the serum of prostate cancer patients. Secreted Cav-1 promotes angiogenesis, cell proliferation and anti-apoptotic activities in prostate cancer patients. Cav-1 upregulation is mainly related to prostate cancer metastasis. Keeping above facts in view, the present study was designed to explore Cav-1 as a target for prostate cancer therapy using computational approach. Molecular docking, structural base molecular modelling and molecular dynamics simulations were performed to investigate Cav-1 inhibitors. A predictive model was generated and validated to establish a stable structure. ZINC database of biogenic compounds was used for induced fit docking (IFD) and high throughput virtual screening. The H-bond interactions of the compounds with active site residues of Cav-1 were estimated by IFD and 100 ns long molecular dynamic simulations. The reported compounds showed significant binding and thus can be considered as potent therapeutic inhibitors of Cav-1. This study provides a valuable insight into biochemical interactions of Cav-1 for therapeutic applications and warrants for experimental validation of the predicted ‘active(s)’.

In this paper, the boiling flow inside a nanochannel with 700000 argon particle has been simulated by molecular dynamics (MD) simulation. This approach has been employed to analysis the superheated flow and its heat transfer pattern as well. For all simulations an external thrust force varying from 1 PN to 12 PN is exerted on inlet nanoparticles along the channel to have the forced annular boiling flow. Computations reveal that saturation condition and superheat degree have significant impacts on the liquid-vapor interface. Furthermore, because of the major influence of surface tension throughout a nanochannel, the x-velocity of liquid film and vapor core has not considerable fluctuations and stay smooth. All provided results show the behaviors completely similar to the available outcomes in the literature.

Since crystals are made of periodic structures in space, predicting their three period vectors starting from any values based on the inside interactions is a basic theoretical physics problem. For the general situation where crystals are under constant external stress, we derived dynamical equations of the period vectors in the framework of Newtonian dynamics, for pair potentials recently (doi:/10.1139/cjp-2014-0518). The derived dynamical equations show that the period vectors are driven by the imbalance between the internal and external stresses. This presents a physical process where when the external stress changes, the crystal structure changes accordingly, since the original internal stress can not balance the external stress. The internal stress has both a full kinetic energy term and a full interaction term. It is extended to many-body interactions in this paper. As a result, all conclusions in the pair-potential case also apply for many-body potentials.

Plant life-long organogenesis involves sequential, time and tissue specific expression of developmental genes. This requires activities of Polycomb Group (PcG) and trithorax Group complexes, respectively responsible for repressive Histone 3 trimethylation at lysine 27 (H3K27me3) and activation-related H3K4me3. However, the genome-wide dynamics in histone modifications that occur during developmental processes have remained elusive. Here, we report the distributions of H3K27me3 and H3K4me3 along with transcriptional changes, in a developmental series including Arabidopsis leaf and three stages of flower development. We found that chromatin mark levels are highly dynamic over the time series on nearly half of all Arabidopsis genes. Moreover, during early flower morphogenesis, changes in H3K4me3 prime over changes in H3K27me3 and quantitatively correlate with transcription changes, while H3K27me3 changes occur after prolonged expression changes. Notably, early activation of PcG target genes is dominated by increases in H3K4me3 while H3K27me3 remains present at the locus. Our results reveal H3K4me3 as greater predictor over H3K27me3 for transcription dynamics, unveil unexpected chromatin mechanisms at gene activation and underline the relevance of tissue-specific temporal epigenomics.

The emergence of SARS/MERS drug resistant COVID-19 with high transmission and mortality has recently been declared a deadly pandemic causing economic chaos and significant health problems. Like all coronaviruses, SARS-CoV-2 is a large virus that has many druggable components within its proteome. In this study, we focused on repurposing approved and investigational drugs by identifying potential drugs that are predicted to effectively inhibit critical enzymes within SARS-CoV-2. We shortlisted seven target proteins with enzymatic activities known to be essential at different stages of the virus life cycle. For virtual screening, the energy minimization of a crystal structure or modeled protein was carried out using Protein Preparation Wizard (Schrödinger LLC, 2020-1). Following active site selection based on data mining and COACH predictions, we performed a high-throughput virtual screen of drugs (n=5903) that are already approved by worldwide regulatory bodies including the FDA, using the ZINC database. Screening was performed against viral targets using three sequential docking modes (i.e. HTVS, SP and XP). Our in-silico virtual screening identified ~290 potential drugs based on the criteria of energy, docking parameters, ligand and binding site strain and score. Drugs specific to each target protein were further analyzed for binding free energy perturbation by molecular mechanics (prime MM-GBSA) and pruning the hits to the top 32 candidates. A top lead from each target group was further subjected to molecular dynamics simulation (MDS) using the Desmond module to validate the efficacy of the screening pipeline. All of the simulated hit-target complexes were predicted to strongly interact and with highly stable binding. Thus, we have identified a number of approved and investigational drugs with high likelihood of inhibiting a variety of key SARS-CoV-2 proteins. Follow-up studies will continue to identify inhibitors suitable for combination therapy based on drug-drug synergy to thwart resistance. In addition, the screening hits that we have identified provide excellent probes for understanding the binding properties of the active sites of all seven targets, further enabling us to derive consensus molecules through computer-aided drug design (CADD). While infections are expanding at a rampant pace, it must be recognized that resistance will grow commensurately through either genetic shift and/or genetic drift to all small molecule drugs identified. Vaccines should provide a more permanent solution through prevention, but resistivity is still a possible scenario. Nevertheless, a persistent multi-target drug development program is essential to curb this ongoing pandemic and to keep reemergence in check.

In the promptness of the COVID-19 outbreak, it would be very important to observe and estimate the pattern of diseases to reduce the contagious infection. To study this effect, we developed a COVID-19 analytical epidemic framework that combines with isolation and lockdown effect by incorporating five various groups of individuals. Then we analyze the model by evaluating the equilibrium points and analyzing their stability as well as determining the basic reproduction number. The extensive numerical simulations show the dynamics of a different group of the population over time. Thus, our findings based on the sensitivity analysis and the reproduction number highlight the role of outbreak of the virus that can be useful to avoid a massive collapse in Bangladesh and rest of the world. The outcome of this study concludes that outbreak will be in control which ensures the social and economic stability.

The current outbreak of coronavirus disease (COVID-19) caused by SARS-CoV-2 in Wuhan, China has killed more than 2600 people since December 2019. Currently there is no effective treatment for this epidemic. Drug for anti SARS-CoV-2 are urgently needed. In this study we evaluated two compound libraries containing launched drugs and compounds from 300 kinds of Traditional Chinese Medicine in order to find anti SARS-CoV-2 drugs. Docking and then calculating binding free energy were performed as workflow against four key anti-SARS-CoV-2 drug targets, 3CLpro, PLpro and RdRp from SARSCoV-2, and AAK1 from human as well. As a result, drugs launched with potential for antiviral usage were selected in the hope of providing some knowledge for future drug discovery.

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.

There must exist a reformulation of quantum field theory which does not refer to classical time. We propose a pre-quantum, pre-spacetime theory, which is a matrix-valued Lagrangian dynamics for gravity, Yang-Mills fields, and fermions. The definition of spin in this theory leads us to an eight dimensional octonionic space-time. The algebra of the octonions reveals the standard model; model parameters are determined by roots of the cubic characteristic equation of the exceptional Jordan algebra. We derive the asymptotic low energy value 1/137 of the fine structure constant, and predict the existence of universally interacting spin one Lorentz bosons, which replace the hypothesised graviton. Gravity is not to be quantized, but is an emergent four-dimensional classical phenomenon, precipitated by the spontaneous localisation of highly entangled fermions.

Background: Previous papers have developed a statistical mechanics of neocortical interactions (SMNI) fit to short-term memory and EEG data. Adaptive Simulated Annealing (ASA) has been developed to perform fits to such nonlinear stochastic systems. An N-dimensional path-integral algorithm for quantum systems, qPATHINT, has been developed from classical PATHINT. Both fold short-time propagators (distributions or wave functions) over long times. Previous papers applied qPATHINT to two systems, in neocortical interactions and financial options. Objective: In this paper the quantum path-integral for Calcium ions is used to derive a closed-form analytic solution at arbitrary time that is used to calculate interactions with classical-physics SMNI interactions among scales. Using fits of this SMNI model to EEG data, including these effects, will help determine if this is a reasonable approach. Method: Methods of mathematical-physics for optimization and for path integrals in classical and quantum spaces are used for this project. Studies using supercomputer resources tested various dimensions for their scaling limits. In this paper the quantum path-integral is used to derive a closed-form analytic solution at arbitrary time that is used to calculate interactions with classical-physics SMNI interactions among scales. Results: The mathematical-physics and computer parts of the study are successful, in that there is modest improvement of cost/objective functions used to fit EEG data using these models. Conclusions: This project points to directions for more detailed calculations using more EEG data and qPATHINT at each time slice to propagate quantum calcium waves, synchronized with PATHINT propagation of classical SMNI.

Statistical mechanics investigates how emergent properties of macroscopic systems (such as temperature and pressure) relate to microscopic state fluctuations. The underlying idea is that global statistical descriptors of order and variability can monitor the relevant dynamics of the whole system at hand. Here we test the possibility of extending such an approach to psychotherapy research investigating the possibility of predicting the outcome of psychotherapy on the sole basis of coarse-grained empirical macro-parameters. Four good-outcome and four poor-outcome brief psychotherapies were recorded, and their transcripts coded in terms of standard psychological categories (abstract, positive emotional and negative emotional language pertaining to patient and therapist). Each patient-therapist interaction is considered as a discrete multivariate time series made of subsequent word-blocks of 150-word length, defined in terms of the above categories. Static analysis (Principal Component Analysis) highlighted a substantial difference between good-outcome and poor-outcome cases in terms of mutual correlations among those descriptors. In the former, the patient’s use of abstract language correlated with therapist’s emotional negative language, while in the latter it co-varied with therapist’s emotional positive language, thus showing the different judgment of the therapists regarding the same variable (abstract language) in poor and good outcome cases. On the other hand, the dynamic analysis, based on five coarse-grained descriptors related to variability, degree of order and complexity of the series, demonstrated a relevant case-specific effect, pointing to the possibility of deriving a consistent picture of any single psychotherapeutic process. Overall, the results showed that the systemic approach to psychotherapy (an old tenet of psychology) is mature to shift from a metaphorical to a fully quantitative status.

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.

A tailored model for the assessment of environmental benefits achievable by “light-weighting” in the automotive field is presented. The model is based on the Fuel Reduction Value (FRV) coefficient, which expresses the Fuel Consumption (FC) saving involved by a 100 kg mass reduction. The work is composed of two main sections: simulation and environmental modelling. Simulation modelling performs an in-depth calculation of weight-induced FC whose outcome is the FRV evaluated for a wide range of Diesel Turbocharged (DT) vehicle case studies. Environmental modelling converts fuel saving to impact reduction basing on the FRVs obtained by simulations. Results show that for the considered case studies, FRV is within the range 0.115–0.143 and 0.142–0.388 L/100 km × 100 kg, respectively, for mass reduction only and powertrain adaptation (secondary effects). The implementation of FRVs within the environmental modelling represents the added value of the research and makes the model a valuable tool for application to real case studies of automotive lightweight LCA.

This paper focuses on the problem of finite-time stabilization of homogeneous, non-Lipschitz systems with dilations. A key contribution of this paper is the design of a virtual recursive Hölder, non-Lipschitz state feedback which renders the non-Lipschitz systems in the special case dominated by a lower-triangular nonlinear system, finite-time stable. The proof is based on a recursive design algorithm developed recently, to construct the virtual Hölder continuous, finite-time stabilizer as well as a C¹ positive definite and proper Lyapunov function that guarantees finite-time stability of the non-Lipschitz non-linear systems.

There is a misconception about using the terms photon and electron. When the electron of the outer ring in the silicon atom executes interstate dynamics for only one cycle, it generates force and energy for the unit photon. The unit photon has a shape similar to a Gaussian distribution with turned ends. When a photon of suitable length interacts with the side of the laterally orientated electron of a semisolid or solid atom, it converts into heat. At an approximate angle of 90°, when a photon interacts with the tip of a laterally orientated electron, it divides into bits of energy with shapes similar to integral symbols. Solid or semisolid element atoms can reveal the phenomenon of heat energy if their electrons address the interactions of photons. In the neutral state silicon atom, the center acts as the reference point for electrons executing interstate dynamics, and the north-south tips of the electrons remain along the north-south poles. The energy shapes around the force tracing along the trajectory of electron dynamics. Two forces are exerted on the electron at one time. In interstate dynamics, the electron of the outer ring first reaches the maximum limit point, where the one-bit energy is shaped. In the remaining half cycle, that electron also generates energy of one bit. When there is an uninterrupted supply of heat energy to the silicon atom, electron dynamics generate photons with a shape-like wave. Path-independent but interstate-dependent forces assume the control of an electron. This electron executes dynamics nearly at the speed of light. In dynamics, conservative forces are exerted on position-acquiring electrons. A photon can be unending in length if the electron dynamics remain uninterrupted. The changing aspect of the electron reflects the auxiliary moment of inertia at each turning point. Atoms of suitable elements generate differently shaped photons when executing dynamics for the outer ring electrons. Thus, they can also reveal the phenomenon of photon energy.

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.

It has been more than a year since the first case of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was found. This coronavirus has infected more than 110 million people worldwide by the end of February, 2021, and several virulent as well as more spreadable mutant forms of SARS-CoV-2 have emerged subsequently. In the latter group, three variants B.1.1.7, B.1.351, and P1 lineages, have been reported. Using computer simulation, the present paper investigates the structural differences between the wild type SARS-CoV-2 spike protein and its Asn501Tyr (N501Y) mutant variant. Time-based structural changes between the receptor binding domains of these two species are also examined. The N501Y mutation is common to all the three aforesaid mutant variants.

Food systems will need to undergo considerable transformation. To be better prepared for and resilient to uncertainty and disturbances in the future, resource users and managers need to further develop knowledge about the food and farming system, with its dominating feedback structures and complexities, and to test robust and integrated system-based solutions. This paper investigates how participatory system dynamics modeling can be adapted to groups at the community level with low or no formal educational background. The paper also analyzes the refinement of workshop participants’ mental models as a consequence of a participatory system dynamics intervention. For this purpose, we ran two workshops with small-scale farmers in Zambia. Analysis of workshop data and post-workshop interviews shows that participatory system dynamics is well adaptable to support an audience-specific learning-by-doing approach. The use of pictures, objects and water glasses in combination with the basic aspects of causal loop diagramming makes for a well-balanced toolbox. Participants acquire understanding that is also relevant beyond systems thinking in that is offers a range of practical insights such as a critical evaluation of common food security strategies.

In this paper, we introduce a simple yet powerful and working version of the molecular dynamics code using the Python 3.9 language. The code contents are published in the link given in the appendix 1. The structure and components of the program is given in detail using flowcharts and code snippets. The program consists of major features like velocity verlet integrator, thermostats, COM removal, input and output modules, virial, pressure, and other thermodynamic quantities estimation etc. The author believes that this program will be helpful to graduate students who perform research in molecular dynamics simulations who intend to write their own code instead of the sophisticated open source packages.

Distant Retrograde Orbits in the Earth-Moon system are gaining in popularity as stable `parking' orbits for various conceptual missions. To investigate the stability of potential Distant Retrograde Orbits, simulations were executed, with propagation running over a thirty-year period. Initial conditions for the vehicle state were limited such that the position and velocity vectors were in the Earth-Moon orbital plane, with the velocity oriented such that it would produce retrograde motion about Moon. The resulting trajectories were investigated for stability against the eccentric relative orbits of Earth and Moon in an environment that also included gravitational perturbations from Sun, Jupiter, and Venus, and the effects of radiation pressure. The results appear to indicate that stability is enhanced for certain resonant states within the Earth-Moon system.

The intracellular calcium dynamics in vascular endothelial cells (VECs) in response to wall shear stress (WSS) and/or adenosine triphosphate (ATP) have been commonly regarded as an important factor in regulating VEC function and behavior including proliferation, migration and apoptosis. However, the effects of time-varying ATP signals have been usually neglected in the past investigations in the field of VEC mechanobiology. In order to investigate the combined effects of WSS and dynamic ATP signals on the intracellular calcium dynamic in VECs, a Y-shaped microfluidic device, which can provide the cultured cells on the bottom of its mixing micro-channel with stimuli of WSS signal alone and different combinations of WSS and ATP signals in one single micro-channel, is proposed. Both numerical simulation and experimental studies verify the feasibility of its application. Cellular experimental results also suggest that a combination of WSS and ATP signals rather than a WSS signal alone might play a more significant role in VEC Ca²⁺ signal transduction induced by blood flow.

We have recently proposed a pre-quantum, pre-space-time theory as a matrix-valued Lagrangian dynamics on an octonionic space-time. This pre-theory offers the prospect of unifying the internal symmetries of the standard model with gravity. It can also predict the values of free parameters of the standard model, because these parameters arising in the Lagrangian are related to the algebra of the octonions which define the underlying non-commutative space-time on which the dynamical degrees of freedom evolve. These free parameters are related to the algebra $J_3(\mathbb O)$ [exceptional Jordan algebra] which in turn is related to the three fermion generations. The exceptional Jordan algebra [also known as the Albert algebra] is the finite dimensional algebra of 3x3 Hermitean matrices with octonionic entries. Its automorphism group is the exceptional Lie group $F_4$. These matrices admit a cubic characteristic equation whose eigenvalues are real and depend on the invariant trace, determinant, and an inner product made from the Jordan matrix. Also, there is some evidence in the literature that the groups $F_4$ and $E_6$ could play a role in the unification of the standard model symmetries, including the Lorentz symmetry. The octonion algebra is known to correctly yield the electric charge values (0, 1/3, 2/3, 1) for standard model fermions, via the eigenvalues of a $U(1)$ number operator, identified with $U(1)_{em}$. In the present article, we use the same octonionic representation of the fermions to compute the eigenvalues of the characteristic equation of the Albert algebra, and compare the resulting eigenvalues with the known mass ratios for quarks and leptons. We find that the ratios of the eigenvalues correctly reproduce the [square root of the] known mass ratios for quarks and charged leptons. We also propose a diagrammatic representation of the standard model bosons, Higgs and three fermion generations, in terms of the octonions, exhibiting an $F_4$ and $E_6$ symmetry. In conjunction with the trace dynamics Lagrangian, the Jordan eigenvalues also provide a first principles theoretical derivation of the low energy value of the fine structure constant, yielding the value $1/137.04006$. The Karolyhazy correction to this value gives an exact match with the measured value of the constant, after assuming a specific value for the electro-weak symmetry breaking energy scale.

First, this paper broaches the epistemological status of scientific tenets and approaches: phenomenological (descriptive only), well-founded (solid first principles, conducive to deep understanding), provisional (can be falsified if universal and verified if existential), and imaginary (fictitious entities or processes, conducive to empirically unsupported beliefs). The ΛCDM “concordance model” involves such beliefs: the emanation of the universe out of a non-physical stage, cosmic inflation (invented ad hoc), Λ (fictitious energy), and exotic dark matter. Big Bang cosmology further faces conceptual and pragmatic problems in delimiting what expands from what does not. The problems dissolve after untying inertia from space. The cosmology that emerges appears immediately compatible with the considered observations and the ‘perfect cosmological principle’. Waves and field perturbations that propagate at c expand exponentially with distance (a gravitational effect). The cosmic web of galaxies does not. Potential -Φ varies as H/(cz) instead of 1/r. Inertial forces arise from the gravitational action of the rest of the universe. Due to dilatation, they are reduced disproportionately at low accelerations. A cut-off value a₀ = 0.168 cH is deduced. This explains the successful description of galaxy rotation curves by MoND. A fully elaborated physical theory is still pending. Wider implications are briefly discussed.

Background: Selectively targeting dopamine receptors has been a persistent challenge in the last years for the development of new treatments to combat the large variety of diseases evolving these receptors. Although, several drugs have been successfully brought to market, the subtype-specific binding mode on a molecular basis has not been fully elucidated. Methods: Homology modeling and molecular dynamics were applied to construct robust conformational models of all dopamine receptor subtypes (D1-like and D2-like receptors). Fifteen structurally diverse ligands were docked to these models. Contacts at the binding pocket were fully described in order to reveal new structural findings responsible for DR sub-type specificity. Results: We showed that the number of conformations for a receptor:ligand complex was associated to unspecific interactions > 2.5 Å and hydrophobic contacts, while the decoys binding energy was influenced by specific electrostatic interactions. Known residues such as 3.32Asp, the serine microdomain and the aromatic microdomain were found interacting in a variety of modes (HB, SB, π-stacking). Purposed TM2-TM3-TM7 microdomain was found to form a hydrophobic network involving Orthosteric Binding Pocket (OBP) and Secondary Binding Pocket (SBP). T-stacking interactions revealed as especially relevant for some large ligands such as apomorphine, risperidone or aripiprazole. Conclusions: This in silico approach was successful in showing known receptor-ligand interactions as well as in determining unique combinations of interactions, key for the design of more specific ligands.

Arylamine N-acetyltransferases (NATs) are cytosolic enzymes, highly polymorphic, present in both eukaryotes and prokaryotes. These enzymes play an important role in the detoxification and activation of xenobiotics as well as in the synthesis of endogenous compounds. Specific NATs have been pointed out in the literature as possible therapeutic targets. In particular, the human NAT1, for the treatment of certain cancers, and the NAT from M. tuberculosis (TBNAT), for the treatment of tuberculosis. This paper describes an in silico approach to prospect and select potentially inhibitors of NAT1 and TBNAT from the Traditional Chinese Medicine (TCM) using free available tools. A library with ligands from TCM was previously screened in order to select only compounds with optimal pharmacological properties. The affinity of the selected ligands with respect to NAT enzymes was then evaluated by virtual screening (VS). Subsequently, the complexes with the best ligands were submitted to molecular dynamics (MD) simulations aiming to obtain better quality information on affinity and selectivity. The results for one specific ligand, ZINC14690579, indicated its potential for affinity and selectivity. ZINC14690579 structure may represent the discovery of a new scaffold for future development of NAT inhibitors.

G protein-coupled receptors (GPCRs) play a key role in many cellular signaling mechanisms, and must select among multiple coupling possibilities in a ligand-specific manner in order to carry out a myriad of functions in diverse cellular contexts. Much has been learned about the molecular mechanisms of ligand-GPCR complexes from Molecular Dynamics (MD) simulations. However, to explore ligand-specific differences in the response of a GPCR to diverse ligands, as is required to understand ligand bias and functional selectivity, necessitates creating very large amounts of data from the needed large-scale simulations. This becomes a Big Data problem for the high dimensionality analysis of the accumulated trajectories. Here we describe a new machine learning (ML) approach to the problem that is based on transforming the analysis of GPCR function-related, ligand-specific differences encoded in the MD simulation trajectories into a representation recognizable by state-of-the-art deep learning object recognition technology. We illustrate this method by applying it to recognize the pharmacological classification of ligands bound to the 5-HT2A and D2 subtypes of class A GPCRs from the serotonin and dopamine families. The ML-based approach is shown to perform the classification task with high accuracy, and we identify the molecular determinants of the classifications in the context of GPCR structure and function. This study builds a framework for the efficient computational analysis of MD Big Data collected for the purpose of understanding ligand-specific GPCR activity.

There must exist a reformulation of quantum field theory, even at low energies, which does not depend on classical time. The octonionic theory proposes such a reformulation, leading to a pre-quantum pre-spacetime theory. The ingredients for constructing such a theory, which is also a unification of the standard model with gravitation, are : (i) the pre-quantum theory of trace dynamics – a matrix-valued Lagrangian dynamics, (ii) the spectral action principle of non-commutative geometry, (iii) the number system known as the octonions, for constructing a non-commutative manifold and for defining elementary particles via Clifford algebras, (iv) a Lagrangian with $E_8 \times E_8$ symmetry. The split bioctonions define a sixteen dimensional space (with left-right symmetry) whose geometry (evolving in Connes time) relates to the four known fundamental forces, while predicting two new forces, $SU(3)_{grav}$ and $U(1)_{grav}$. This latter interaction is possibly the theoretical origin of MOND. Coupling constants of the standard model result from left-right symmetry breaking, and their values are theoretically determined by the characteristic equation of the exceptional Jordan algebra of the octonions. The quantum-to-classical transition, precipitated by the entanglement of a critical number of fermions, is responsible for the emergence of classical spacetime, and also for the familiar formulation of quantum theory on a spacetime background.

An important stage in SARS-CoV-2 life cycle is the fusion of spike(S) protein with the ACE2 host-cell receptor. Therefore, to explore conserved features in S protein dynamics and to identify potentially novel regions for drugging, we measured variability derived from 791 viral genomes and studied its properties by MD simulation. The findings indicated that S2 subunit (HR1, CH, and CD domains) showed low variability, low fluctuations in MD, and displayed a trimer cavity. By contrast, the RBD domain, which is typically targeted in drug discovery programmes, exhibits more sequence variability and flexibility. Interpretations from MD suggest that the monomer is in constant motion showing transitions up-to-down state, and the trimer cavity may function as a 'bouncing spring' that may facilitates S protein interactions with ACE2. Feasibility of trimer cavity for potential drug target was examined by SBVS screening. Several hits that have already been validated or suggested to inhibit the SARS-CoV-2 virus in cell systems were identified; in particular, the data suggest an action mechanism for such molecules including Chitosan and macrolide types. These findings identify a novel binding-site formed by the S protein, that might assist in future drug discovery programmes aimed at targeting the CoV family of viruses.

After a local pollution control plant significantly reduced phosphorus loading into a phytoplankton-dominated tributary of the Potomac River in the early 1980’s, water quality and biological communities were monitored bi-weekly from April-September. After a 10-year time-lag, submerged aquatic vegetation (SAV), once abundant in this freshwater tidal embayment, returned to the area in 1993. After additional reductions in nitrogen load starting in 2000, the system completely switched to a SAV-dominated state in 2005. Fish abundance didn’t change during these distinct phase changes, but the fish community structure did. Increases in SAV provided refuge and additional spawning substrate for species with adhesive eggs such as Banded Killifish (Fundulus diaphanus), which is now the most abundant species in the embayment. Other changes seen were a decrease in the relative contribution of open water dwelling species such as White Perch (Morone americana), and an increase of visual predators such as Largemouth Bass (Micropterus salmoides). The 30-year record of data from this Potomac River tributary has revealed many important long-term trends that validate the effectiveness of initiatives to improve water treatment, and will aid in the continued management of the watershed and point source inputs.

Assessment of the trends of land cover and vegetation dynamics (VD) using remote sensing (RS) and indicators such as anthropogenic activities and the socio-demographic information is essential in order to make proper planning for sustainable management. This paper attempts to evaluate land cover change (LCC) and VD in Kebbi State, Nigeria using historical Landsat data from 1986-2016 by means of remote sensing. The Driver-Pressure-State-Impact-Response (DPSIR) framework was later employed using both primary and secondary data for a better understanding of the drivers, the state of the environmental condition, the causes as well as the impact of the change. The images were classified into five thematic land cover classes as Dense Vegetation, shrubs/built area, farmland, bare/grassland and water body by means of Maximum likelihood supervised classification technique in accordance with Anderson classification scheme level 1, with acceptable accuracy. Pre-classification and post-classification change detection (CD) methodologies were executed using Normalized difference vegetation index (NDVI) and Image differencing respectively. The study illustrates a steady decline in dense vegetation and shrubs/build areas while farmland and bare/grassland increases, however, water bodies remain unchanged. The DPSIR pin-point that the major drivers of change in the study area have been the pressing need for farming land as the population grows and socioeconomic demands including fuelwood consumption and endemic poverty. Expansion of Farming land, fuelwood consumption and the need for construction materials are identified as the main key elements exerting pressure for the change. The state of the condition indicates a steady decline in dense vegetation and shrubs areas while farmland and bare/grassland are increasing significantly. The impacts include land degradation, the decline in the provision ecosystem goods and services, biodiversity loss through loss of habitats. The study, however, noted that many international and national policies in response to land degradation are channelled toward land restoration and remediating of the environment, through afforestation programs and improving the livelihood of the rural people through providing alternative income sources since they depend heavily on land for sustenance. However, the state governments, communities and individual commonly organized annual tree planting campaign with the main purpose of environmental protection.

As paradigm of complex behavior, Self-organized Criticality (SOC) reflects the ability of nonequilibrium dynamical systems to self-adjust into metastable states that are scale independent. The goal of this report is to tentatively show that the hierarchy of Standard Model masses and mixing angles follows from the universal scaling attributes of SOC.

Waddington’s epigenetic landscape, a famous metaphor in developmental biology, depicts how a stem cell progresses from an undifferentiated phenotype to a differentiated one. The concept of “landscape” in the context of dynamical system theory represents a high-dimensional cell state space, in which each cell phenotype is considered as an “attractor” that is determined by interactions among multiple variables (molecular players), and is buffered against environmental fluctuations. Further, biological noise is thought to play an important role during these cell-fate decisions and in fact controls transitions between different phenotypes. Here, we discuss these phenotypic transitions in cancer from a dynamical systems perspective and invoke the concept of “cancer attractors” – hidden stable states of the underlying regulatory network that are not occupied by normal cells. Using epithelial-to-mesenchymal transition (EMT), cancer stem-like properties, metabolic reprogramming and the emergence of drug/hormone resistance as examples, we illustrate how phenotypic plasticity in cancer cells enables them to acquire hybrid phenotypes (such as hybrid epithelial/mesenchymal and hybrid metabolic phenotypes) that tend to be more aggressive and notoriously resilient to drug/hormone treatment. Furthermore, we highlight multiple factors that may give rise to phenotypic plasticity in cancer cells, such as (a) multi-stability or oscillatory behaviors governed by underlying regulatory networks involved in cell-fate decisions in cancer cells, and (b) network rewiring due to conformational dynamics of intrinsically disordered proteins (IDPs) that are highly enriched in cancer cells. We conclude by discussing why a therapeutic approach that promotes ‘recanalization’, i.e. the exit from “cancer attractors” and re-entry into “normal attractors”, is more likely to succeed rather than a conventional approach that targets individual molecules/pathways.

G protein-coupled receptors (GPCRs) are critical drug targets. GPCRs convey signals from the extracellular to the intracellular environment through G proteins. There is evidence that some ligands that bind to the GPCRs activate different downstream signaling pathways. G protein activation or -arrestin biased signaling involves ligands binding to receptors and stabilizing conformations that trigger a specific pathway. Molecular dynamics (MD) simulations are especially valuable for obtaining detailed mechanistic information, including identification of allosteric sites and understanding modulators' interactions between receptors and ligands. Here, we highlight recent simulation studies and methods used to study biased G protein-coupled receptor signaling and their conformational dynamics. We also highlight applications of MD simulations to drug discovery.

In this study, we present case studies to explore the dynamics of implementation and maintenance of health interventions. We analyze how specific interventions are built and eroded, how the building and erosion mechanisms are interconnected, and why we can see significantly different erosion rates across otherwise similar organizations. We use multiple comparative obesity prevention case studies to provide empirical information on the mechanisms of interest, and use qualitative systems modeling to integrate our evolving understanding into an internally consistent and transparent theory of the phenomenon. Our preliminary results identify reinforcing feedback mechanisms, including design of organizational processes, motivation of stakeholders, and communication among stakeholders, which influence implementation and maintenance of intervention components. Over time, these feedback mechanisms may drive a wedge between otherwise similar organizations, leading to distinct configurations of implementation and maintenance processes.

The aim of this paper is to attempt to understand why the popular academic and policy field of promoting, studying and evangelising “entrepreneurship” should have been associated with great success but, in the past twenty years or more in many advanced economies, so much failure. From the US to lesser and developing countries, emerging economies and the European Union, entrepreneurship, especially in regard to start-ups and particularly high-tech start-ups have been in constant more or less recent decline. This is seldom registered in the mainstream literature where a positive and benign profile is generally presented. The paper examines this phenomenon, ties it partly with the “productivity paradox” and seeks tentative hypotheses in relation to the apparent illusions if not delusions regarding “entrepreneurship”.

Starting with a hypothetical geostrophic zonal current in an unbounded ocean, the investigation points out the response of this simple system to a thermal forcing, applied to the free surface and consistent with the maintenance of the geostrophic balance. The main result is the formation of a meridional component of the current, according to the Sverdrup relation, such that the full velocity vector rotates clockwise for heating and anticlockwise for cooling to adjust eventually in the initial zonal direction for large depths.

In this paper, a two-phase model of air shock wave induced by rock-fall was described. The model was made up of the uniform motion phase (velocity was close to 0 m·s^-1) and the acceleration movement phase. The uniform motion phase was determined by experience, meanwhile the acceleration movement phase was derived by the theoretical analysis. A series of experiments were performed to verify the two-phase model and obtained the law of the uniform motion phase. The acceleration movement phase was taking a larger portion when height of rock-fall was higher with the observations. Experimental results of different falling heights showed good agreements with theoretical analysis values. Computational fluid dynamics (CFD) numerical simulation had been carried out to study the variation velocity with different falling height. As a result of this, the two-phase model could accurately and convenient estimating the velocity of air shock wave induced by rock-fall. The two-phase model could provide a reference and basis for estimating the air shock waves' velocity and designing the protective measures.

Background: Since March 2020, Spain is severely hit by the ongoing pandemic of coronavirus disease 19 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Understanding and disrupting the early transmission dynamics of the infection is crucial for impeding sustained transmission. Methods: We recorded all COVID-19 cases and traced their contacts in an isolated rural community. We also sampled 10 households, 6 public service sites and the wastewater from the village sewage for environmental SARS-CoV-2 RNA. Results: The first village patient diagnosed with COVID-19-compatible symptoms occurred on March 3, 2020, twelve days before lockdown. A peak of 39 cases occurred on March 30. By May 15, the accumulated number of symptomatic cases was 53 (6% of the population), of which only 22 (41%) had been tested and confirmed by RT-PCR as SARS-CoV-2 infected, including 16 hospitalized patients. Contacts (n=144) were six times more likely to develop symptoms. Environmental sampling detected SARS-CoV-2 RNA in two households with known active cases and in two public service sites: the petrol station and the pharmacy. Samples from other sites and the wastewater tested negative. Conclusions: The low proportion of patients tested by RT-PCR calls for urgent changes in disease management. We propose that early testing of all cases and their close contacts would reduce infection spread, reducing the disease burden and fatalities. In a context of restricted testing, environmental RNA surveillance might prove useful for early warning and to identify high-risk settings enabling a targeted resource deployment.

The concept of smart morphing blades, which can control themselves to reduce or eliminate the need for active control systems, is a highly attractive solution in blade technology. In this paper an innovative passive control system based on Shape Memory Alloys (SMAs) is proposed. On the basis of previous thermal and shape characterization of a single morphing blade for a heavy-duty automotive cooling axial fan, this study deals with the numerical analysis of the aerodynamic loads acting on the fan. By coupling CFD and FEM approaches it is possible to analyze the actual blade shape resulting from both the aerodynamic and centrifugal loads. The numerical results indicate that the polymeric blade structure ensures proper resistance and enables shape variation due to the action of the SMA strips.

Self-organized criticality (SOC) is a universal mechanism for self-sustained critical behavior in large-scale systems evolving outside equilibrium. Our report explores a tentative link between SOC and Lagrangian field theory, with the long-term goal of bridging the gap between complex dynamics and the non-perturbative behavior of quantum fields.

An evolutionary game is introduced which considers game-theoretic strategies in the context of non-linear population matrix models. This game considers the states and actions of the organisms of the evolving population, and a notion of dynamic equilibrium between strategies is described. The game’s formalism is expounded and a proof about equilibrium is given; specifically that any stable equilibrium can be described by proportions of pure strategies; particularly when population matrices are not defective.

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.

Tree-killing bark beetles are the most economically important insects in conifer forests worldwide. Yet despite >200 years of research, the drivers of population eruptions or crashes are still not fully understood, precluding reliable predictions of the effects of global change on beetle population dynamics and impacts on ecosystems and humans. We critically analyze potential biotic and abiotic drivers of population dynamics of the European spruce bark beetle (Ips typographus) and present a novel ecological framework that integrates the multiple drivers governing this bark beetle system. We call for large-scale collaborative research efforts to improve our understanding of the population dynamics of this important pest; an approach that might serve as a blueprint for other eruptive forest insects.

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.

Fractional-time Schrödinger equation (FTSE) describes the evolution of quantum processes endowed with memory effects. FTSE manifestly breaks all consistency requirements of quantum field theory (unitarity, locality and compliance with the clustering theorem), unless the order of fractional differentiation and integration ( ) falls close to the standard index . Working in the context of the minimal fractal manifold (where , ), we confirm here that FTSE approximates the attributes of gravitational metric and provides an unforeseen generation mechanism for massive fields.

Nuclear export of viral mRNAs, is an essential step in the HIV replication cycle. This role is played by a small regulatory protein of HIV-1 called Rev.The N-terminal region of Rev contains an arginine-rich sequence. The arginine-rich motif (ARM) is located between amino acids 38-50 and forms an alpha-helical secondary structure. Expression of the structural proteins of human immunodeficiency virus type 1 requires the direct interaction of multiple copies of the viral Rev protein with its highly structured RNA target sequence, the Rev Response element (RRE). The major viral proteins are not produced if this transport of RNA is stopped. Therefore, knowledge of Rev structure is essential for understanding of its cooperative binding to the RRE, for understanding the mechanism of HIV infection and for the development of antiviral drugs that interfere with Rev’s essential functions and for acknowledgment of good candidate drugs for treatment of AIDS. To understand how REV interact with RRE element of HIV-RNA and its formation of oligomeric complex it is better to characterize the domain wise structure of REV with regard in function of each domain. Due to lack of structural data on Rev no single compound is reported as inhibitor of REV expect antiviral drugs. Identification of a high-affinity RNA-binding site for the human immunodeficiency virus type 1 Rev Protein is much more important. The ARM is a highly specific sequence which allows for the multimerization of Rev and also binding of REV with RNA. Here we are first time exploring the structural characteristics of REV protein both in free form and in complex with RNA at domain function level especially explore the role of ARM motif in REV HIV-1 protein as RNA binding sites by molecular dynamics (MD) simulation and homology modeling studies. Results indicate that the arginine-rich motif (ARM) is crucial in stability of this complex. The residues ARG38, 39, 41, 43, 44, 48, 50, and ASN40 are most interacting with nucleobases of RRE in Crystal structure of Rev and Rev-response-element RNA complex. Our study plays a major role in elaboration of binding of RNA with REV and pave the way for further investigation for therapeutically agent for HIV.

Contemporary climate change is exposing plant populations to novel combinations of temperatures, drought stress, [CO₂] and other abiotic and biotic conditions. These changes are rapidly disrupting the evolutionary dynamics of plants. Despite the multifactorial nature of climate change, most studies typically manipulate only one climatic factor. In this opinion piece, we seek to explore how climate change factors interact with each other and with biotic pressures to alter evolutionary processes. We first explore the ramifications of climate change for key life history stages (germination, growth and reproduction). We then examine how mating system variation influences population persistence under rapid environmental change and propose that mixed mating could be advantageous in future climates. Furthermore, we discuss how spatial and temporal mismatches between plants and their mutualists and antagonists could promote or constrain adaptive responses to climate change. For example, plant-virus interactions vary from highly pathogenic to mildly facilitative, and are partly mediated by temperature, moisture availability and [CO₂]. Will host plants exposed to novel, stressful abiotic conditions be more susceptible to viral pathogens? Finally, we propose novel experimental approaches that could illuminate how plants will cope with unprecedented global change, such as resurrection studies combined with experimental evolution, genomics or epigenetics.

As the interface between human and machine becomes blurred, hydrogel incorporated electronics and devices have emerged to be a new class of flexible/stretchable electronic and ionic devices due to their extraordinary properties, such as soft, mechanically robust and biocompatible. However, heat dissipation in these devices could be a critical issue and remains unexplored. Here, we report the experimental measurements and equilibrium molecular dynamics simulations of thermal conduction in polyacrylamide (PAAm) hydrogels. The thermal conductivity of PAAm hydrogels can be modulated by both the crosslinking density and water content in hydrogels. The crosslinking density dependent thermal conductivity in hydrogels varies from 0.33 to 0.51 Wm^-1K^-1, giving a 54% enhancement. We attribute the crosslinking effect to the competition between the increased conduction pathways and the enhanced phonon scattering effect. Moreover，water content can act as filler in polymers which lead to nearly 40% enhancement in thermal conductivity in PAAm hydrogels with water content vary from 23 to 88 wt%. Furthermore，we find the thermal conductivity of PAAm hydrogel is insensitive to temperature in the range of 25 ^oC – 40 ^oC. Our study offers fundamental understanding of thermal transport in soft materials and provides design guidance for hydrogel-based devices.

Many studies discuss carbon-based materials because of the versatility of its element. They include different opinions for scientific problems and discuss fairly at convincing and compelling levels within the scope and application. A gas-state carbon atom converts into various states depending on its conditions of processing. The electron transfer mechanism in the gas-state carbon atom is responsible to convert it into various states, namely, graphite, nanotube, fullerene, diamond, lonsdaleite and graphene. The shape of ‘energy trajectory’ enables transferring electrons from the left- and right-sides of an atom is like a parabola. That ‘energy trajectory’ is linked to states (filled state and suitable nearby unfilled state) where force-exertion along the poles of transferring electrons is remained balance. So, the mechanism of originating different states of a gas-state carbon atom is under the involvement of energy first. This is not the case for atoms executing confined inter-state electron-dynamics as the force is involved first. Graphite-, nanotube- and fullerene-state atoms ‘partially evolve partially develop’ (form) their structures. These possess one-dimensional, two-dimensional and four-dimensional ordering of atoms, respectively. Their structural formation also comprises ‘energy curve’ having a shape-like parabola. Transferring suitable filled state electron to suitable nearby unfilled state is under a balance force exerting along the poles. The graphite structure under only attained-dynamics of atoms can also be formed but in two-dimension. Here, binding energy between graphite-state carbon atoms is for a small difference of exerting forces along their opposite poles. Structural formation in diamond, lonsdaleite and graphene atoms involve energy to gain required infinitesimal displacements of electrons through which they maintain orientationally-controlled exerting forces along dedicated poles. In this study, the growth of diamond is found to be south to east-west (ground) where atoms bound ground to south. Thus, diamond atoms merge for a tetra-electron ground to south topological structure. Lonsdaleite atoms merge for a bi-electron ground to just-south topological structure. The growth of graphene is found to be north to ground where atoms bound ground to north. Thus, graphene atoms merge for a tetra-electron ground to north topological structure. Glassy carbon exhibits layered-topological structure where, tri-layers of gas-, graphite- and lonsdaleite-state atoms successively bind in repetitive order. Nanoscale hardness is also sketched based on different force-energy behaviors of different state carbon atoms. Here, structure evolution in each carbon state atom explores its own science.

This article attempts to formulate a mathematical model for a potential explanation regarding the unavoidable impact of a rigid body's peculiar shape on the seamless flow over it. The solid body completely immersed in a Newtonian fluid and respectively has a relative open circuit flow on it will typically experience various observable phenomena like flow separation, flow transition, down-wash, stalling at the higher angle of attack, stalling velocity and how cambered airfoil can typically generate lift at a zero incidence angle. This article respectively represents an understanding of the laminar flow over a rigid body's external surface with due respect to its distinctive shape and size. This working paper formulates a more realistic and simplified mathematical model for open circuit laminar flow over a body, based on the historical data of aerodynamics and theoretical mechanics. This is intended to properly estimate forces on the continuous surface of the body in a laminar flow, to properly explain, understand and predict mentioned phenomena. Most of all the mechanism of streamline formation and its deformation with due regards to flow, shape and size of the body in an open-circuit laminar are formulated mathematically to enhance better design theory which can reduce experimentation while designing a streamlined body.

In the present work we propose a dynamical mathematical model of the lung cells inflammation process in response to SARS-CoV-2 infection. In this scenario the main protease Nsp5 enhances the inflammatory process, increasing the levels of NF kB, IL-6, Cox2, and PGE2 with respect to a reference state without the virus. In presence of the virus the translation rates of NF kB and IkB arise to a high constant value, and when the translation rate of IL-6 also increases above the threshold value of 7 pg mL^-1 s^-1 the model predicts a persistent over stimulated immune state with high levels of the cytokine IL-6. Our model shows how such over stimulated immune state becomes autonomous of the signals from other immune cells such as macrophages and lymphocytes, and does not shut down by itself. We also show that in the context of the dynamical model presented here, Dexamethasone or Nimesulide have little effect on such inflammation state of the infected lung cell, and the only form to suppress it is with the inhibition of the activity of the viral protein Nsp5.To that end, our model suggest that drugs like Saquinavir may be useful. In this form, our model suggests that Nsp5 is effectively a central node underlying the severe acute lung inflammation during SARS-CoV-2 infection. The persistent production of IL-6 by lung cells can be one of the causes of the cytokine storm observed in critical patients with COVID19. Nsp5 seems to be the switch to start inflammation, the consequent overproduction of the ACE2 receptor, and an important underlying cause of the most severe cases of COVID19.

In recent years, applying deep learning to solve physics problems has attracted much attention. Data-driven deep learning methods produce operators that can learn solutions to the whole system of partial differential equations. However, the existing methods are only evaluated on simple flow equations (e.g., Burger's equation), and only consider the generalization ability on different initial conditions. In this paper, we construct CFDBench, a benchmark with four classic problems in computational fluid dynamics (CFD): lid-driven cavity flow, laminar boundary layer flow in circular tubes, dam flows through the steps, and periodic Karman vortex street. Each flow problem includes data with various boundary conditions, fluid physical properties, and domain geometries. Compared to existing datasets, the advantages of CFDBench are (1) comprehensive. It contains common physical parameters such as velocity, pressure, and cavity fraction. (2) realistic. It is very suitable for deep learning solutions of fluid mechanics equations. (3) challenging. It has a certain learning difficulty, prompting to find models with strong learning ability. (4) standardized. CFDBench facilitates a comprehensive and fair comparison of different deep learning methods for CFD. We make appropriate modifications to popular deep neural networks to apply them to CFDBench and enable the accommodation of more changing inputs. The evaluation on CFDBench reveals some new shortcomings of existing works and we propose possible directions for solving such problems.

Modeling brain dynamics to better understand and control complex behaviors underlying various cognitive brain functions are of interests to engineers, mathematicians, and physicists from the last several decades. With a motivation of developing computationally efficient models of brain dynamics to use in designing control-theoretic neurostimulation strategies, we have developed a novel data-driven approach in a long short-term memory (LSTM) neural network architecture to predict the temporal dynamics of complex systems over an extended long time-horizon in future. In contrast to recent LSTM-based dynamical modeling approaches that make use of multi-layer perceptrons or linear combination layers as output layers, our architecture uses a single fully connected output layer and reversed-order sequence-to-sequence mapping to improve short time-horizon prediction accuracy and to make multi-timestep predictions of dynamical behaviors. We demonstrate the efficacy of our approach in reconstructing the regular spiking to bursting dynamics exhibited by an experimentally-validated 9-dimensional Hodgkin-Huxley model of hippocampal CA1 pyramidal neurons. Through simulations, we show that our LSTM neural network can predict the multi-time scale temporal dynamics underlying various spiking patterns with reasonable accuracy. Moreover, our results show that the predictions improve with increasing predictive time-horizon in the multi-timestep deep LSTM neural network.

Transport of hydrophobic drugs in the human body exhibits complications due to the low solubility of these compounds. With the purpose of enhancing the bioavailability and biodistribution of such drugs, recent studies have reported the use of amphiphilic molecules, such as phospholipids, for synthesis of nanoparticles or nanocapsules. Given that phospholipids can self–assemble in liposomes or micellar structures, they are ideal candidates to function as vehicles of hydrophobic molecules. In this work, we report mesoscopic simulations of nanoliposomes, constituted by lecithin and coated with a shell of chitosan. The stability of such structure and the efficiency of encapsulation of capsaicin, as well as the internal and superficial distribution of capsaicin and chitosan inside the nanoliposome were analyzed. The characterization of the system was carried out through density maps and the potentials of mean force for the lecithin–capsaicin, lecithin–chitosan and capsaicin–chitosan interactions. The results of these simulations show that chitosan is deposited on the surface of the nanoliposome, as has been reported in some experimental works. It was also observed that a nanoliposome of approximately 18 nm in diameter is stable during the simulation. The deposition behavior was found to be influenced by pattern of N-acetylation of chitosan.

The details of antigen-antibody interactions and the identification of epitopes are critical for the development of monoclonal antibody drugs. Ab42 is a native human-derived anti-CFH monoclonal antibody. In this study, the interaction between antigen pCFH and antibody (Ab42) was theoretically demonstrated by molecular docking and MD simulation, combined with free energy calculation and computational alanine scanning (CAS), and key amino acids and epitopes were identified. Experimental alanine scanning (EAS) was then carried out to verify the results of the calculation, and our results indicated that Ab42 antibody forms hydrogen bonds and interacts hydrophobically with pCFH through the Tyr315, Ser100, Gly33, and Tyr53 residues on its CDR, while the main pCFH epitopes are located at the six sites of Pro441, Ile442, Asp443, Asn444, Ile447, and Thr448. In conclusion, this study has explored the mechanism of antigen-antibody interaction from both theoretical and experimental aspects, and our results have important theoretical significance for the design and development of relevant antibody drugs.

This work is a top-level summary of several contributions published in the last three decades. It makes the case that complex dynamics of nonlinear systems lies at the heart of foundational physics.

Osteoarthritis (OA) is a common joint disorder found mostly in elderly people. The role of mechanical behavior in the progression of OA is complex and remains unclear. The stress-relaxation behavior of human articular cartilage in clinically defined osteoarthritic stages may have importance in diagnosis and prognosis of OA. In this study we investigated differences in the biomechanical responses among human cartilage of ICRS grades I, II and III using polymer dynamics theory. We collected 24 explants of human articular cartilage (eight each of ICRS grade I, II and III) and acquired stress-relaxation data applying a continuous load on the articular surface of each cartilage explant for 1180 s. We observed a significant decrease in Young’s modulus, stress-relaxation time, and stretching exponent in advanced stages of OA (ICRS grade III). The stretch exponential model indicated that significant loss in hyaluronic acid polymer might be the reason for the loss of proteoglycan in advanced OA. This work encourages further biomechanical modelling of osteoarthritic cartilage utilizing these data as input parameters to enhance the fidelity of computational models aimed at revealing how mechanical behaviors play a role in pathogenesis of OA.

The active power-assist function greatly expands the potential applications of exoskeleton robots, yet the motion intention estimation (MIE) for active power-assist strategy is quite problematic. Through the analysis of the conduction path and the different stage manifestations of motion intention in human body, we confirmed that the joint torque of human body meets the basic requirements of MIE for the active power-assist that we suggest, namely: (i) it reflects the direction and intensity of the wearer’s efforts; (ii) it precedes the human limb motion; (iii) it generates real-time and continuous output. Thus, an online calculation method of human joint torque was proposed. The sensing system integrated in exoskeleton robots was designed to perceive motion data and foot contact force of a human body. A special inverse dynamics with a parameterized model of the human body was proposed. Contrast experiments were carried out with the motion capture system, which results’ accuracy and similarity were evaluated via the root mean square error and correlation coefficient. The comparative analysis of two synchronous results shows good accuracy of the proposed MIE method, which lays the foundation for the realization of active power-assist.

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.

Our earlier work has tentatively shown that the hierarchy of fermion masses and mixing angles follows from the universal behavior of nonlinear dynamics. In this brief sequel we survey a similar scenario in which the Cabibbo angle arises from the nonlinear dynamics of charged-current interactions.

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.

In this research, atomistic molecular dynamics simulations are combined with mesoscopic phase-field computational methods in order to investigate phase-transformation in polycrystalline Aluminum microstructure. In fact, microstructural computational modeling of engineering materials could help to optimize their mechanical properties for industrial applications (e.g. directional solidification for turbine blades). As a result, a multiscale modeling approach is developed to find a relation between manufacturing variables (e.g. temperature) and microstructural properties of crystalline materials (e.g. grain size), which could be used to develop an advanced manufacturing process for sensitive applications. The results show that atomistic modeling of grain growth could be used as a first-principle approach in order to study phase transformation's kinetics, which could capture morphology of polycrystalline materials more accurately. On the other hand, phase-field mesoscopic approach needs less computational efforts, but still it relies on semi-empirical data to capture accurate phase transformation regimes, which makes this approach suitable for rapid examining of new manufacturing conditions as well as its effects on microstructural properties of polycrystalline materials.

Markov chain Monte Carlo sampling propagators, including numerical integrators for stochastic dynamics, are central to the calculation of thermodynamic quantities and determination of structure for molecular systems. Efficiency is paramount, and to a great extent, this is determined by the integrated autocorrelation time (IAcT). This quantity varies depending on the observable that is being estimated. It is suggested that it is the maximum of the IAcT over all observables that is the relevant metric. Reviewed here is a method for estimating this quantity. For reversible propagators (which are those that satisfy detailed balance), the maximum IAcT is determined by the spectral gap in the forward transfer operator, but for irreversible propagators, the maximum IAcT can be far less than or greater than what might be inferred from the spectral gap. This is consistent with recent theoretical results (not to mention past practical experience) suggesting that irreversible propagators generally perform better if not much better than reversible ones. Typical irreversible propagators involve a parameter controlling the mix of ballistic and diffusive movement. To gain insight into the effect of the damping parameter for Langevin dynamics, its optimal value is obtained here for a multidimensional quadratic potential energy function.

Three cases of Day-Night Band (DNB) observations of the Visible Infrared Imaging Radiometer Suite (VIIRS) onboard the Suomi National Polar-orbiting Partnership (SNPP) are explored for applications to assess the ocean environment and monitor ocean dynamics. An approach to use the DNB radiance ratio was developed in order to better continuously assess the ocean diurnal and short-term environmental changes with VIIRS DNB observations. In the La Plata River Estuary, the sediment fronts showed 20–25 km diurnal inshore-offshore movements on March 13, 2017. In the waters off the Argentina coast in the South Atlantic, VIIRS DNB measurements provided continuous observations and monitoring of the algae bloom development and migration between 24–26 March 2016. This algae bloom generally kept the same spatial patterns, but moved nearly 20 km eastward in the three-day period. In the Yangtze River Estuary and Hangzhou Bay region in China’s east coast, VIIRS DNB observations also revealed the complicated coastal dynamic changes between 12–14 April 2017. Even though there are still some challenges and limitations for monitoring the ocean environment with VIIRS DNB observations, this study shows that satellite DNB observations can provide additional data sources for ocean observations, especially observations during the nighttime.

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.

Borneol and menthol were two terpenes wildly used as penetrate enhancer in transdermal drug delivery. To explore their penetration enhancement effect towards hydrophilic drug, 5-FU was selected as model drug. A method combined vitro permeation studies and coarse grain molecular dynamics was used to investigate their penetration enhancement effect towards 5-FU. As a result, although both borneol and menthol showed a penetration enhancement effect towards 5-FU, they differed a lot in the penetration enhancement mechanism, which was also thought to account for their different penetration enhancement effect. As for menthol, SC bilayer disrupting effect seemed to be its main mechanism. While for borneol, its mechanism seemed to be more complicated. Except for disrupting the SC bilayer, it could also increase the permeation of 5-FU by enhancing the diffusion rate of 5-FU or inducing the formation of transient pore. All of this enable us a molecular understanding of borneol and menthol’s penetration enhancement effect towards hydrophilic drug, which might provide some guidance in the latter research and application.

The fundamental question in the field of sensor research is new directions of scientific fields, which play a vital role in the progress of science and technology. This study confronts this question here by developing a bibliometric analysis, which endeavors to explain the evolution of sensor research and new technologies that are critical to science and society. The database of Scopus concerning scientific documents and patents is used for statistical and computational analyses in these topics. Results suggest that emerging technological trajectories in sensors are wireless sensor networks, wearable sensors and biosensors. Main characteristics of these growing research fields and technologies in sensors are described for fruitful implications of research and innovation policy directed to science advances and technological change in society.

Relationships among the population abundance of four predator groups for rice insect pests, namely: carabid beetles, staphylinid beetles, green mirid bugs, and spiders in three landscape categories were evaluated. Both rice plots and the associated bund margins of these rice plots found among three Bangladesh landscape categories were sampled by sweep net. The results revealed that the abundance significantly varied across landscapes. The rice landscape of one location harbored higher numbers of a specific predator than other location in other regions of Bangladesh. The results also showed a dependency on the width of the rice bund margins of the rice plots, where spiders populations increased with increased bund widths, but the population abundance of these predators did not depend on the diversity of the number of weed species found on the rice bund margins. The relative abundance of predator populations also significantly differed among the three landscapes, with the green mirid bug having the highest number among the four predators. This study indicates that predators of rice insect pests are highly landscape specific. In order to design integrated pest management systems for different Bangladeshi rice production locales, considerations unique to the characteristics of each locale are necessary. Preliminary efforts to apply variography analyses to the RED spectral band of LANDSAT 8 imagery from December 2016 are presented as first step toward learning a suite of methods which describe useful local characteristics affecting rice pest predators.

The spread of ideas is a fundamental concern of today’s news ecology. Understanding the dynamics of the spread of information and its co-option by interested parties is of critical importance. Research on this topic has shown that individuals tend to cluster in echo-chambers and are driven by confirmation bias. In this paper, we leverage the active inference framework to provide an in silico model of confirmation bias and its effect on echo-chamber formation. We build a model based on active inference, where agents tend to sample information in order to justify their own view of reality, which eventually leads to them to have a high degree of certainty about their own beliefs. We show that, once agents have reached a certain level of certainty about their beliefs, it becomes very difficult to get them to change their views. This system of self-confirming beliefs is upheld and reinforced by the evolving relationship between agent's beliefs and its observations, which over time will continue to provide evidence for their ingrained ideas about the world. The epistemic communities that are consolidated by these shared beliefs, in turn, tend to produce perceptions of reality that reinforce those shared beliefs. We provide an active inference account of this community formation mechanism. We postulate that agents are driven by the epistemic value that they obtain from sampling or observing the behaviors of other agents. Inspired by digital social networks like Twitter, we build a generative model in which agents generate observable social claims or posts (e.g. `tweets') while reading the socially-observable claims of other agents, that lend support towards one of two mutually-exclusive abstract topics. Agents can choose which other agent they pay attention to at each timestep, and crucially who they attend to and what they choose to read influences their beliefs about the world. Agents also assess their local network’s perspective, influencing which kinds of posts they expect to see other agents making. The model was built and simulated simulated using the freely-available Python package pymdp. The proposed active inference model can reproduce the formation of echo-chambers over social networks, and gives us insight into the cognitive processes that lead to this phenomenon.

When gravity is quantum, the point structure of space-time should be replaced by a non-commutative geometry. This is true even for quantum gravity in the infra-red. Using the octonions as space-time coordinates, we construct a pre-spacetime, pre-quantum Lagrangian dynamics. We show that the symmetries of this non-commutative space unify the standard model of particle physics with $SU(2)_R$ chiral gravity. The algebra of the octonionic space yields spinor states which can be identified with three generations of quarks and leptons. The geometry of the space implies quantisation of electric charge, and leads to a theoretical derivation of the mysterious mass ratios of quarks and the charged leptons. Quantum gravity is quantisation not only of the gravitational field, but also of the point structure of space-time.

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 this work, a tool for estimating the contact angle from the molecular dynamics simulations is developed and presented. The tool (Achilles) can detect water droplet on hydrophobic and hydrophilic surfaces. The tool can reconstruct the droplets broken across the periodic boundaries. Further a neighbor density based accurate filter is used to find the droplet liquid vapor interface and a circle is fitted using it after removing the dense layers of water next to solid surface. This fitted circle is solved for contact angle and results are outputted in the form of graphical images and text. The entire content of the internal computations of the tool is broken down into 4 phases and users can monitor the outcomes at every phase through output images. The tool is tested using sample molecular dynamics results of water droplet on hydrophobic and hydrophilic surfaces. We believe this tool can be a good addition to the molecular dynamics simulation community who work on the interfacial physics, droplet evaporation, super hydrophobic surfaces, and wettability etc.

The Autopoiesis and Cognition Theory (ACT), by Maturana and Varela, based on the notions of Biological Closure and Structural Coupling, is a well-known theory on how to understand biological organization [1, 2, 3]. Although, for example, the Free Energy Principle framework evokes some entailments of autopoiesis in a more formal setting [4, 5]; and ACT has been used in many fields, its impact has been restricted because it lacks quantitative analysis. Here we present a theoretical framework grounded in accepted and well-developed ideas from Mathematics and Physics which advance the understanding of the Principles of Biological Organization under the guidance of Biological Closure and Structural Coupling. The disciplines of Differential Geometry/Topology, Mechanics and Complex Dynamical Systems provide a powerful, elegant, and well-established body of knowledge to support our Biological Organization Principles (BOP) framework. In particular, Stochastic Mechanics and KAM theory (from Kolmogorov, Arnold and Moser theorem) allow us to develop, using the notions of Biological Closure and Structural Coupling, a central core of BOP termed Dynamical Closure Mechanism. Under the proposed framework, a wide variety of bio- logical phenomena can be understood, shedding new light on biological explanations. However, an understanding of biological organization may require the re-evaluation of dogmas on how we think on biology as it seems inescapable that what is needed is an integration of analysis and notions derived from mathematics, physics, and biology to generate a new landscape of ideas.

Micro-pollutants such as 17β-Estradiol (E2) have been detected in different water resources and their negative effects on the environment and organisms have been observed. Aptamers are established as a possible detection tool, but the underlying ligand binding is largely unexplored. In this study, a previously described 35-mer E2-specific aptamer was used to analyse the binding characteristics between E2 and the aptamer with a MD simulation in an aqueous medium. Because there is no 3D structure information available for this aptamer, it was modeled using coarse-grained modeling method. The E2 ligand was positioned inside a potential binding area of the predicted aptamer structure, the complex was used for an 25 ns MD simulation, and the interactions were examined for each time step. We identified E2-specific bases within the interior loop of the aptamer and also demonstrated the influence of frequently underestimated water-mediated hydrogen bonds. The study contributes to the understanding of the behavior of ligands binding with aptamer structure in an aqueous solution. The developed workflow allows generating and examining further appealing ligand-aptamer complexes.

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.

In this paper we develop from first principles a unique law pertaining to the flow of fluids through closed conduits. This law, which we call “Quinn’s Law”, may be described as follows: When fluids are forced to flow through closed conduits under the driving force of a pressure gradient, there is a linear relationship between the normalized dimensionless pressure gradient, P_Q, and the normalized dimensionless fluid current, C_Q. The relationship is expressed mathematically as: P_Q = k₁ +k₂C_Q. This linear relationship remains the same whether the conduit is filled with or devoid of solid obstacles. The law differentiates, however, between a packed and an empty conduit by virtue of the tortuosity of the fluid path, which is seamlessly accommodated within the normalization framework of the law itself. When movement of the fluid is very close to being at rest, i.e., very slow, this relationship has the unique minimum constant value of k₁, and as the fluid acceleration increases, it varies with a slope of k₂ as a function of normalized fluid current. Quinn’s Law is validated herein by applying it to the data from published classical studies of measured permeability in both packed and empty conduits, as well as to the data generated by home grown experiments performed in the author’s own laboratory.

Various distance metrics and their induced norms are employed in the quantitative modeling of evolutionary dynamics. Minimization of these distance metrics when applied to evolutionary optimization are hypothesized to result in different outcomes. Here, we apply the different distance metrics to the evolutionary trait dynamics brought about by the interaction between two competing species infected by parasites (exploiters). We present deterministic cases showing the distinctive selection outcomes under the Manhattan, Euclidean and Chebyshev norms. Specifically, we show how they differ in the time of convergence to the desired optima (e.g., no disease), and in the egalitarian sharing of carrying capacity between the competing species. However, when randomness is introduced to the population dynamics of parasites and to the trait dynamics of the competing species, the distinctive characteristics of the outcomes under the three norms become indistinguishable. Our results provide theoretical cases when evolutionary dynamics using different distance metrics exhibit similar outcomes.

The problem of fluid dynamics can be greatly simplified if, for every point in space, the strain-rate tensor is diagonalized. This tensor is introduced into the Navier-Stokes equations via material law and divergence of the stress tensor. This article shows that local SO(3)xU(1) gauge fields can be used to locally diagonalize the diffusion components of the strain-rate tensor. The gauge fields resulting from the connection can be interpreted as convection components of the flow, they show properties of quasiparticles and can be interpreted as elementary vortices. Thus, the proposed approach not only offers new insights for the solution and situative simplification of the Navier-Stokes equations, it also uncovers hidden symmetries within the flow convection, allowing - depending on boundary conditions - further interpretation.

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.

Histone Lysine Specific Demethylase 1 (LSD1) is overexpressed in many cancers and become a new target for anticancer drugs. In recent years, the small molecule inhibitors with various structures targeting LSD1 have been reported. Here we report the binding interaction modes of a series of thieno[3,2-b]pyrrole-5-carboxamides LSD1 inhibitors using molecular docking, three dimensional quantitative structure-activity relationship (3D-QSAR). Comparative molecular field analysis (CoMFA q²=0.783, r²=0.944, r²_pred=0.851) and Comparative molecular similarity indices analysis (CoMSIA q²=0.728, r²=0.982, r²_pred=0.814) were used to establish 3D-QSAR models, which had good verification and prediction capabilities. Based on the contour maps and the information of molecular docking, 8 novel small molecules were designed in silico, among which compounds D4, D5 and D8 with high predictive activity were subjected to further molecular dynamics simulations (MD), and their possible binding modes were explored. It was found that Asn535 plays a crucial role in stabilizing the inhibitors. Furthermore, the ADME and bioavailability prediction for D4, D5 and D8 were carried out. The results would provide valuable guidance for designing new reversible LSD1 inhibitors in the future.

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.

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.

Flavonoids are widely recognized as natural polydrugs, given their anti-inflammatory, antioxidant, sedative and antineoplastic activities. Recently, different studies have shown that flavonoid have the potential to inhibit BET bromodomains. Previous reports suggest that flavonoids are putative inhibitors of the ZA channel due to their orientation and interactions with P86, V87, L92, L94 and N140. Herein, a comprehensive characterization of the binding mode of the biflavonoid amentoflavone and fisetin is discussed. To this end, both compounds were docked with BRD4 using four docking programs. Results were post-processed with protein-ligand interaction fingerprints. To gain further insights into the binding mode of the two natural products, docking results were further analyzed with molecular dynamics. Results showed that amentoflavone makes numerous contacts in the ZA channel, as previously described for flavonoids and kinase inhibitors. It was also found that amentoflavone can potentially make contacts with non-canonical residues for BET inhibition. Most of these contacts were not observed with fisetin. Based on these results, amentoflavone was tested for BRD4 inhibition, showing activity in the micromolar range. This work may serve as basis for scaffold optimization and further characterization of flavonoids as BET inhibitors.

The maximization of entropy S within a closed system is accepted as an inevitability (as the second law of thermodynamics) by statistical inference alone. The Maximum Entropy Production Principle (MEPP) states that not only S maximizes, but $\dot{S}$ as well: a system will dissipate as fast as possible. There is still no consensus on the general validity of this MEPP, even though it shows remarkable explanatory power (both qualitatively and quantitatively), and has been empirically demonstrated for many domains. In this theoretical paper I provide a generalization of entropy gradients, to show that the MEPP actually follows from the same statistical inference, as that of the 2nd law of thermodynamics. For this generalization I only use the concepts of super-statespaces and microstate-density. These concepts also allow for the abstraction of 'Self Organizing Criticality' to a bifurcating local difference in this density, and allow for a generalization of the fundamentally unresolved concepts of 'chaos' and 'order'.

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.