A stochastic thermostat algorithm for coarse-grained thermomechanical modeling of large-scale soft matters : theory and application to microfilaments

As all-atom molecular dynamics method is limited by its enormous computational cost, various coarse-grained strategies have been developed to extend the length scale of soft matters in the modeling of mechanical behaviors. However, the classical thermostat algorithm in highly coarse-grained molecular dynamics method would underestimate the thermodynamic behaviors of soft matters (e.g. microfilaments in cells), which can weaken the ability of materials to overcome local energy traps in granular modeling. Based on all-atom molecular dynamics modeling of microfilament fragments (G-actin clusters), a new stochastic thermostat algorithm is developed to retain the representation of thermodynamic properties of microfilaments at extra coarse-grained level. The accuracy of this stochastic thermostat algorithm is validated by all-atom MD simulation. This new stochastic thermostat algorithm provides an efficient way to investigate the thermomechanical properties of large-scale soft matters.

A Microtube Viscometer with a Thermostat

The viscometer presented in this paper is suit-able for measuring the viscosity of liquids in micro-litre quantities. It consists of a micro-flow experimental system with a thermostat. Using the measurements of the flow rates and pressure drops of a liquid passing through a microtube, the liquid's viscosity can be calculated from on Hagen-Poiseuille theory. After calibration, the viscometer was used to measure viscosities of deionized water and ethyl alcohol at temperatures ranging from 0 to 40 "C. For both test liquids, the relative deviation of the measured values from those quoted in the literature (obtained using other viscometers) was less than 2.6o/o. The relative uncertainty of the experimental system was reduced to +-l.8% using the relative measuring method. Due to the micro-scale of the test section, only a micro-litre quantity of liquid is needed for a test, this is a potential advantage for measurement of bio-liquid viscosities.

An isothermal-isobaric Langevin thermostat for simulating nanoparticles under pressure: Application to Au clusters

We present a method for simulating clusters or, molecules subjected to an external pressure, which is exerted by a pressure-transmitting medium. It is based on the canoninical Langevin thermostat, but extended in such a way that the Brownian forces are allowed to operate only from the region exterior to the cluster. We show that the frictional force of the Langevin thermostat is linked to the pressure of the reservoir in a unique way, and that this property manifests itself when the particle it acts upon is not pointlike but has finite dimensions. By choosing appropriately the strength of the random forces and the friction coefficient, both temperature and pressure can be controlled independently. We illustrate the capabilities of this new method by calculating the compressibility of small gold clusters under pressure.

Effect of the thermostat in the molecular dynamics simulation on the folding of the model protein chignolin

Molecular dynamics simulations of the model protein chignolin with explicit solvent were carried out, in order to analyze the influence of the Berendsen thermostat on the evolution and folding of the peptide. The dependence of the peptide behavior on temperature was tested with the commonly employed thermostat scheme consisting of one thermostat for the protein and another for the solvent. The thermostat coupling time of the protein was increased to infinity, when the protein is not in direct contact with the thermal bath, a situation known as minimally invasive thermostat. In agreement with other works, it was observed that only in the last situation the instantaneous temperature of the model protein obeys a canonical distribution. As for the folding studies, it was shown that, in the applications of the commonly utilized thermostat schemes, the systems are trapped in local minima regions from which it has difficulty escaping. With the minimally invasive thermostat the time that the protein needs to fold was reduced by two to three times. These results show that the obstacles to the evolution of the extended peptide to the folded structure can be overcome when the temperature of the peptide is not directly controlled.

Adjusting building thermostats for environmental gains : understanding the issues

There has been increasing reliance on mechanical heating, ventilation and air-conditioning (HVAC) systems to achieve thermal comfort in office buildings. The use of universal standards for thermal comfort adopted in air-conditioned spaces often results in a large disparity between mean daily external summer temperatures and temperatures experienced indoors. The extensive overuse of air-conditioning in warm climates not only isolates us from the vagaries of the external environment, but is generally dependent on non-renewable energy. A pilot study conducted at the Queensland University of Technology (QUT) involved altering the thermostat set-points to two or three degrees above the normal summer setting in two air-conditioned buildings during the subtropical summer. This paper presents the findings of the research that led to the formulation of the test study. The findings of the test study are printed in the companion paper DES 72: Adjusting Building Thermastats for Environmental Gains – a Pilot Study.

The ECOS Project. An interactive simulation to promote an ECO-friendly universe.

Custom designed for display on the Cube Installation situated in the new Science and Engineering Centre (SEC) at QUT, the ECOS project is a playful interface that uses real-time weather data to simulate how a five-star energy building operates in climates all over the world. In collaboration with the SEC building managers, the ECOS Project incorporates energy consumption and generation data of the building into an interactive simulation, which is both engaging to users and highly informative, and which invites play and reflection on the roles of green buildings. ECOS focuses on the principle that humans can have both a positive and negative impact on ecosystems with both local and global consequence. The ECOS project draws on the practice of Eco-Visualisation, a term used to encapsulate the important merging of environmental data visualization with the philosophy of sustainability. Holmes (2007) uses the term Eco-Visualisation (EV) to refer to data visualisations that ‘display the real time consumption statistics of key environmental resources for the goal of promoting ecological literacy’. EVs are commonly artifacts of interaction design, information design, interface design and industrial design, but are informed by various intellectual disciplines that have shared interests in sustainability. As a result of surveying a number of projects, Pierce, Odom and Blevis (2008) outline strategies for designing and evaluating effective EVs, including ‘connecting behavior to material impacts of consumption, encouraging playful engagement and exploration with energy, raising public awareness and facilitating discussion, and stimulating critical reflection.’ Consequently, Froehlich (2010) and his colleagues also use the term ‘Eco-feedback technology’ to describe the same field. ‘Green IT’ is another variation which Tomlinson (2010) describes as a ‘field at the juncture of two trends… the growing concern over environmental issues’ and ‘the use of digital tools and techniques for manipulating information.’ The ECOS Project team is guided by these principles, but more importantly, propose an example for how these principles may be achieved. The ECOS Project presents a simplified interface to the very complex domain of thermodynamic and climate modeling. From a mathematical perspective, the simulation can be divided into two models, which interact and compete for balance – the comfort of ECOS’ virtual denizens and the ecological and environmental health of the virtual world. The comfort model is based on the study of psychometrics, and specifically those relating to human comfort. This provides baseline micro-climatic values for what constitutes a comfortable working environment within the QUT SEC buildings. The difference between the ambient outside temperature (as determined by polling the Google Weather API for live weather data) and the internal thermostat of the building (as set by the user) allows us to estimate the energy required to either heat or cool the building. Once the energy requirements can be ascertained, this is then balanced with the ability of the building to produce enough power from green energy sources (solar, wind and gas) to cover its energy requirements. Calculating the relative amount of energy produced by wind and solar can be done by, in the case of solar for example, considering the size of panel and the amount of solar radiation it is receiving at any given time, which in turn can be estimated based on the temperature and conditions returned by the live weather API. Some of these variables can be altered by the user, allowing them to attempt to optimize the health of the building. The variables that can be changed are the budget allocated to green energy sources such as the Solar Panels, Wind Generator and the Air conditioning to control the internal building temperature. These variables influence the energy input and output variables, modeled on the real energy usage statistics drawn from the SEC data provided by the building managers.

Cross-scale biophysics modelling of F-actin cytoskeleton in cell

This thesis is a comparative study of the modelling of mechanical behaviours of F-actin cytoskeleton which is an important structural component in living cells. A new granular model was developed for F-actin cytoskeleton based on the concept of multiscale modelling. This framework overcomes difficulties encountered in physical modelling of cytoskeleton in conventional continuum mechanics modelling, and the computational challenges in all-atom molecular dynamics simulation. The thermostat algorithm was further modified to better predict the thermodynamic properties of F-actin cytoskeleton in modelling. This multiscale modelling framework was applied in explaining the physical mechanisms of cytoskeleton responses to external mechanical loads.

Effects of Isothermal and Adiabatic Thermal Loadings on Size and Strain Rate Dependence of Copper Nanowire

In the present paper, the size and strain rate effects on ultra-thin < 100 >/{100} Cu nanowires at an initial temperature of 10 K have been discussed. Extensive molecular dynamics (MD) simulations have been performed using Embedded atom method (EAM) to investigate the structural behaviours and properties under high strain rate. Velocity-Verlet algorithm has been used to solve the equation of motions. Two different thermal loading cases have been considered: (i) Isothermal loading, in which Nose-Hoover thermostat is used to maintain the constant system temperature, and (ii) Adiabatic loading, i.e., without any thermostat. Five different wire cross-sections were considered ranging from 0.723 x 0.723 nm(2) to 2.169 x 2.169 nm(2) The strain rates used in the present study were 1 x 10(9) s(-1), 1 x 10(8) s(-1), and 1 x 10(7) s(-1). The effect of strain rate on the mechanical properties of copper nanowires was analysed, which shows that elastic properties are independent of thermal loading for a given strain rate and cross-sectional dimension of nanowire. It showed a decreasing yield stress and yield strain with decreasing strain rate for a given cross- section. Also, a decreasing yield stress and increasing yield strain were observed for a given strain rate with increasing cross-sectional area. Elastic modulus was found to be similar to 100 GPa, which was independent of processing temperature, strain rate, and size for a given initial temperature. Reorientation of < 100 >/{100} square cross-sectional copper nanowire into a series of stable ultra-thin Pentagon copper nanobridge structures with dia of similar to 1 nm at 10 K was observed under high strain rate tensile loading. The effect of isothermal and adiabatic loading on the formation of such pentagonal nanobridge structure has been discussed.

The changes in various iron compounds in culture media for algae. [Translation from: Archiv fur Hydrobiologie Supplemente Band 38(1/2) 151-169, 1970.]

The purpose of this work is a contribution to the quantitative record of the use of iron by planktonic algae. Preliminary experiments with Chlorella to determine the rate of iron intake in the presence of inorganic sources of iron did not produce the desired result. The crucial point of this work is the investigation of the influence of various external factors on the stability of FeEDTA (FeEDTA = Ferric(III)-compound of ethylene-diamine tetra-acetic acid), since this compound appears to be particularly well-suited as a source of iron for planktonic algae (e.g. TAMIYA et al. 1953). Cultures of Chlorella fusca in a light thermostat were used in experimental research. Methods and results are discussed.

Generalized Langevin equation: An efficient approach to nonequilibrium molecular dynamics of open systems

The generalized Langevin equation (GLE) has been recently suggested to simulate the time evolution of classical solid and molecular systems when considering general nonequilibrium processes. In this approach, a part of the whole system (an open system), which interacts and exchanges energy with its dissipative environment, is studied. Because the GLE is derived by projecting out exactly the harmonic environment, the coupling to it is realistic, while the equations of motion are non-Markovian. Although the GLE formalism has already found promising applications, e. g., in nanotribology and as a powerful thermostat for equilibration in classical molecular dynamics simulations, efficient algorithms to solve the GLE for realistic memory kernels are highly nontrivial, especially if the memory kernels decay nonexponentially. This is due to the fact that one has to generate a colored noise and take account of the memory effects in a consistent manner. In this paper, we present a simple, yet efficient, algorithm for solving the GLE for practical memory kernels and we demonstrate its capability for the exactly solvable case of a harmonic oscillator coupled to a Debye bath.

Applications of the generalized Langevin equation: Towards a realistic description of the baths

The generalized Langevin equation (GLE) method, as developed previously [L. Stella et al., Phys. Rev. B 89, 134303 (2014)], is used to calculate the dissipative dynamics of systems described at the atomic level. The GLE scheme goes beyond the commonly used bilinear coupling between the central system and the bath, and permits us to have a realistic description of both the dissipative central system and its surrounding bath. We show how to obtain the vibrational properties of a realistic bath and how to convey such properties into an extended Langevin dynamics by the use of the mapping of the bath vibrational properties onto a set of auxiliary variables. Our calculations for a model of a Lennard-Jones solid show that our GLE scheme provides a stable dynamics, with the dissipative/relaxation processes properly described. The total kinetic energy of the central system always thermalizes toward the expected bath temperature, with appropriate fluctuation around the mean value. More importantly, we obtain a velocity distribution for the individual atoms in the central system which follows the expected canonical distribution at the corresponding temperature. This confirms that both our GLE scheme and our mapping procedure onto an extended Langevin dynamics provide the correct thermostat. We also examined the velocity autocorrelation functions and compare our results with more conventional Langevin dynamics.

Molecular dynamics study on the equilibrium magnetization properties and structure of ferrofluids

We investigate in detail the initial susceptibility, magnetization curves, and microstructure of ferrofluids in various concentration and particle dipole moment ranges by means of molecular dynamics simulations. We use the Ewald summation for the long-range dipolar interactions, take explicitly into account the translational and rotational degrees of freedom, coupled to a Langevin thermostat. When the dipolar interaction energy is comparable with the thermal energy, the simulation results on the magnetization properties agree with the theoretical predictions very well. For stronger dipolar couplings, however, we find systematic deviations from the theoretical curves. We analyze in detail the observed microstructure of the fluids under different conditions. The formation of clusters is found to enhance the magnetization at weak fields and thus leads to a larger initial susceptibility. The influence of the particle aggregation is isolated by studying ferro-solids, which consist of magnetic dipoles frozen in at random locations but which are free to rotate. Due to the artificial suppression of clusters in ferrosolids the observed susceptibility is considerably lowered when compared to ferrofluids.

Passive solar options for reducing heating demand and maintain indoor climate in a multifamily house in Sweden

This research was carried out by studying possible renovation of a two-storey detached multifamily building by using passive solar design options in a cold climate in Borlänge, Sweden where the heating Degree Days are 4451 (base 20°C). Borlänge`s housing company, Tunabyggen, plans to renovate the project house located inthe multicultural district, Jakobsgårdarna. The goal of the thesis was to suggest a redesign of the current building, decrease the heating energy use, by applying passive solar design and control strategies, in a most reasonable way. In addition ensure a better thermal comfort for the tenants in the dwellings. Literatures have been studied, from which can be inferred that passive design should be abasic design consideration for all housing constructions, because it has advantages to ensure thermal comfort, and reduce the energy use. In addition further savings can be achieved applying different types of control strategies, from which the house will be more personalized, and better adapted to the user’s needs.The proposed method is based on simulations by using TRNSYS software. First a proper building model was set up, which represents the current state of the project building. Then the thermal insulation and the windows were upgraded, based on today's building regulations. The developments of the passive solar options were accomplished in two steps. First of all the relevant basic passive design elements were considered, then those advantages were compared to the advantages of applying new conventional thermostat, and shading control strategies.The results show that there is significant potential with the different types of passive solar design; their usage depends primarily on the location of the site as well as the orientation of the project building. Applying the control strategies, such as thermostat, and shading control, along the thermal insulation upgrade, may lead to significant energy savings (around 40 %), by comparison to the reference building without any upgrade.

Colina em rações para tilápia do Nilo: desempenho produtivo e respostas hematológicas antes e após o estímulo a frio

Pós-graduação em Zootecnia - FMVZ