Avalanche dynamics of fiber bundle models

We present a detailed analytical and numerical study of the avalanche distributions of the continuous damage fiber bundle model CDFBM . Linearly elastic fibers undergo a series of partial failure events which give rise to a gradual degradation of their stiffness. We show that the model reproduces a wide range of mechanical behaviors. We find that macroscopic hardening and plastic responses are characterized by avalanche distributions, which exhibit an algebraic decay with exponents between 5/2 and 2 different from those observed in mean-field fiber bundle models. We also derive analytically the phase diagram of a family of CDFBM which covers a large variety of potential avalanche size distributions. Our results provide a unified view of the statistics of breaking avalanches in fiber bundle models

Electronic couplings and on-site energies for hole transfer in DNA: systematic quantum mechanical/molecular dynamic study

The electron hole transfer (HT) properties of DNA are substantially affected by thermal fluctuations of the π stack structure. Depending on the mutual position of neighboring nucleobases, electronic coupling V may change by several orders of magnitude. In the present paper, we report the results of systematic QM/molecular dynamic (MD) calculations of the electronic couplings and on-site energies for the hole transfer. Based on 15 ns MD trajectories for several DNA oligomers, we calculate the average coupling squares 〈 V2 〉 and the energies of basepair triplets X G+ Y and X A+ Y, where X, Y=G, A, T, and C. For each of the 32 systems, 15 000 conformations separated by 1 ps are considered. The three-state generalized Mulliken-Hush method is used to derive electronic couplings for HT between neighboring basepairs. The adiabatic energies and dipole moment matrix elements are computed within the INDO/S method. We compare the rms values of V with the couplings estimated for the idealized B -DNA structure and show that in several important cases the couplings calculated for the idealized B -DNA structure are considerably underestimated. The rms values for intrastrand couplings G-G, A-A, G-A, and A-G are found to be similar, ∼0.07 eV, while the interstrand couplings are quite different. The energies of hole states G+ and A+ in the stack depend on the nature of the neighboring pairs. The X G+ Y are by 0.5 eV more stable than X A+ Y. The thermal fluctuations of the DNA structure facilitate the HT process from guanine to adenine. The tabulated couplings and on-site energies can be used as reference parameters in theoretical and computational studies of HT processes in DNA

Charge transfer in DNA: hole charge is confined to a single base pair due to solvation effects

We include solvation effects in tight-binding Hamiltonians for hole states in DNA. The corresponding linear-response parameters are derived from accurate estimates of solvation energy calculated for several hole charge distributions in DNA stacks. Two models are considered: (A) the correction to a diagonal Hamiltonian matrix element depends only on the charge localized on the corresponding site and (B) in addition to this term, the reaction field due to adjacent base pairs is accounted for. We show that both schemes give very similar results. The effects of the polar medium on the hole distribution in DNA are studied. We conclude that the effects of polar surroundings essentially suppress charge delocalization in DNA, and hole states in (GC)n sequences are localized on individual guanines

Soil water dynamics related to the degree of compaction of two brazilian oxisols under no-tillage

Soil water properties are related to crop growth and environmental aspects and are influenced by the degree of soil compaction. The objective of this study was to determine the water infiltration and hydraulic conductivity of saturated soil under field conditions in terms of the compaction degree of two Oxisols under a no-tillage (NT). Two commercial fields were studied in the state of Rio Grande do Sul, Brazil: one a Haplortox after 14 years under NT; the other a Hapludox after seven years under NT. Maps (50 x 30 m) of the levels of mechanical penetration resistance (PR) were drawn based on the kriging method, differentiating three compaction degrees (CD): high, intermediate and low. In each CD area, the infiltration rate (initial and steady-state) and cumulative water infiltration were measured using concentric rings, with six replications, and the saturated hydraulic conductivity (K(θs)) was determined using the Guelph permeameter. Statistical evaluation was performed based on a randomized design, using the least significant difference (LSD) test and regression analysis. The steady-state infiltration rate was not influenced by the compaction degree, with mean values of 3 and 0.39 cm h-1 in the Haplortox and the Hapludox, respectively. In the Haplortox, saturated soil hydraulic conductivity was 26.76 cm h-1 at a low CD and 9.18 cm h-1 at a high CD, whereas in the Hapludox, this value was 5.16 cm h-1 and 1.19 cm h-1 for the low and high CD, respectively. The compaction degree did not affect the initial and steady-state water infiltration rate, nor the cumulative water infiltration for either soil type, although the values were higher for the Haplortox than the Hapludox.

Force volume and stiffness tomography investigation on the dynamics of stiff material under bacterial membranes.

The determination of the characteristics of micro-organisms in clinical specimens is essential for the rapid diagnosis and treatment of infections. A thorough investigation of the nanoscale properties of bacteria can prove to be a fundamental tool. Indeed, in the latest years, the importance of high resolution analysis of the properties of microbial cell surfaces has been increasingly recognized. Among the techniques available to observe at high resolution specific properties of microscopic samples, the Atomic Force Microscope (AFM) is the most widely used instrument capable to perform morphological and mechanical characterizations of living biological systems. Indeed, AFM can routinely study single cells in physiological conditions and can determine their mechanical properties with a nanometric resolution. Such analyses, coupled with high resolution investigation of their morphological properties, are increasingly used to characterize the state of single cells. In this work, we exploit the capabilities and peculiarities of AFM to analyze the mechanical properties of Escherichia coli in order to evidence with a high spatial resolution the mechanical properties of its structure. In particular, we will show that the bacterial membrane is not mechanically uniform, but contains stiffer areas. The force volume investigations presented in this work evidence for the first time the presence and dynamics of such structures. Such information is also coupled with a novel stiffness tomography technique, suggesting the presence of stiffer structures present underneath the membrane layer that could be associated with bacterial nucleoids.

Relaxation dynamics in suspensions of ferromagnetic particles

We have studied the relaxation dynamics of a dilute assembly of ferromagnetic particles in suspension. A formalism based on the Smoluchowski equation, describing the evolution of the probability density for the directions of the magnetic moment and of the axis of easy magnetization of the particles, has been developed. We compute the rotational viscosity from a Green-Kubo formula and give an expression for the relaxation time of the particles which comes from the dynamic equations of the correlation functions. Concerning the relaxation time for the particles, our results agree quite well with experiments performed on different samples of ferromagnetic particles for which the magnetic energy, associated with the interaction between the magnetic moments and the external field, or the energy of anisotropy plays a dominant role.

MMPBSA decomposition of the binding energy throughout a molecular dynamics simulation of Amyloid-Beta (Aß10-35) aggregation

Recent experiments with amyloid-beta (Aß) peptides indicate that the formation of toxic oligomers may be an important contribution to the onset of Alzheimer's disease. The toxicity of Aß oligomers depend on their structure, which is governed by assembly dynamics. However, a detailed knowledge of the structure of at the atomic level has not been achieved yet due to limitations of current experimental techniques. In this study, replica exchange molecular dynamics simulations are used to identify the expected diversity of dimer conformations of Aß10-35 monomers. The most representative dimer conformation has been used to track the dimer formation process between both monomers. The process has been characterized by means of the evolution of the decomposition of the binding free energy, which provides an energetic profile of the interaction. Dimers undergo a process of reorganization driven basically by inter-chain hydrophobic and hydrophilic interactions and also solvation/desolvation processes.

Atomic Force Microscopy Stiffness Tomography on Living Arabidopsis thaliana Cells Reveals the Mechanical Properties of Surface and Deep Cell-Wall Layers during Growth.

Cell-wall mechanical properties play a key role in the growth and the protection of plants. However, little is known about genuine wall mechanical properties and their growth-related dynamics at subcellular resolution and in living cells. Here, we used atomic force microscopy (AFM) stiffness tomography to explore stiffness distribution in the cell wall of suspension-cultured Arabidopsis thaliana as a model of primary, growing cell wall. For the first time that we know of, this new imaging technique was performed on living single cells of a higher plant, permitting monitoring of the stiffness distribution in cell-wall layers as a function of the depth and its evolution during the different growth phases. The mechanical measurements were correlated with changes in the composition of the cell wall, which were revealed by Fourier-transform infrared (FTIR) spectroscopy. In the beginning and end of cell growth, the average stiffness of the cell wall was low and the wall was mechanically homogenous, whereas in the exponential growth phase, the average wall stiffness increased, with increasing heterogeneity. In this phase, the difference between the superficial and deep wall stiffness was highest. FTIR spectra revealed a relative increase in the polysaccharide/lignin content.

Real-time monitoring of protein conformational changes using a nano-mechanical sensor.

Proteins can switch between different conformations in response to stimuli, such as pH or temperature variations, or to the binding of ligands. Such plasticity and its kinetics can have a crucial functional role, and their characterization has taken center stage in protein research. As an example, Topoisomerases are particularly interesting enzymes capable of managing tangled and supercoiled double-stranded DNA, thus facilitating many physiological processes. In this work, we describe the use of a cantilever-based nanomotion sensor to characterize the dynamics of human topoisomerase II (Topo II) enzymes and their response to different kinds of ligands, such as ATP, which enhance the conformational dynamics. The sensitivity and time resolution of this sensor allow determining quantitatively the correlation between the ATP concentration and the rate of Topo II conformational changes. Furthermore, we show how to rationalize the experimental results in a comprehensive model that takes into account both the physics of the cantilever and the dynamics of the ATPase cycle of the enzyme, shedding light on the kinetics of the process. Finally, we study the effect of aclarubicin, an anticancer drug, demonstrating that it affects directly the Topo II molecule inhibiting its conformational changes. These results pave the way to a new way of studying the intrinsic dynamics of proteins and of protein complexes allowing new applications ranging from fundamental proteomics to drug discovery and development and possibly to clinical practice.

Nuclear magnetic resonance investigation of structure and dynamics of some electrolyte solutions

Electrolyte solutions are of importance in a wide range of scientific contexts and as such have attracted considerable theoretical and experimental effort over many years. Nuclear Magnetic resonance provides a precise and versatile tool for investigation of electrolyte solutions, both in water and in organic solvents. Many structural and dynamic properties can be obtained through NMR experiments. The solution of aluminum chloride in water was studied. Different concentrations were taken for investigation. Independence of maximum line shift from concentration and acidity was shown. Six-coordinated structure of solvation shell was confirmed by experiments on 'H and 27A1 nuclei. Diffusion coefficients were studied. The solution of nickel chloride in methanol was studied. Lines, corresponding to coordinated and bulk methanol were found. Four-, five- and six-coordinated structures were found in different temperatures. The line for coordinated -OD group of deuterated methanol was observed on 2H spectrum for the first time. Partial deuteration of CH3 group was detected. Inability to observe coordinated -OH group was explained.

Validation of plant dynamic model by online and laboratory measurements – a tool to predict online COD loads out of production of mechanical printing papers

COD discharges out of processes have increased in line with elevating brightness demands for mechanical pulp and papers. The share of lignin-like substances in COD discharges is on average 75%. In this thesis, a plant dynamic model was created and validated as a means to predict COD loading and discharges out of a mill. The assays were carried out in one paper mill integrate producing mechanical printing papers. The objective in the modeling of plant dynamics was to predict day averages of COD load and discharges out of mills. This means that online data, like 1) the level of large storage towers of pulp and white water 2) pulp dosages, 3) production rates and 4) internal white water flows and discharges were used to create transients into the balances of solids and white water, referred to as “plant dynamics”. A conversion coefficient was verified between TOC and COD. The conversion coefficient was used for predicting the flows from TOC to COD to the waste water treatment plant. The COD load was modeled with similar uncertainty as in reference TOC sampling. The water balance of waste water treatment was validated by the reference concentration of COD. The difference of COD predictions against references was within the same deviation of TOC-predictions. The modeled yield losses and retention values of TOC in pulping and bleaching processes and the modeled fixing of colloidal TOC to solids between the pulping plant and the aeration basin in the waste water treatment plant were similar to references presented in literature. The valid water balances of the waste water treatment plant and the reduction model of lignin-like substances produced a valid prediction of COD discharges out of the mill. A 30% increase in the release of lignin-like substances in the form of production problems was observed in pulping and bleaching processes. The same increase was observed in COD discharges out of waste water treatment. In the prediction of annual COD discharge, it was noticed that the reduction of lignin has a wide deviation from year to year and from one mill to another. This made it difficult to compare the parameters of COD discharges validated in plant dynamic simulation with another mill producing mechanical printing papers. However, a trend of moving from unbleached towards high-brightness TMP in COD discharges was valid.

Advances in modelling of mechatronic systems: the toolset FASIM_C++ for the simulation of vehicle dynamics

The increasing complexity of controller systems, applied in modern passenger cars, requires adequate simulation tools. The toolset FASIM_C++, described in the following, uses complex vehicle models in three-dimensional vehicle dynamics simulation. The structure of the implemented dynamic models and the generation of the equations of motion applying the method of kinematic differentials is explained briefly. After a short introduction in methods of event handling, several vehicle models and applications like controller development, roll-over simulation and real-time-simulation are explained. Finally some simulation results are presented.

Dynamics of coherent vortices in mixing layers using direct numerical and large-eddy simulations

Coherent vortices in turbulent mixing layers are investigated by means of Direct Numerical Simulation (DNS) and Large-Eddy Simulation (LES). Subgrid-scale models defined in spectral and physical spaces are reviewed. The new "spectral-dynamic viscosity model", that allows to account for non-developed turbulence in the subgrid-scales, is discussed. Pseudo-spectral methods, combined with sixth-order compact finite differences schemes (when periodic boundary conditions cannot be established), are used to solve the Navier- Stokes equations. Simulations in temporal and spatial mixing layers show two types of pairing of primary Kelvin-Helmholtz (KH) vortices depending on initial conditions (or upstream conditions): quasi-2D and helical pairings. In both cases, secondary streamwise vortices are stretched in between the KH vortices at an angle of 45° with the horizontal plane. These streamwise vortices are not only identified in the early transitional stage of the mixing layer but also in self-similar turbulence conditions. The Re dependence of the "diameter" of these vortices is analyzed. Results obtained in spatial growing mixing layers show some evidences of pairing of secondary vortices; after a pairing of the primary Kelvin-Helmholtz (KH) vortices, the streamwise vortices are less numerous and their diameter has increased than before the pairing of KH vortices.

Nonlinear dynamics and control of multibody systems

The dynamics of flexible systems, such as robot manipulators , mechanical chains or multibody systems in general, is becoming increasingly important in engineering. This article deals with some nonlinearities that arise in the study of dynamics and control of multibody systems in connection to large rotations. Specifically, a numerical scheme that adresses the conservation of fundamental constants is presented in order to analyse the control-structure interaction problems.

On the numerical integration of rigid body nonlinear dynamics in presence of parameters singularities

One of the main complexities in the simulation of the nonlinear dynamics of rigid bodies consists in describing properly the finite rotations that they may undergo. It is well known that, to avoid singularities in the representation of the SO(3) rotation group, at least four parameters must be used. However, it is computationally expensive to use a four-parameters representation since, as only three of the parameters are independent, one needs to introduce constraint equations in the model, leading to differential-algebraic equations instead of ordinary differential ones. Three-parameter representations are numerically more efficient. Therefore, the objective of this paper is to evaluate numerically the influence of the parametrization and its singularities on the simulation of the dynamics of a rigid body. This is done through the analysis of a heavy top with a fixed point, using two three-parameter systems, Euler's angles and rotation vector. Theoretical results were used to guide the numerical simulation and to assure that all possible cases were analyzed. The two parametrizations were compared using several integrators. The results show that Euler's angles lead to faster integration compared to the rotation vector. An Euler's angles singular case, where representation approaches a theoretical singular point, was analyzed in detail. It is shown that on the contrary of what may be expected, 1) the numerical integration is very efficient, even more than for any other case, and 2) in spite of the uncertainty on the Euler's angles themselves, the body motion is well represented.