Proposal of three thermodynamic variables to discriminate between storms associated with hail and storms with intense rainfall in Catalonia (Proposta de tres variables termodinàmiques per discriminar entre tempestes amb pedregades i tempestes amb pluges intenses)

Between late spring and early fall, the development of storms is common in Catalonia. Despite the fact that they usually produce heavy showers of short duration, they can also involve severe weather with ice pellets or hail. While the latter usually affect inland regions, and there are numerous publications on these cases; the analysis of events affecting the coast and causing damage to public and private properties is not so well developed. The aim of this study is to provide additional thermodynamic indicators that help differentiate storms with hail from storms without hail, considering cases that have affected various regions of Catalonia, mainly coastal areas. The aim is to give more information to improve prognosis and the ability to detail information in these situations. The procedure developed involved the study of several episodes of heavy rainfall and hail that hit Catalonia during the 2003-2009 period, mainly in the province of Girona, and validated the proposal during the campaign of late summer and fall of 2009, as well as 2012. For each case, several variables related to temperature, humidity and wind were analyzed at different levels of the atmosphere, while the information provided by the radio sounding in Barcelona was also taken into account. From this study, it can be concluded that the temperature difference between 500 hPa and 850 hPa, the humidity in the lower layers of the atmosphere and the LI index are good indicators for the detection of storms with associated hail.

Microstructural evolution of primary crystallization in diffusion-controlled grain growth

A model has been developed for evaluating grain size distributions in primary crystallizations where the grain growth is diffusion controlled. The body of the model is grounded in a recently presented mean-field integration of the nucleation and growth kinetic equations, modified conveniently in order to take into account a radius-dependent growth rate, as occurs in diffusion-controlled growth. The classical diffusion theory is considered, and a modification of this is proposed to take into account interference of the diffusion profiles between neighbor grains. The potentiality of the mean-field model to give detailed information on the grain size distribution and transformed volume fraction for transformations driven by nucleation and either interface- or diffusion-controlled growth processes is demonstrated. The model is evaluated for the primary crystallization of an amorphous alloy, giving an excellent agreement with experimental data. Grain size distributions are computed, and their properties are discussed.

Fourier transform convolution integrals applied to generalized Born molecular volume.

Generalized Born methods are currently among the solvation models most commonly used for biological applications. We reformulate the generalized Born molecular volume method initially described by (Lee et al, 2003, J Phys Chem, 116, 10606; Lee et al, 2003, J Comp Chem, 24, 1348) using fast Fourier transform convolution integrals. Changes in the initial method are discussed and analyzed. Finally, the method is extensively checked with snapshots from common molecular modeling applications: binding free energy computations and docking. Biologically relevant test systems are chosen, including 855-36091 atoms. It is clearly demonstrated that, precision-wise, the proposed method performs as good as the original, and could better benefit from hardware accelerated boards.

Energy relaxation in nonlinear one-dimensional lattices

We study energy relaxation in thermalized one-dimensional nonlinear arrays of the Fermi-Pasta-Ulam type. The ends of the thermalized systems are placed in contact with a zero-temperature reservoir via damping forces. Harmonic arrays relax by sequential phonon decay into the cold reservoir, the lower-frequency modes relaxing first. The relaxation pathway for purely anharmonic arrays involves the degradation of higher-energy nonlinear modes into lower-energy ones. The lowest-energy modes are absorbed by the cold reservoir, but a small amount of energy is persistently left behind in the array in the form of almost stationary low-frequency localized modes. Arrays with interactions that contain both a harmonic and an anharmonic contribution exhibit behavior that involves the interplay of phonon modes and breather modes. At long times relaxation is extremely slow due to the spontaneous appearance and persistence of energetic high-frequency stationary breathers. Breather behavior is further ascertained by explicitly injecting a localized excitation into the thermalized arrays and observing the relaxation behavior.

Cluster glasses of ultrasoft particles

We present molecular dynamics (MD) simulations results for dense fluids of ultrasoft, fully penetrable particles. These are a binary mixture and a polydisperse system of particles interacting via the generalized exponential model, which is known to yield cluster crystal phases for the corresponding monodisperse systems. Because of the dispersity in the particle size, the systems investigated in this work do not crystallize and form disordered cluster phases. The clusteringtransition appears as a smooth crossover to a regime in which particles are mostly located in clusters, isolated particles being infrequent. The analysis of the internal cluster structure reveals microsegregation of the big and small particles, with a strong homo-coordination in the binary mixture. Upon further lowering the temperature below the clusteringtransition, the motion of the clusters" centers-of-mass slows down dramatically, giving way to a cluster glass transition. In the cluster glass, the diffusivities remain finite and display an activated temperature dependence, indicating that relaxation in the cluster glass occurs via particle hopping in a nearly arrested matrix of clusters. Finally we discuss the influence of the microscopic dynamics on the transport properties by comparing the MD results with Monte Carlo simulations.

Sugarcane crop efficiency in two growing seasons in São Paulo State, Brazil

A conceptual framework for crop production efficiency was derived using thermodynamic efficiency concept, in order to generate a tool for performance evaluation of agricultural systems and to quantify the interference of determining factors on this performance. In Thermodynamics, efficiency is the ratio between the output and input of energy. To establish this relationship in agricultural systems, it was assumed that the input energy is represented by the attainable crop yield, as predicted through simulation models based on environmental variables. The method of FAO's agroecological zones was applied to the assessment of the attainable sugarcane yield, while Instituto Brasileiro de Geografia e Estatística (IBGE) data were used as observed yield. Sugarcane efficiency production in São Paulo state was evaluated in two growing seasons, and its correlation with some physical factors that regulate production was calculated. A strong relationship was identified between crop production efficiency and soil aptitude. This allowed inferring the effect of agribusiness factors on crop production efficiency. The relationships between production efficiency and climatic variables were also quantified and indicated that solar radiation, annual rainfall, water deficiency, and maximum air temperature are the main factors affecting the sugarcane production efficiency.

Viral self-assembly as a thermodynamic process

The protein shells, or capsids, of nearly all spherelike viruses adopt icosahedral symmetry. In the present Letter, we propose a statistical thermodynamic model for viral self-assembly. We find that icosahedral symmetry is not expected for viral capsids constructed from structurally identical protein subunits and that this symmetry requires (at least) two internal switching configurations of the protein. Our results indicate that icosahedral symmetry is not a generic consequence of free energy minimization but requires optimization of internal structural parameters of the capsid proteins

Thermostatistical description of small systems in nonequilibrium conditions: Energy conversion and harvesting

Hysteresis cycles are very important features of energy conversion and harvesting devices, such as batteries. The efficiency of these may be strongly affected by the physical size of the system. Here, we show that in systems which are small enough, the existence of physical boundaries which produce nonhomogeneities of the interaction potential gives rise to inflections and barriers in the associated free energy. This in turn brings on irreversible processes which can be triggered under suitable external conditions imposed by a heat bath. As an example, by controlling the temperature, the state of a small system may be impelled to oscillate between two different structural configurations or aggregation states avoiding equilibrium coexistence and therefore dissipating energy. This cyclical behavior associated with a hysteresis cycle may be prototypical of energy conversion, storage, or generating nanodevices, as exemplified by Li-ion insertion batteries.

Entropic transport: kinetics, scaling, and control mechanisms

We show that transport in the presence of entropic barriers exhibits peculiar characteristics which makes it distinctly different from that occurring through energy barriers. The constrained dynamics yields a scaling regime for the particle current and the diffusion coefficient in terms of the ratio between the work done to the particles and available thermal energy. This interesting property, genuine to the entropic nature of the barriers, can be utilized to effectively control transport through quasi-one-dimensional structures in which irregularities or tortuosity of the boundaries cause entropic effects. The accuracy of the kinetic description has been corroborated by simulations. Applications to different dynamic situations involving entropic barriers are outlined.

Kinetic equations for difussion in the presence of entropic barriers

We use the mesoscopic nonequilibrium thermodynamics theory to derive the general kinetic equation of a system in the presence of potential barriers. The result is applied to a description of the evolution of systems whose dynamics is influenced by entropic barriers. We analyze in detail the case of diffusion in a domain of irregular geometry in which the presence of the boundaries induces an entropy barrier when approaching the exact dynamics by a coarsening of the description. The corresponding kinetic equation, named the Fick-Jacobs equation, is obtained, and its validity is generalized through the formulation of a scaling law for the diffusion coefficient which depends on the shape of the boundaries. The method we propose can be useful to analyze the dynamics of systems at the nanoscale where the presence of entropy barriers is a common feature.

A Thermokinetic approach to radiative heat transfer at the nanoscale

Radiative heat exchange at the nanoscale presents a challenge for several areas due to its scope and nature. Here, we provide a thermokinetic description of microscale radiative energy transfer including phonon-photon coupling manifested through a non-Debye relaxation behavior. We show that a lognormal-like distribution of modes of relaxation accounts for this non-Debye relaxation behavior leading to the thermal conductance. We also discuss the validity of the fluctuation-dissipation theorem. The general expression for the thermal conductance we obtain fits existing experimental results with remarkable accuracy. Accordingly, our approach offers an overall explanation of radiative energy transfer through micrometric gaps regardless of geometrical configurations and distances.

Algebraic decay of velocity fluctuations near a wall

Computer simulations of the dynamics of a colloidal particle suspended in a fluid confined by an interface show that the asymptotic decay of the velocity correlation functions is algebraic. The exponents of the long-time tails depend on the direction of motion of the particle relative to the surface, as well as on the specific nature of the boundary conditions. In particular, we find that for the angular velocity correlation function, the decay in the presence of a slip surface is faster than the one corresponding to a stick one. An intuitive picture is introduced to explain the various long-time tails, and the simulations are compared with theoretical expressions where available.

Entropic transport in confined media: a challenge for computational studies in biological and soft-matter systems

Transport in small-scale biological and soft-matter systems typically occurs under confinement conditions in which particles proceed through obstacles and irregularities of the boundaries that may significantly alter their trajectories. A transport model that assimilates the confinement to the presence of entropic barriers provides an efficient approach to quantify its effect on the particle current and the diffusion coefficient. We review the main peculiarities of entropic transport and treat two cases in which confinement effects play a crucial role, with the appearance of emergent properties. The presence of entropic barriers modifies the mean first-passage time distribution and therefore plays a very important role in ion transport through micro- and nano-channels. The functionality of molecular motors, modeled as Brownian ratchets, is strongly affected when the motor proceeds in a confined medium that may constitute another source of rectification. The interplay between ratchet and entropic rectification gives rise to a wide variety of dynamical behaviors, not observed when the Brownian motor proceeds in an unbounded medium. Entropic transport offers new venues of transport control and particle manipulation and new ways to engineer more efficient devices for transport at the nanoscale.

Binding free energy differences in a TCR-peptide-MHC complex induced by a peptide mutation: a simulation analysis.

Recognition by the T-cell receptor (TCR) of immunogenic peptides presented by class I major histocompatibility complexes (MHCs) is the determining event in the specific cellular immune response against virus-infected cells or tumor cells. It is of great interest, therefore, to elucidate the molecular principles upon which the selectivity of a TCR is based. These principles can in turn be used to design therapeutic approaches, such as peptide-based immunotherapies of cancer. In this study, free energy simulation methods are used to analyze the binding free energy difference of a particular TCR (A6) for a wild-type peptide (Tax) and a mutant peptide (Tax P6A), both presented in HLA A2. The computed free energy difference is 2.9 kcal/mol, in good agreement with the experimental value. This makes possible the use of the simulation results for obtaining an understanding of the origin of the free energy difference which was not available from the experimental results. A free energy component analysis makes possible the decomposition of the free energy difference between the binding of the wild-type and mutant peptide into its components. Of particular interest is the fact that better solvation of the mutant peptide when bound to the MHC molecule is an important contribution to the greater affinity of the TCR for the latter. The results make possible identification of the residues of the TCR which are important for the selectivity. This provides an understanding of the molecular principles that govern the recognition. The possibility of using free energy simulations in designing peptide derivatives for cancer immunotherapy is briefly discussed.

Atomic Ordering and Martensitic Transitions in Cu-Zn-Al Shape Memory Alloys

We present results from both, calorimetric and dilatometric studies of the isothermal ordering process taking place in a Cu-Zn-Al shape memory alloy after quenches from Tq temperatures ranging from 350 K to 1200 K. The dissipated energy and the length variations of the system are obtained during the process. The change of these quantities in the whole process have been compared with the difference [MATH] between Ms, measured after the relaxation and Ms measured just after the quench. We obtain that these three quantities present, as a function of Tq, the same qualitative behaviour. These changes are then associated with changes of the L21 ordering after the quench in the system. The relaxational process does not follow a single exponential decay. Instead, a continuous slowing down is observed. A relaxation time [MATH] has been defined to characterize the relaxation rate. We show that [MATH] depends on both the annealing and the quenching (Tq [MATH] 800 K) temperatures through an Arrhenius law.