Asymptotic analysis of combined breather-kink modes in a Fermi-Pasta-Ulam chain

We find approximations to travelling breather solutions of the one-dimensional Fermi-Pasta-Ulam (FPU) lattice. Both bright breather and dark breather solutions are found. We find that the existence of localised (bright) solutions depends upon the coefficients of cubic and quartic terms of the potential energy, generalising an earlier inequality derived by James [CR Acad Sci Paris 332, 581, (2001)]. We use the method of multiple scales to reduce the equations of motion for the lattice to a nonlinear Schr{\"o}dinger equation at leading order and hence construct an asymptotic form for the breather. We show that in the absence of a cubic potential energy term, the lattice supports combined breathing-kink waveforms. The amplitude of breathing-kinks can be arbitrarily small, as opposed to traditional monotone kinks, which have a nonzero minimum amplitude in such systems. We also present numerical simulations of the lattice, verifying the shape and velocity of the travelling waveforms, and confirming the long-lived nature of all such modes.

MOLECULAR MODELING OF EPON 862-DETDA / CARBON COMPOSITES

The thermoset epoxy resin EPON 862, coupled with the DETDA hardening agent, are utilized as the polymer matrix component in many graphite (carbon fiber) composites. Because it is difficult to experimentally characterize the interfacial region, computational molecular modeling is a necessary tool for understanding the influence of the interfacial molecular structure on bulk-level material properties. The purpose of this research is to investigate the many possible variables that may influence the interfacial structure and the effect they will have on the mechanical behavior of the bulk level composite. Molecular models are established for EPON 862-DETDA polymer in the presence of a graphite surface. Material characteristics such as polymer mass-density, residual stresses, and molecular potential energy are investigated near the polymer/fiber interface. Because the exact degree of crosslinking in these thermoset systems is not known, many different crosslink densities (degrees of curing) are investigated. It is determined that a region exists near the carbon fiber surface in which the polymer mass density is different than that of the bulk mass density. These surface effects extend ~10 Å into the polymer from the center of the outermost graphite layer. Early simulations predict polymer residual stress levels to be higher near the graphite surface. It is also seen that the molecular potential energy in polymer atoms decreases with increasing crosslink density. New models are then established in order to investigate the interface between EPON 862-DETDA polymer and graphene nanoplatelets (GNPs) of various atomic thicknesses. Mechanical properties are extracted from the models using Molecular Dynamics techniques. These properties are then implemented into micromechanics software that utilizes the generalized method of cells to create representations of macro-scale composites. Micromechanics models are created representing GNP doped epoxy with varying number of graphene layers and interfacial polymer crosslink densities. The initial micromechanics results for the GNP doped epoxy are then taken to represent the matrix component and are re-run through the micromechanics software with the addition of a carbon fiber to simulate a GNP doped epoxy/carbon fiber composite. Micromechanics results agree well with experimental data, and indicate GNPs of 1 to 2 atomic layers to be highly favorable. The effect of oxygen bonded to the surface of the GNPs is lastly investigated. Molecular Models are created for systems with varying graphene atomic thickness, along with different amounts of oxygen species attached to them. Models are created for graphene containing hydroxyl groups only, epoxide groups only, and a combination of epoxide and hydroxyl groups. Results show models of oxidized graphene to decrease in both tensile and shear modulus. Attaching only epoxide groups gives the best results for mechanical properties, though pristine graphene is still favored.

Ab Initio Quantum Chemical Studies on Neutral-Radical Reactions of Ethynyl (C2H) and Cyano (CN) with Unsaturated Hydrocarbons

An Ab Initio/RRKM study of the reaction mechanism and product branching ratios of neutral-radical ethynyl (C2H) and cyano (CN) radical species with unsaturated hydrocarbons is performed. The reactions studied apply to cold conditions such as planetary atmospheres including Titan, the Interstellar Medium (ISM), icy bodies and molecular clouds. The reactions of C2H and CN additions to gaseous unsaturated hydrocarbons are an active area of study. NASA’s Cassini/Huygens mission found a high concentration of C2H and CN from photolysis of ethyne (C2H2) and hydrogen cyanide (HCN), respectively, in the organic haze layers of the atmosphere of Titan. The reactions involved in the atmospheric chemistry of Titan lead to a vast array of larger, more complex intermediates and products and may also serve as a chemical model of Earth’s primordial atmospheric conditions. The C2H and CN additions are rapid and exothermic, and often occur barrierlessly to various carbon sites of unsaturated hydrocarbons. The reaction mechanism is proposed on the basis of the resulting potential energy surface (PES) that includes all the possible intermediates and transition states that can occur, and all the products that lie on the surface. The B3LYP/6-311g(d,p) level of theory is employed to determine optimized electronic structures, moments of inertia, vibrational frequencies, and zero-point energy. They are followed by single point higher-level CCSD(T)/cc-vtz calculations, including extrapolations to complete basis sets (CBS) of the reactants and products. A microcanonical RRKM study predicts single-collision (zero-pressure limit) rate constants of all reaction paths on the potential energy surface, which is then used to compute the branching ratios of the products that result. These theoretical calculations are conducted either jointly or in parallel to experimental work to elucidate the chemical composition of Titan’s atmosphere, the ISM, and cold celestial bodies.

Modelling UV response of biomolecules: from electronic structure to ab-initio dynamics.

The interaction of organic chromophores with light initiates ultrafast processes in the timescale of femtoseconds. An atomistic understanding of the mechanism driving such photoinduced reactions opens up the door to exploit them for our benefit. This thesis studies the interactions of ultraviolet light with the DNA/RNA molecules and the amino-acid tryptophan. Using some of the most accurate electronic structure methods and sophisticated environmental modelling, the works documented herein enable quantitative comparisons with cutting-edge experimental data. The relaxation pathways undertaken by the excited molecule are revealed through static and dynamical investigations of the excited-state potential energy surface. The profound role played by the dynamic response of the environment to guide the excitation in these timescales is addressed thoroughly.

Molecular modelling of organic functional materials

The investigation of the mechanisms lying behind the (photo-)chemical processes is fundamental to address and improve the design of new organic functional materials. In many cases, dynamics simulations represent the only tool to capture the system properties emerging from complex interactions between many molecules. Despite the outstanding progresses in calculation power, the only way to carry out such computational studies is to introduce several approximations with respect to a fully quantum mechanical (QM) description. This thesis presents an approach that combines QM calculations with a classical Molecular Dynamics (MD) approach by means of accurate QM-derived force fields. It is based on a careful selection of the most relevant molecular degrees of freedom, whose potential energy surface is calculated at QM level and reproduced by the analytic functions of the force field, as well as by an accurate tuning of the approximations introduced in the model of the process to be simulated. This is made possible by some tools developed purposely, that allow to obtain and test the FF parameters through comparison with the QM frequencies and normal modes. These tools were applied in the modelling of three processes: the npi* photoisomerisation of azobenzene, where the FF description was extended to the excited state too and the non-adiabatic events were treated stochastically with Tully fewest switching algorithm; the charge separation in donors-acceptors bulk heterojunction organic solar cells, where a tight-binding Hamiltonian was carefully parametrised and solved by means of a code, also written specifically; the effect of the protonation state on the photoisomerisation quantum yield of the aryl-azoimidazolium unit of the axle molecule of a rotaxane molecular shuttle. In each case, the QM-based MD models that were specifically developed gave noteworthy information about the investigated phenomena, proving to be a fundamental key for a deeper comprehension of several experimental evidences.

How does basis set superposition error change the potential surfaces for hydrogen-bonded dimers?

We describe a simple method to automate the geometric optimization of molecular orbital calculations of supermolecules on potential surfaces that are corrected for basis set superposition error using the counterpoise (CP) method. This method is applied to the H-bonding complexes HF/HCN, HF/H2O, and HCCH/H2O using the 6-31G(d,p) and D95 + + (d,p) basis sets at both the Hartree-Fock and second-order Møller-Plesset levels. We report the interaction energies, geometries, and vibrational frequencies of these complexes on the CP-optimized surfaces; and compare them with similar values calculated using traditional methods, including the (more traditional) single point CP correction. Upon optimization on the CP-corrected surface, the interaction energies become more negative (before vibrational corrections) and the H-bonding stretching vibrations decrease in all cases. The extent of the effects vary from extremely small to quite large depending on the complex and the calculational method. The relative magnitudes of the vibrational corrections cannot be predicted from the H-bond stretching frequencies alone

How does basis set superposition error change the potential surfaces for hydrogen-bonded dimers?

We describe a simple method to automate the geometric optimization of molecular orbital calculations of supermolecules on potential surfaces that are corrected for basis set superposition error using the counterpoise (CP) method. This method is applied to the H-bonding complexes HF/HCN, HF/H2O, and HCCH/H2O using the 6-31G(d,p) and D95 + + (d,p) basis sets at both the Hartree-Fock and second-order Møller-Plesset levels. We report the interaction energies, geometries, and vibrational frequencies of these complexes on the CP-optimized surfaces; and compare them with similar values calculated using traditional methods, including the (more traditional) single point CP correction. Upon optimization on the CP-corrected surface, the interaction energies become more negative (before vibrational corrections) and the H-bonding stretching vibrations decrease in all cases. The extent of the effects vary from extremely small to quite large depending on the complex and the calculational method. The relative magnitudes of the vibrational corrections cannot be predicted from the H-bond stretching frequencies alone

Analysis of confinement potential fluctuation and band-gap renormalization effects on excitonic transition in GaAs/AlGaAs multiquantum wells grown on (100) and (311)A GaAs surfaces

The competition between confinement potential fluctuations and band-gap renormalization (BGR) in GaAs/AlxGa1-xAs quantum wells grown on [1 0 0] and [3 1 1]A GaAs substrates is evaluated. The results clearly demonstrate the coexistence of the band-tail states filling related to potential fluctuations and the band-gap renormalization caused by an increase in the density of photogenerated carriers during the photoluminescence (PL) experiments. Both phenomena have strong influence on temperature dependence of the PL-peak energy (E-PL(T)). As the photon density increases, the E-PL can shift to either higher or lower energies, depending on the sample temperature. The temperature at which the displacement changes from a blueshift to a redshift is governed by the magnitude of the potential fluctuations and by the variation of BGR with excitation density. A simple band-tail model with a Gaussian-like distribution of the density of state was used to describe the competition between the band-tail filling and the BGR effects on E-PL(T). (C) 2012 Elsevier B.V. All rights reserved.

Self-supporting CVD diamond charge state conversion surfaces for high resolution imaging of low-energy neutral atoms in space plasmas

Two polycrystalline diamond surfaces, manufactured by chemical vapour deposition (CVD) technique, are investigated regarding their applicability as charge state conversion surfaces (CS) for use in a low energy neutral atom imaging instrument in space research. The capability of the surfaces for converting neutral atoms into negative ions via surface ionisation processes was measured for hydrogen and oxygen with particle energies in the range from 100 eV to 1 keV and for angles of incidence between 6 deg and 15 deg. We observed surface charging during the surface ionisation processes for one of the CVD samples due to low electrical conductivity of the material. Measurements on the other CVD diamond sample resulted in ionisation efficiencies of ~2 % for H and up to 12 % for O. Analysis of the angular scattering revealed very narrow and almost circular scattering distributions. Comparison of the results with the data of the CS of the IBEX-Lo sensor shows that CVD diamond has great potential as CS material for future space missions.

On the energy harvesting potential of piezoaeroelastic systems

This paper investigates the concept of piezoaeroelasticity for energy harvesting. The focus is placed on mathematical modeling and experimental validations of the problem of generating electricity at the flutter boundary of a piezoaeroelastic airfoil. An electrical power output of 10.7 mW is delivered to a 100 k load at the linear flutter speed of 9.30 m/s (which is 5.1% larger than the short-circuit flutter speed). The effect of piezoelectric power generation on the linear flutter speed is also discussed and a useful consequence of having nonlinearities in the system is addressed. (C) 2010 American Institute of Physics. [doi:10.1063/1.3427405]

Chemical potential and the nature of dark energy: The case of a phantom field

The influence of a possible nonzero chemical potential mu on the nature of dark energy is investigated by assuming that the dark energy is a relativistic perfect simple fluid obeying the equation of state, p=omega rho (omega < 0, constant). The entropy condition, S >= 0, implies that the possible values of omega are heavily dependent on the magnitude, as well as on the sign of the chemical potential. For mu > 0, the omega parameter must be greater than -1 (vacuum is forbidden) while for mu < 0 not only the vacuum but even a phantomlike behavior (omega <-1) is allowed. In any case, the ratio between the chemical potential and temperature remains constant, that is, mu/T=mu(0)/T(0). Assuming that the dark energy constituents have either a bosonic or fermionic nature, the general form of the spectrum is also proposed. For bosons mu is always negative and the extended Wien's law allows only a dark component with omega <-1/2, which includes vacuum and the phantomlike cases. The same happens in the fermionic branch for mu < 0. However, fermionic particles with mu > 0 are permitted only if -1 energy (omega less than or similar to-10/21).

Three-dimensional analysis of fracture, corrosion and wear surfaces

A brief look at the history of fractography has shown a recent trend in the quantification of topographic parameters through the use of three-dimensional reconstruction techniques, which associate SEM stereoscopy and stereophotogrammetry software, allowing the calculation of the elevation measurement at numerous points of the topography due to the parallax that takes place during the tilting of the sample along the microscope eucentric plane. Several investigators have used reconstruction techniques to correlate some fractographic parameters, such as fractal dimension and fractured to projected area ratio, to the mechanical properties of materials, such as fracture toughness and tensile strength. So far, the search for a clear relationship between the fracture topography and mechanical properties has provided ambiguous results. The present work applied a surface metrology software to reconstruct three-dimensionally fracture surfaces (transgranular cleavage, intergranular and dimple fracture), corrosion pits and tribo-surfaces in order to explore the potential of this stereophotogrammetry technique. The existence of a variation in the calculated topographic parameters with the conditions of SEM image acquisition reinforces the importance of both good image acquisition and accurate calibration methods in order to validate this 3D reconstruction technique in metrological terms. Preliminary results did not indicate the existence of a clear relationship between either the true to project area ratio and CVN absorbed energy or the fractal dimension and CVN absorbed energy. It is likely that each fracture mechanism presents a proper relationship between the fractographic parameters and mechanical properties. (C) 2009 Elsevier Ltd. All rights reserved.

Alterations in energy properties of eucalyptus wood and bark subjected to torrefaction: The potential of mass loss as a synthetic indicator

Torrefaction is a mild pyrolysis process (usually up to 300 degrees C) that changes the chemical and physical properties of biomass. This process is a possible pre-treatment prior to further processes (transport, grinding, combustion, gasification, etc) to generate energy or biofuels. In this study, three eucalyptus wood species and bark were subjected to different torrefaction conditions to determine the alterations in their structural and energy properties. The most severe treatment (280 degrees C, 5 h) causes mass losses of more than 35%, with severe damage to anatomical structure, and an increase of about 27% in the specific energy content. Bark is more sensitive to heat than wood. Energy yields are always higher than mass yields, thereby demonstrating the benefits of torrefaction in concentrating biomass energy. The overall mass loss is proposed as a relevant parameter to synthesize the effect of torrefaction conditions (temperature and duration). Accordingly, all results are summarised by analytical expressions able to predict the energy properties as a function of the overall mass loss. These expressions are intended to be used in any optimization procedure, from production in the field to the final use. (c) 2010 Elsevier Ltd. All rights reserved.

Variation in the dispersal potential of non-feeding invertebrate larvae: the desperate larva hypothesis and larval size

For many species of marine invertebrates, variability in larval settlement behaviour appears to be the rule rather than the exception. This variability has the potential to affect larval dispersal, because settlement behaviour will influence the length of time larvae are in the plankton. Despite the ubiquity and importance of this variability, relatively few sources of variation in larval settlement behaviour have been identified. One important factor that can affect larval settlement behaviour is the nutritional state of larvae. Non-feeding larvae often become less discriminating in their 'choice' of settlement substrate, i.e. more desperate to settle, when energetic reserves run low. We tested whether variation in larval size (and presumably in nutritional reserves) also affects the settlement behaviour of 3 species of colonial marine invertebrate larvae, the bryozoans Bugula neritina and Watersipora subtorquata and the ascidian Diplosoma listerianum. For all 3 species, larger larvae delayed settlement for longer in the absence of settlement cues, and settlement of Bugula neritina larvae was accelerated by the presence of settlement cues, independently of larval size. In the field, larger W subtorquata larvae also took longer to settle than smaller larvae and were more discriminating towards settlement surfaces. These differences in settlement time are likely to result in differences in the distance that larvae disperse in the field. We suggest that species that produce non-feeding larvae can affect the dispersal potential of their offspring by manipulating larval size and thus larval desperation.