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).

Signatures of quantum phase transitions in parallel quantum dots: Crossover from local moment to underscreened spin-1 Kondo physics

We study a strongly interacting "quantum dot 1" and a weakly interacting "dot 2" connected in parallel to metallic leads. Gate voltages can drive the system between Kondo-quenched and non-Kondo free-moment phases separated by Kosterlitz-Thouless quantum phase transitions. Away from the immediate vicinity of the quantum phase transitions, the physical properties retain signatures of first-order transitions found previously to arise when dot 2 is strictly noninteracting. As interactions in dot 2 become stronger relative to the dot-lead coupling, the free moment in the non-Kondo phase evolves smoothly from an isolated spin-one-half in dot 1 to a many-body doublet arising from the incomplete Kondo compensation by the leads of a combined dot spin-one. These limits, which feature very different spin correlations between dot and lead electrons, can be distinguished by weak-bias conductance measurements performed at finite temperatures.

Accurate ab initio potential energy surfaces for the (3)A '' and (3)A ' electronic states of the O(P-3) plus HBr system

In this work, we report the construction of potential energy surfaces for the (3)A '' and (3)A' states of the system O(P-3) + HBr. These surfaces are based on extensive ab initio calculations employing the MRCI+Q/CBS+SO level of theory. The complete basis set energies were estimated from extrapolation of MRCI+Q/aug-cc-VnZ(-PP) (n = Q, 5) results and corrections due to spin-orbit effects obtained at the CASSCF/aug-cc-pVTZ(-PP) level of theory. These energies, calculated over a region of the configuration space relevant to the study of the reaction O(P-3) + HBr -> OH + Br, were used to generate functions based on the many-body expansion. The three-body potentials were interpolated using the reproducing kernel Hilbert space method. The resulting surface for the (3)A '' electronic state contains van der Waals minima on the entrance and exit channels and a transition state 6.55 kcal/mol higher than the reactants. This barrier height was then scaled to reproduce the value of 5.01 kcal/mol, which was estimated from coupled cluster benchmark calculations performed to include high-order and core-valence correlation, as well as scalar relativistic effects. The (3)A' surface was also scaled, based on the fact that in the collinear saddle point geometry these two electronic states are degenerate. The vibrationally adiabatic barrier heights are 3.44 kcal/mol for the (3)A '' and 4.16 kcal/mol for the (3)A' state. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4705428]

Theoretical luminescence spectra in p-type superlattices based on InGaAsN

In this work, we present a theoretical photoluminescence (PL) for p-doped GaAs/InGaAsN nanostructures arrays. We apply a self-consistent method in the framework of the effective mass theory. Solving a full 8 x 8 Kane's Hamiltonian, generalized to treat different materials in conjunction with the Poisson equation, we calculate the optical properties of these systems. The trends in the calculated PL spectra, due to many-body effects within the quasi-two-dimensional hole gas, are analyzed as a function of the acceptor doping concentration and the well width. Effects of temperature in the PL spectra are also investigated. This is the first attempt to show theoretical luminescence spectra for GaAs/InGaAsN nanostructures and can be used as a guide for the design of nanostructured devices such as optoelectronic devices, solar cells, and others.

Structural properties and energetics of diffuse Rb-87 clusters in three-dimension

A correlated two-body basis function is used to describe the three-dimensional bosonic clusters interacting via two-body van der Waals potential. We calculate the ground state and the zero orbital angular momentum excited states for Rb-N clusters with up to N = 40. We solve the many-particle Schrodinger equation by potential harmonics expansion method, which keeps all possible two-body correlations in the calculation and determines the lowest effective many-body potential. We study energetics and structural properties for such diffuse clusters both at dimer and tuned scattering length. The motivation of the present study is to investigate the possibility of formation of N-body clusters interacting through the van der Waals interaction. We also compare the system with the well studied He, Ne, and Ar clusters. We also calculate correlation properties and observe the generalised Tjon line for large cluster. We test the validity of the shape-independent potential in the calculation of the ground state energy of such diffuse cluster. These are the first such calculations reported for Rb clusters. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4730972]

Junction Field Effect on the Retention Time for One-Transistor Floating-Body RAM

One-transistor floating-body random access memory retention time distribution is investigated on silicon-on-insulator UTBOX devices. It is shown that the average retention time can be improved by two to three orders of magnitude by reducing the body-junction electric field. However, the retention time distribution, which is mainly caused by the generation-recombination center density variation, remains similar.

New Isomers in the Full Seniority Scheme of Neutron-Rich Lead Isotopes: The Role of Effective Three-Body Forces

The neutron-rich lead isotopes, up to Pb-216, have been studied for the first time, exploiting the fragmentation of a primary uranium beam at the FRS-RISING setup at GSI. The observed isomeric states exhibit electromagnetic transition strengths which deviate from state-of-the-art shell-model calculations. It is shown that their complete description demands the introduction of effective three-body interactions and two-body transition operators in the conventional neutron valence space beyond Pb-208.

Modelling Mathematical Reasoning in Physics Education

Many findings from research as well as reports from teachers describe students' problem solving strategies as manipulation of formulas by rote. The resulting dissatisfaction with quantitative physical textbook problems seems to influence the attitude towards the role of mathematics in physics education in general. Mathematics is often seen as a tool for calculation which hinders a conceptual understanding of physical principles. However, the role of mathematics cannot be reduced to this technical aspect. Hence, instead of putting mathematics away we delve into the nature of physical science to reveal the strong conceptual relationship between mathematics and physics. Moreover, we suggest that, for both prospective teaching and further research, a focus on deeply exploring such interdependency can significantly improve the understanding of physics. To provide a suitable basis, we develop a new model which can be used for analysing different levels of mathematical reasoning within physics. It is also a guideline for shifting the attention from technical to structural mathematical skills while teaching physics. We demonstrate its applicability for analysing physical-mathematical reasoning processes with an example.

Body weight, metabolism and clock genes

Biological rhythms are present in the lives of almost all organisms ranging from plants to more evolved creatures. These oscillations allow the anticipation of many physiological and behavioral mechanisms thus enabling coordination of rhythms in a timely manner, adaption to environmental changes and more efficient organization of the cellular processes responsible for survival of both the individual and the species. Many components of energy homeostasis exhibit circadian rhythms, which are regulated by central (suprachiasmatic nucleus) and peripheral (located in other tissues) circadian clocks. Adipocyte plays an important role in the regulation of energy homeostasis, the signaling of satiety and cellular differentiation and proliferation. Also, the adipocyte circadian clock is probably involved in the control of many of these functions. Thus, circadian clocks are implicated in the control of energy balance, feeding behavior and consequently in the regulation of body weight. In this regard, alterations in clock genes and rhythms can interfere with the complex mechanism of metabolic and hormonal anticipation, contributing to multifactorial diseases such as obesity and diabetes. The aim of this review was to define circadian clocks by describing their functioning and role in the whole body and in adipocyte metabolism, as well as their influence on body weight control and the development of obesity.