Records of sedimentary dynamics in the continental shelf and upper slope between Aveiro-Espinho (N Portugal)

The sedimentary unconsolidated cover of the Aveiro-Espinho continental shelf and upper slope (NW Portugal) records a complex interplay of processes including wave energy and currents, fluvial input, sediment transport alongshore and cross-shelf, geological and oceanographic processes and sediment sources and sinks. In order to study this record, a set of surface sediment samples was studied. Sediment grain size and composition, as well as the mineralogical composition (by XRD) of the fine (<63 mu m) and clay (<2 mu m) fractions and benthic microfaunal (foraminifera) data were analysed. Cluster analysis applied to the sedimentological data (grain size, sediment composition and mineralogy) allowed the establishment of three main zones corresponding to the: inner-, mid- and outer-shelf/upper slope. On the inner-shelf, the sedimentary coverture is composed of siliciclastic fine to very fine sand, essentially comprising modern (immature) terrigenous particles. The sediment grain size, as well as mineralogical and microfaunal composition, denote the high energetic conditions of this sector in which the alongshore transport of sand is predominantly southward and occurs mostly during the spring-summer oceanographic regime, when the main river providing sediments to this area, the River Douro, undergoes periods of drought. This effect may emphasize the erosive character of this coastal sector at present, since the Ria de Aveiro provides the shelf with few sediments. On the mid-shelf, an alongshore siliciclastic band of coarse sand and gravel can be found between the 40 m and 60 m isobaths. This gravelly deposit includes relic sediments deposited during lower sea-level stands. This structure stays on the surface due to the high bottom energy, which promotes the remobilization of the fine-grained sediments, and/or events of sediments bypassing. Benthic foraminifera density and "Benthic Foraminifera High Productivity" (BFHP) proxy values are in general low, which is consistent with the overall small supply of organic matter to the oceanic bottom in the inner- and mid-shelf. However, the Ria de Aveiro outflow, which delivers organic matter to the shelf, leaves its imprint mainly on the mid-shelf, identifiable by the increase in foraminifera density and BFHP values in front of the lagoon mouth. The higher values of BFHP along the 100 m isobath trace the present position of an oceanic thermal front whose situation may have changed in the last 3/5 ka BP. This zone marks a clear difference in the density, diversity and composition of benthic foraminifera assemblages. Here, in addition, sediment composition changes significantly, giving rise to carbonate-rich fine to medium sand in the deeper sector. The low bottom energy and the small sedimentation rate of the outer-shelf contributed to the preservation of a discontinuous carbonate-rich gravel band, between the 100 m and 140 m isobaths, also related to paleo-littorals, following the transgression that has occurred since the Last Glacial Maximum. The winter oceanographic regime favours the transport of fine grained sediments to the outer-shelf and upper slope. The inner- and mid-shelf, however, have low amounts of this kind of sediment and the Cretacic carbonated complexes Pontal da Galega and Pontal da Cartola, rocky outcrops located at the mid- and outer-shelf, act as morphological barriers to the cross-shelf transport of sediments. Thus a reduced sedimentation rate occurs in these deeper sectors, as indicated by the lower abundance of detrital minerals, which is compensated for the high sedimentary content of biogenic carbonates. The relatively high BFHP and Shannon Index values indicate water column stratification, high supply of organic matter and environmental stability, which provide favourable conditions for a diversified benthic fauna to flourish. These conditions also encourage authigenic chemical changes, favourable to glauconite formation, as well as illite and kaolinite degradation. Benthic foraminifera and clay mineral assemblages also reveal the effect of the internal waves pushing upward, and downslope losses of the sediments on the outer-shelf and upper slope.

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

Study of surface free energy on a lattice-gas model applied to Liquid bridge

The pulmonary crackling and the formation of liquid bridges are problems that for centuries have been attracting the attention of scientists. In order to study these phenomena, it was developed a canonical cubic lattice-gas­ like model to explain the rupture of liquid bridges in lung airways [A. Alencar et al., 2006, PRE]. Here, we further develop this model and add entropy analysis to study thermodynamic properties, such as free energy and force. The simulations were performed using the Monte Carlo method with Metropolis algorithm. The exchange between gas and liquid particles were performed randomly according to the Kawasaki dynamics and weighted by the Boltzmann factor. Each particle, which can be solid (s), liquid (l) or gas (g), has 26 neighbors: 6 + 12 + 8, with distances 1, √2 and √3, respectively. The energy of a lattice's site m is calculated by the following expression: Em = ∑k=126 Ji(m)j(k) in witch (i, j) = g, l or s. Specifically, it was studied the surface free energy of the liquid bridge, trapped between two planes, when its height is changed. For that, was considered two methods. First, just the internal energy was calculated. Then was considered the entropy. It was fond no difference in the surface free energy between this two methods. We calculate the liquid bridge force between the two planes using the numerical surface free energy. This force is strong for small height, and decreases as the distance between the two planes, height, is increased. The liquid-gas system was also characterized studying the variation of internal energy and heat capacity with the temperature. For that, was performed simulation with the same proportion of liquid and gas particle, but different lattice size. The scale of the liquid-gas system was also studied, for low temperature, using different values to the interaction Jij.

Dynamics of photoexcited carriers in the presence of disorder in radial heterostructured GaAs/AlGaAs/GaAs nanowires

Time-resolved photoluminescence was employed to study electron-hole dynamics in radial heterostructured GaAs/AlGaAs/GaAs core/inner shell/outer shell nanowires. It was found that impurity random potential results in a red shift of the recombination time maximum with respect to the photoluminescence peak energy.