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

Records of migration in the exoplanet configurations

When compared to our Solar System, many exoplanet systems exhibit quite unusual planet configurations; some of these are hot Jupiters, which orbit their central stars with periods of a few days, others are resonant systems composed of two or more planets with commensurable orbital periods. It has been suggested that these configurations can be the result of a migration processes originated by tidal interactions of the planets with disks and central stars. The process known as planet migration occurs due to dissipative forces which affect the planetary semi-major axes and cause the planets to move towards to, or away from, the central star. In this talk, we present possible signatures of planet migration in the distribution of the hot Jupiters and resonant exoplanet pairs. For this task, we develop a semi-analytical model to describe the evolution of the migrating planetary pair, based on the fundamental concepts of conservative and dissipative dynamics of the three-body problem. Our approach is based on an analysis of the energy and the orbital angular momentum exchange between the two-planet system and an external medium; thus no specific kind of dissipative forces needs to be invoked. We show that, under assumption that dissipation is weak and slow, the evolutionary routes of the migrating planets are traced by the stationary solutions of the conservative problem (Birkhoff, Dynamical systems, 1966). The ultimate convergence and the evolution of the system along one of these modes of motion are determined uniquely by the condition that the dissipation rate is sufficiently smaller than the roper frequencies of the system. We show that it is possible to reassemble the starting configurations and migration history of the systems on the basis of their final states, and consequently to constrain the parameters of the physical processes involved.

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.

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.

Multi-planet extrasolar systems - detection and dynamics

20 years after the discovery of the first planets outside our solar system, the current exoplanetary population includes more than 700 confirmed planets around main sequence stars. Approximately 50% belong to multiple-planet systems in very diverse dynamical configurations, from two-planet hierarchical systems to multiple resonances that could only have been attained as the consequence of a smooth large-scale orbital migration. The first part of this paper reviews the main detection techniques employed for the detection and orbital characterization of multiple-planet systems, from the (now) classical radial velocity (RV) method to the use of transit time variations (TTV) for the identification of additional planetary bodies orbiting the same star. In the second part we discuss the dynamical evolution of multi-planet systems due to their mutual gravitational interactions. We analyze possible modes of motion for hierarchical, secular or resonant configurations, and what stability criteria can be defined in each case. In some cases, the dynamics can be well approximated by simple analytical expressions for the Hamiltonian function, while other configurations can only be studied with semi-analytical or numerical tools. In particular, we show how mean-motion resonances can generate complex structures in the phase space where different libration islands and circulation domains are separated by chaotic layers. In all cases we use real exoplanetary systems as working examples.

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]

Stability and dynamics of terrestrial planets in binary star systems: application on Alpha Centauri

Binary stars are frequent in the universe, with about 50% of the known main sequence stars being located at a multiple star system (Abt, 1979). Even though, they are universally thought as second rate sites for the location of exo-planets and the habitable zone, due to the difficulties of detection and high perturbation that could prevent planet formation and long term stability. In this work we show that planets in binary star systems can have regular orbits and remain on the habitable zone. We introduce a stability criterium based on the solution of the restricted three body problem and apply it to describe the short period planar and three-dimentional stability zones of S-type orbits around each star of the Alpha Centauri system. We develop as well a semi-analytical secular model to study the long term dynamics of fictional planets in the habitable zone of those stars and we verify that planets on the habitable zone would be in regular orbits with any eccentricity and with inclination to the binary orbital plane up until 35 degrees. We show as well that the short period oscillations on the semi-major axis is 100 times greater than the Earth's, but at all the time the planet would still be found inside the Habitable zone.

Means to an end: Neotropical parrots manage to pull strings to meet their goals

Although parrots share with corvids and primates many of the traits believed to be associated with advanced cognitive processing, knowledge of parrot cognition is still limited to a few species, none of which are Neotropical. Here we examine the ability of three Neotropical parrot species (Blue-Fronted Amazons, Hyacinth and Lear`s macaws) to spontaneously solve a novel physical problem: the string-pulling test. The ability to pull up a string to obtain out-of-reach food has been often considered a cognitively complex task, as it requires the use of a sequence of actions never previously assembled, along with the ability to continuously monitor string, food and certain body movements. We presented subjects with pulling tasks where we varied the spatial relationship between the strings, the presence of a reward and the physical contact between the string and reward to determine whether (1) string-pulling is goal-oriented in these parrots, (2) whether the string is recognized as a means to obtain the reward and (3) whether subjects can visually determine the continuity between the string and the reward, selecting only those strings for which no physical gaps between string and reward were present. Our results show that some individuals of all species were able to use the string as a means to reach a specific goal, in this case, the retrieval of the food treat. Also, subjects from both macaw species were able to visually determine the presence of physical continuity between the string and reward, making their choices consistently with the recognition that no gaps should be present between the string and the reward. Our findings highlight the potential of this taxonomic group for the understanding of the underpinnings of cognition in evolutionarily distant groups such as birds and primates.

On bifurcation of solutions of the Yamabe problem in product manifolds

We study local rigidity and multiplicity of constant scalar curvature metrics in arbitrary products of compact manifolds. Using (equivariant) bifurcation theory we determine the existence of infinitely many metrics that are accumulation points of pairwise non-homothetic solutions of the Yamabe problem. Using local rigidity and some compactness results for solutions of the Yamabe problem, we also exhibit new examples of conformal classes (with positive Yamabe constant) for which uniqueness holds. (C) 2011 Elsevier Masson SAS. All rights reserved.

Effects of excess body mass on strength and fatigability of quadriceps in postmenopausal women

Objective: Obesity is a major public health problem leading to, among other things, reduced functional capacity. Moreover, obesity-related declines in functional capacity may be compounded by the detrimental consequences of menopause. The aim of this study was to understand the potential effects of excess body mass on measures of functional capacity in postmenopausal women. Methods: Forty-five postmenopausal women aged 50 to 60 years were divided into two groups according to body mass index (BMI): obese (BMI, >= 30 kg/m(2); n = 19) and nonobese (BMI, 18.5-29.9 kg/m(2); n = 26). To determine clinical characteristics, body composition, bone mineral density, and maximal exercise testing was performed, and a 3-day dietary record was estimated. To assess quadriceps function, isokinetic exercise testing at 60 degrees per second (quadriceps strength) and at 300 degrees per second (quadriceps fatigue) was performed. Results: The absolute value of the peak torque was not significantly different between the groups; however, when the data were normalized by body mass and lean mass, significantly lower values were observed for obese women compared with those in the nonobese group (128% +/- 25% vs 155% +/- 24% and 224% +/- 38% vs 257% +/- 47%, P < 0.05). The fatigue index did not show any significant difference for either group; however, when the data were normalized by the body mass and lean mass, significantly lower values were observed for obese women (69% +/- 16% vs 93% +/- 18% and 120% +/- 25% vs. 135% +/- 23%, P < 0.01). Conclusions: Our results show that despite reduced muscle force, the combination of obesity and postmenopause may be associated with greater resistance to muscle fatigue.

Antecedents of the Intent to Recommend: a problem with fast food restaurants

The fast and strong social and economic transformations in the economies of many countries has raised the competition for consumers. One of the elements required to adapt to such scenario is knowing customers and their perceptions about products or services, mainly regarding word of mouth recommendations. This study adapts, to the fast food business, a model originally designed to analyze the antecedents of the intent to recommend by clients of formal restaurants. Three constructs were considered: service quality, satisfaction, and social well-being, the latter comprised of positive and negative affections. Six hypotheses were considered, three of which relating to social well-being (that it influences satisfaction, service quality, and the intent to recommend), two relating to service quality (that in influences the intent to recommend and satisfaction), and one relating to the influence of satisfaction on the intent to recommend. None was rejected, indicating adherence and adjustment of the simplication and adaptation of the consolidated model. Through a successful empirical application, the main contribution made by this research is the simplification of a model through its application in a similar context, but with a different scope.

Hypergraphs with many Kneser colorings

For fixed positive integers r, k and E with 1 <= l < r and an r-uniform hypergraph H, let kappa(H, k, l) denote the number of k-colorings of the set of hyperedges of H for which any two hyperedges in the same color class intersect in at least l elements. Consider the function KC(n, r, k, l) = max(H epsilon Hn) kappa(H, k, l), where the maximum runs over the family H-n of all r-uniform hypergraphs on n vertices. In this paper, we determine the asymptotic behavior of the function KC(n, r, k, l) for every fixed r, k and l and describe the extremal hypergraphs. This variant of a problem of Erdos and Rothschild, who considered edge colorings of graphs without a monochromatic triangle, is related to the Erdos-Ko-Rado Theorem (Erdos et al., 1961 [8]) on intersecting systems of sets. (C) 2011 Elsevier Ltd. All rights reserved.

Relationship between body size and flying-related structures in Neotropical social wasps (Polistinae, Vespidae, Hymenoptera)

Body size influences wing shape and associated muscles in flying animals which is a conspicuous phenomenon in insects, given their wide range in body size. Despite the significance of this, to date, no detailed study has been conducted across a group of species with similar biology allowing a look at specific relationship between body size and flying structures. Neotropical social vespids are a model group to study this problem as they are strong predators that rely heavily on flight while exhibiting a wide range in body size. In this paper we describe the variation in both wing shape, as wing planform, and mesosoma muscle size along the body size gradient of the Neotropical social wasps and discuss the potential factors affecting these changes. Analyses of 56 species were conducted using geometric morphometrics for the wings and lineal morphometrics for the body; independent contrast method regressions were used to correct for the phylogenetic effect. Smaller vespid species exhibit rounded wings, veins that are more concentrated in the proximal region, larger stigmata and the mesosoma is proportionally larger than in larger species. Meanwhile, larger species have more elongated wings, more distally extended venation, smaller stigmata and a proportionally smaller mesosoma. The differences in wing shape and other traits could be related to differences in flight demands caused by smaller and larger body sizes. Species around the extremes of body size distribution may invest more in flight muscle mass than species of intermediate sizes.

The need of a weight management control program in judo: a proposal based on the successful case of wrestling

Judo competitions are divided into weight classes. However, most athletes reduce their body weight in a few days before competition in order to obtain a competitive advantage over lighter opponents. To achieve fast weight reduction, athletes use a number of aggressive nutritional strategies so many of them place themselves at a high health-injury risk. In collegiate wrestling, a similar problem has been observed and three wrestlers died in 1997 due to rapid weight loss regimes. After these deaths, the National Collegiate Athletic Association had implemented a successful weight management program which was proven to improve weight management behavior. No similar program has ever been discussed by judo federations even though judo competitors present a comparable inappropriate pattern of weight control. In view of this, the basis for a weight control program is provided in this manuscript, as follows: competition should begin within 1 hour after weigh-in, at the latest; each athlete is allowed to be weighed-in only once; rapid weight loss as well as artificial rehydration (i.e., saline infusion) methods are prohibited during the entire competition day; athletes should pass the hydration test to get their weigh-in validated; an individual minimum competitive weight (male athletes competing at no less than 7% and females at no less than 12% of body fat) should be determined at the beginning of each season; athletes are not allowed to compete in any weight class that requires weight reductions greater than 1.5% of body weight per week. In parallel, educational programs should aim at increasing the athletes', coaches' and parents' awareness about the risks of aggressive nutritional strategies as well as healthier ways to properly manage body weight.