Applying many-body expansion aiming to make ab in initio molecular dynamics more efficient

In this present work we present a methodology that aims to apply the many-body expansion to decrease the computational cost of ab initio molecular dynamics, keeping acceptable accuracy on the results. We implemented this methodology in a program which we called ManBo. In the many-body expansion approach, we partitioned the total energy E of the system in contributions of one body, two bodies, three bodies, etc., until the contribution of the Nth body [1-3]: E = E1 + E2 + E3 + …EN. The E1 term is the sum of the internal energy of the molecules; the term E2 is the energy due to interaction between all pairs of molecules; E3 is the energy due to interaction between all trios of molecules; and so on. In Manbo we chose to truncate the expansion in the contribution of two or three bodies, both for the calculation of the energy and for the calculation of the atomic forces. In order to partially include the many-body interactions neglected when we truncate the expansion, we can include an electrostatic embedding in the electronic structure calculations, instead of considering the monomers, pairs and trios as isolated molecules in space. In simulations we made we chose to simulate water molecules, and use the Gaussian 09 as external program to calculate the atomic forces and energy of the system, as well as reference program for analyzing the accuracy of the results obtained with the ManBo. The results show that the use of the many-body expansion seems to be an interesting approach for reducing the still prohibitive computational cost of ab initio molecular dynamics. The errors introduced on atomic forces in applying such methodology are very small. The inclusion of an embedding electrostatic seems to be a good solution for improving the results with only a small increase in simulation time. As we increase the level of calculation, the simulation time of ManBo tends to largely decrease in relation to a conventional BOMD simulation of Gaussian, due to better scalability of the methodology presented. References [1] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 3, 46 (2007). [2] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 4, 1 (2008). [3] R. Rivelino, P. Chaudhuri and S. Canuto; J. Chem. Phys., 118, 10593 (2003).

Annealed Ising model on hierarchical structures: a transfer matrix approach.

Spin systems in the presence of disorder are described by two sets of degrees of freedom, associated with orientational (spin) and disorder variables, which may be characterized by two distinct relaxation times. Disordered spin models have been mostly investigated in the quenched regime, which is the usual situation in solid state physics, and in which the relaxation time of the disorder variables is much larger than the typical measurement times. In this quenched regime, disorder variables are fixed, and only the orientational variables are duly thermalized. Recent studies in the context of lattice statistical models for the phase diagrams of nematic liquid-crystalline systems have stimulated the interest of going beyond the quenched regime. The phase diagrams predicted by these calculations for a simple Maier-Saupe model turn out to be qualitative different from the quenched case if the two sets of degrees of freedom are allowed to reach thermal equilibrium during the experimental time, which is known as the fully annealed regime. In this work, we develop a transfer matrix formalism to investigate annealed disordered Ising models on two hierarchical structures, the diamond hierarchical lattice (DHL) and the Apollonian network (AN). The calculations follow the same steps used for the analysis of simple uniform systems, which amounts to deriving proper recurrence maps for the thermodynamic and magnetic variables in terms of the generations of the construction of the hierarchical structures. In this context, we may consider different kinds of disorder, and different types of ferromagnetic and anti-ferromagnetic interactions. In the present work, we analyze the effects of dilution, which are produced by the removal of some magnetic ions. The system is treated in a “grand canonical" ensemble. The introduction of two extra fields, related to the concentration of two different types of particles, leads to higher-rank transfer matrices as compared with the formalism for the usual uniform models. Preliminary calculations on a DHL indicate that there is a phase transition for a wide range of dilution concentrations. Ising spin systems on the AN are known to be ferromagnetically ordered at all temperatures; in the presence of dilution, however, there are indications of a disordered (paramagnetic) phase at low concentrations of magnetic ions.

Anomalous diffusion evidenced by particle tracking at high frame rates.

Diffusion is a common phenomenon in nature and generally is associated with a system trying to reach a local or a global equilibrium state, as a result of highly irregular individual particle motion. Therefore it is of fundamental importance in physics, chemistry and biology. Particle tracking in complex fluids can reveal important characteristics of its properties. In living cells, we coat the microbead with a peptide (RGD) that binds to integrin receptors at the plasma membrane, which connects to the CSK. This procedure is based on the hypothesis that the microsphere can move only if the structure where it is attached move as well. Then, the observed trajectory of microbeads is a probe of the cytoskeleton (CSK), which is governed by several factors, including thermal diffusion, pressure gradients, and molecular motors. The possibility of separating the trajectories into passive and active diffusion may give information about the viscoelasticity of the cell structure and molecular motors activity. And also we could analyze the motion via generalized Stokes-Einstein relation, avoiding the use of any active techniques. Usually a 12 to 16 Frames Per Second (FPS) system is used to track the microbeads in cell for about 5 minutes. Several factors make this FPS limitation: camera computer communication, light, computer speed for online analysis among others. Here we used a high quality camera and our own software, developed in C++ and Linux, to reach high FPS. Measurements were conducted with samples for 10£ and 20£ objectives. We performed sequentially images with different intervals, all with 2 ¹s exposure. The sequences of intervals are in milliseconds: 4 5 ms (maximum speed) 14, 25, 50 and 100 FPS. Our preliminary results highlight the difference between passive and active diffusion, since the passive diffusion is represented by a Gaussian in the distribution of displacements of the center of mass of individual beads between consecutive frames. However, the active process, or anomalous diffusion, shows as long tails in the distribution of displacements.

Comparison of physical properties of two exogenous surfactants: new parameter.

The pulmonary surfactant has essential physical properties for normal lung function. The most important property is the surface tension. In this work, it was evaluated the surface tension of two commercial exogenous surfactants used in surfactant replacement therapy, poractant alfa (Curosurf, Chiesi Farmaceuticals, Italy) and beractant (Survanta, Abbott Laboratories, USA) using new parameters. A Langmuir film balance (Minitrough, KSV Instruments, Finland) was used to measure surface tension of poractant alfa and beractant samples. For both samples, we prepared a solution of 1 mg/m dissolved in chloroform (100π`), which was applied over a subphase of milli-Q water (175 ml) in the chamber of the balance. The chamber has two moving barriers that can change its surface area between a maximal value of 112.5 cm 2 , and a minimal value of 22.5 cm 2, defining a balance cycle. Each surfactant had its surface tension evaluated during 20 balance cycles for three times. Four quantities were calculated from the experiment: Minimum Surface Tension (MTS), defined as the surface tension at minimal surface area during the first cycle; Mean Work Cycle (MWC), defined as the mean hysteresis area of the measured surface tension curve of the last 16 balance cycles; Critical Active Surface Area in Compression (CASAC) or in Expansion (CASAE), defined as the maximal chamber area where the surfactant is active on the surface in compression or expansion. The t-test was applied to verify for statistical significance of the results. Comproved with the MST is the same reported in literature, the differences between MWC, CASAC, and CASAE were statistically significant (p<0.001). The MWC, CASAC and CASAE were higher for poractant alfa than for beractant. A higher MWC for poractant alfa means higher elastic recoil of the lung in comparison with beractant. Using a different methodology, our results showed that poractant alfa is probably more effective in a surfactant replacement therapy than beractant due the use of poractant alfa in relation to the use of beractant in preterm infants with Respiratory Distress Syndrome (RDS).

Fourier´s law from a chain of coupled anharmonic oscillators under energy conserving shot noise.

We analyze the transport of heat along a chain of particles interacting through anharmonic potentials consisting of quartic terms in addition to harmonic quadratic terms and subject to heat reservoirs at its ends. Each particle is also subject to an impulsive shot noise with exponentially distributed waiting times whose effect is to change the sign of its velocity, thus conserving the energy of the chain. We show that the introduction of this energy conserving stochastic noise leads to Fourier's law. That is for large system size L the heat current J behaves as J ‘approximately’ 1/L, which amounts to say that the conductivity k is constant. The conductivity is related to the current by J = kΔT/L, where ΔT is the difference in the temperatures of the reservoirs. The behavior of heat conductivity k for small intensities¸ of the shot noise and large system sizes L are obtained by assuming a scaling behavior of the type k = ‘L POT a Psi’(L’lambda POT a/b’) where a and b are scaling exponents. For the pure harmonic case a = b = 1, characterizing a ballistic conduction of heat when the shot noise is absent. For the anharmonic case we found values for the exponents a and b smaller then 1 and thus consistent with a superdiffusive conduction of heat without the shot noise. We also show that the heat conductivity is not constant but is an increasing function of temperature.

Heat conduction in a chain of harmonic and anharmonic oscillators under the presence of an energy conserving noise

Reproducing Fourier's law of heat conduction from a microscopic stochastic model is a long standing challenge in statistical physics. As was shown by Rieder, Lebowitz and Lieb many years ago, a chain of harmonically coupled oscillators connected to two heat baths at different temperatures does not reproduce the diffusive behaviour of Fourier's law, but instead a ballistic one with an infinite thermal conductivity. Since then, there has been a substantial effort from the scientific community in identifying the key mechanism necessary to reproduce such diffusivity, which usually revolved around anharmonicity and the effect of impurities. Recently, it was shown by Dhar, Venkateshan and Lebowitz that Fourier's law can be recovered by introducing an energy conserving noise, whose role is to simulate the elastic collisions between the atoms and other microscopic degrees of freedom, which one would expect to be present in a real solid. For a one-dimensional chain this is accomplished numerically by randomly flipping - under the framework of a Poisson process with a variable “rate of collisions" - the sign of the velocity of an oscillator. In this poster we present Langevin simulations of a one-dimensional chain of oscillators coupled to two heat baths at different temperatures. We consider both harmonic and anharmonic (quartic) interactions, which are studied with and without the energy conserving noise. With these results we are able to map in detail how the heat conductivity k is influenced by both anharmonicity and the energy conserving noise. We also present a detailed analysis of the behaviour of k as a function of the size of the system and the rate of collisions, which includes a finite-size scaling method that enables us to extract the relevant critical exponents. Finally, we show that for harmonic chains, k is independent of temperature, both with and without the noise. Conversely, for anharmonic chains we find that k increases roughly linearly with the temperature of a given reservoir, while keeping the temperature difference fixed.

Model for the Dynamics of the Cytoskeleton.

Particle tracking of microbeads attached to the cytoskeleton (CSK) reveals an intermittent dynamic. The mean squared displacement (MSD) is subdiffusive for small Δt and superdiffusive for large Δt, which are associated with periods of traps and periods of jumps respectively. The analysis of the displacements has shown a non-Gaussian behavior, what is indicative of an active motion, classifying the cells as a far from equilibrium material. Using Langevin dynamics, we reconstruct the dynamic of the CSK. The model is based on the bundles of actin filaments that link themself with the bead RGD coating, trapping it in an harmonic potential. We consider a one- dimensional motion of a particle, neglecting inertial effects (over-damped Langevin dynamics). The resultant force is decomposed in friction force, elastic force and random force, which is used as white noise representing the effect due to molecular agitation. These description until now shows a static situation where the bead performed a random walk in an elastic potential. In order to modeling the active remodeling of the CSK, we vary the equilibrium position of the potential. Inserting a motion in the well center, we change the equilibrium position linearly with time with constant velocity. The result found exhibits a MSD versus time ’tau’ with three regimes. The first regime is when ‘tau’ < ‘tau IND 0’, where ‘tau IND 0’ is the relaxation time, representing the thermal motion. At this regime the particle can diffuse freely. The second regime is a plateau, ‘tau IND 0’ < ‘tau’ < ‘tau IND 1’, representing the particle caged in the potential. Here, ‘tau IND 1’ is a characteristic time that limit the confinement period. And the third regime, ‘tau’ > ‘tau IND 1’, is when the particles are in the superdiffusive behavior. This is where most of the experiments are performed, under 20 frames per second (FPS), thus there is no experimental evidence that support the first regime. We are currently performing experiments with high frequency, up to 100 FPS, attempting to visualize this diffusive behavior. Beside the first regime, our simple model can reproduce MSD curves similar to what has been found experimentally, which can be helpful to understanding CSK structure and properties.

BDPI: Biblioteca Digital da Produção Intelectual da USP

Biblioteca Digital da Produção Intelectual da Universidade de São Paulo Objetivo: captura, armazenamento, indexação, preservação, busca e disseminação de material digital.

Infraestrutura para as aplicações corporativa em implantação

BDPI ( Biblioteca Digital da Produção Intelectual). Portal de Revistas USP. Portal de WorkShops de Publicação Científica. Portal de Obras Raras. Portal de Busca Integrada. Dédalus (Sistema Aleph).

LOCKSS - USP: uma iniciativa de Preservação Digital

LOCKSS e preservação digital no BRASIL. Podemos definir Preservação Digital como o conjunto de atividades ou processos responsáveis por garantir o acesso contínuo a longo-prazo à informação científica, tecnológica e cultural existente em formatos digitais. O projeto LOCKSS (Lot Of Copies Keep Stuff Safe) é um software de código aberto desenvolvido pela Universidade de Stanford, com o intuito de preservar dados digitais. Possui uma abordagem peer-to-peer e garante a integridade e disponibilidade do dos acervos digitais preservados. No Brasil, a Rede Brasileira de Serviços de Preservação Digital, ou Rede Cariniana, é uma iniciativa do Instituto Brasileiro de Informação em Ciência e Tecnologia para a salvaguarda do patrimônio científico, tecnológico e cultural, de forma distribuída, por instituições de pesquisa em parceria com o IBICT. Sua implantação foi fundamentada na utilização do LOCKSS em uma infraestrutura descentralizada, utilizando recursos de computação distribuída .

Portal de revistas da USP

O Portal de Revistas da USP apresenta-se como uma estratégia do Sistema Integrado de Bibliotecas, por meio de seu Programa de Apoio às Publicações Científicas Periódicas, para aumentar a visibilidade e acessibilidade das revistas científicas da Universidade de São Paulo. Para fazer parte deste Programa as revistas devem cumprir com critérios de qualidade, definidos por uma Comissão de Credenciamento, formada por professores e bibliotecários, e pautados em padrões internacionais. Atualmente 62 publicações estão credenciadas neste Programa.

Produção e preservação de conteúdos digitais: montagem, operação e gestão de laboratórios de digitalização

Com o objetivo de preservar e disseminar o acervo da Universidade de São Paulo, foi criada a Biblioteca de Obras Raras e Especiais. Contamos atualmente com cerca de 2500 obras digitalizadas em parceria com a Brasiliana USP e apoio da FAPESP. Estas obras estão disponibilizadas online, para consulta aberta pelo amplo público, promovendo assim acesso ao conhecimento. Utilizamos para isto o Corisco, uma plataforma desenvolvida por pesquisadores da Brasiliana USP baseada em DSpace, software de código aberto e com recursos que promovem a interoperabilidade com outras bases. O DSpace é atualmente o aplicativo de Repositório Digital de código aberto mais utilizado nas Universidades ao redor do mundo. As capturas originais em alta resolução das obras são armazenadas e indexadas para fins de preservação. Em outra frente, realiza-se o tratamento, compactação e reconhecimento ótico de caracteres nas imagens, para que então as obras possam ser disponibilizadas na internet com o recurso de busca no conteúdo. A consulta é livre e gratuita, garantindo assim, além da preservação, o acesso ao conhecimento pela comunidade acadêmica e público em geral.

Estratégias para implementar o acesso aberto na Universidade de São Paulo

Objetivo: Compartilhar a experiência da construção da política e do repositório institucional da USP, a partir da descrição das distintas e complexas estratégias desenvolvidas para a mudança de paradigma e inserção da universidade no cenário atual do acesso aberto.

Digitalização e disponibilização de obras raras e especiais

Com o objetivo de preservar e disseminar o acervo da Universidade de São Paulo, foi disponibilizada a Biblioteca Digital de Obras Raras e Especiais, com apoio da FAPESP, utilizando a plataforma Corisco, desenvolvida por pesquisadores da Brasiliana USP, baseada em softwares de código aberto com recursos que permitem a interoperabilidade com outros bancos de dados. A Plataforma Corisco é um produto desenvolvido pela Brasiliana USP, consistindo de um sistema integrado de aplicativos de código aberto e recomendações para sustentar a implantação e gerenciamento de repositórios digitais de acordo com padrões de interoperabilidade. A Plataforma Corisco roda sobre o DSpace, um dos aplicativos de Repositório Digital de código aberto mais utilizados nas Universidades ao redor do mundo. As obras, em domínio público, são digitalizadas e suas imagens tratadas, para então serem disponibilizadas. E graças ao reconhecimento ótico de caracteres é possível realizar a busca no conteúdo, potencializando resultados. A consulta é livre e gratuita, garantindo assim, além da preservação, o acesso ao conhecimento pela comunidade acadêmica e público em geral.