Nuclear Hartree-Fock calculations on parallel computers

Hartree-Fock (HF) calculations have had remarkable success in describing large nuclei at high spin, temperature and deformation. To allow full range of possible deformations, the Skyrme HF equations can be discretized on a three-dimensional mesh. However, such calculations are currently limited by the computational resources provided by traditional supercomputers. To take advantage of recent developments in massively parallel computing technology, we have implemented the LLNL Skyrme-force static and rotational HF codes on Intel's DELTA and GAMMA systems at Caltech.

We decomposed the HF code by assigning a portion of the mesh to each node, with nearest neighbor meshes assigned to nodes connected by communication· channels. This kind of decomposition is well-suited for the DELTA and the GAMMA architecture because the only non-local operations are wave function orthogonalization and the boundary conditions of the Poisson equation for the Coulomb field.

Our first application of the HF code on parallel computers has been the study of identical superdeformed (SD) rotational bands in the Hg region. In the last ten years, many SD rotational bands have been found experimentally. One very surprising feature found in these SD rotational bands is that many pairs of bands in nuclei that differ by one or two mass units have nearly identical deexcitation gamma-ray energies. Our calculations of the five rotational bands in ^(192)Hg and ^(194)Pb show that the filling of specific orbitals can lead to bands with deexcitation gamma-ray energies differing by at most 2 keV in nuclei differing by two mass units and over a range of angular momenta comparable to that observed experimentally. Our calculations of SD rotational bands in the Dy region also show that twinning can be achieved by filling or emptying some specific orbitals.

The interpretation of future precise experiments on atomic parity nonconservation (PNC) in terms of parameters of the Standard Model could be hampered by uncertainties in the atomic and nuclear structure. As a further application of the massively parallel HF calculations, we calculated the proton and neutron densities of the Cesium isotopes from A = 125 to A = 139. Based on our good agreement with experimental charge radii, binding energies, and ground state spins, we conclude that the uncertainties in the ratios of weak charges are less than 10^(-3), comfortably smaller than the anticipated experimental error.

Interplay of martensitic phase transformation and plastic slip in polycrystals

Inspired by key experimental and analytical results regarding Shape Memory Alloys (SMAs), we propose a modelling framework to explore the interplay between martensitic phase transformations and plastic slip in polycrystalline materials, with an eye towards computational efficiency. The resulting framework uses a convexified potential for the internal energy density to capture the stored energy associated with transformation at the meso-scale, and introduces kinetic potentials to govern the evolution of transformation and plastic slip. The framework is novel in the way it treats plasticity on par with transformation.

We implement the framework in the setting of anti-plane shear, using a staggered implicit/explict update: we first use a Fast-Fourier Transform (FFT) solver based on an Augmented Lagrangian formulation to implicitly solve for the full-field displacements of a simulated polycrystal, then explicitly update the volume fraction of martensite and plastic slip using their respective stick-slip type kinetic laws. We observe that, even in this simple setting with an idealized material comprising four martensitic variants and four slip systems, the model recovers a rich variety of SMA type behaviors. We use this model to gain insight into the isothermal behavior of stress-stabilized martensite, looking at the effects of the relative plastic yield strength, the memory of deformation history under non-proportional loading, and several others.

We extend the framework to the generalized 3-D setting, for which the convexified potential is a lower bound on the actual internal energy, and show that the fully implicit discrete time formulation of the framework is governed by a variational principle for mechanical equilibrium. We further propose an extension of the method to finite deformations via an exponential mapping. We implement the generalized framework using an existing Optimal Transport Mesh-free (OTM) solver. We then model the $\alpha$--$\gamma$ and $\alpha$--$\varepsilon$ transformations in pure iron, with an initial attempt in the latter to account for twinning in the parent phase. We demonstrate the scalability of the framework to large scale computing by simulating Taylor impact experiments, observing nearly linear (ideal) speed-up through 256 MPI tasks. Finally, we present preliminary results of a simulated Split-Hopkinson Pressure Bar (SHPB) experiment using the $\alpha$--$\varepsilon$ model.

A model for energy and morphology of crystalline grain boundaries with arbitrary geometric character

It has been well-established that interfaces in crystalline materials are key players in the mechanics of a variety of mesoscopic processes such as solidification, recrystallization, grain boundary migration, and severe plastic deformation. In particular, interfaces with complex morphologies have been observed to play a crucial role in many micromechanical phenomena such as grain boundary migration, stability, and twinning. Interfaces are a unique type of material defect in that they demonstrate a breadth of behavior and characteristics eluding simplified descriptions. Indeed, modeling the complex and diverse behavior of interfaces is still an active area of research, and to the author's knowledge there are as yet no predictive models for the energy and morphology of interfaces with arbitrary character. The aim of this thesis is to develop a novel model for interface energy and morphology that i) provides accurate results (especially regarding "energy cusp" locations) for interfaces with arbitrary character, ii) depends on a small set of material parameters, and iii) is fast enough to incorporate into large scale simulations.

In the first half of the work, a model for planar, immiscible grain boundary is formulated. By building on the assumption that anisotropic grain boundary energetics are dominated by geometry and crystallography, a construction on lattice density functions (referred to as "covariance") is introduced that provides a geometric measure of the order of an interface. Covariance forms the basis for a fully general model of the energy of a planar interface, and it is demonstrated by comparison with a wide selection of molecular dynamics energy data for FCC and BCC tilt and twist boundaries that the model accurately reproduces the energy landscape using only three material parameters. It is observed that the planar constraint on the model is, in some cases, over-restrictive; this motivates an extension of the model.

In the second half of the work, the theory of faceting in interfaces is developed and applied to the planar interface model for grain boundaries. Building on previous work in mathematics and materials science, an algorithm is formulated that returns the minimal possible energy attainable by relaxation and the corresponding relaxed morphology for a given planar energy model. It is shown that the relaxation significantly improves the energy results of the planar covariance model for FCC and BCC tilt and twist boundaries. The ability of the model to accurately predict faceting patterns is demonstrated by comparison to molecular dynamics energy data and experimental morphological observation for asymmetric tilt grain boundaries. It is also demonstrated that by varying the temperature in the planar covariance model, it is possible to reproduce a priori the experimentally observed effects of temperature on facet formation.

Finally, the range and scope of the covariance and relaxation models, having been demonstrated by means of extensive MD and experimental comparison, future applications and implementations of the model are explored.

Petrografia, litogeoquímica e datação Ar-Ar dos montes submarinos e dos rochedos de Martim Vaz - Cadeia Vitória-Trindade

Este trabalho de mestrado estudou a ilha de Martin Vaz e cinco montes submarinos da Cadeia Vitória-Trindade Columbia, Dogaressa, Davis, Jaseur e Montague. Martin Vaz é um conjunto de ilhas formado pela ilha principal - Martin Vaz, duas ilhotas íngremes e inacessíveis - a Ilha do Norte e a Ilha do Sul além de vários rochedos menores, como o Rochedo Agulha, espalhados a 48 km a leste de Trindade, perfazendo uma área total de 0,3 km. Martin Vaz, assim como os montes submarinos, pertencem a chamada Cadeia Vitória-Trindade, estão inseridos em um contexto tectônico regional cujo trend W-E sugere representar o track da pluma mantélica de Trindade quando da passagem da Placa Sul Americana sobre ela desde o Terciário (CROUGH et al., 1980; OCONNOR & DUNCAN, 1990, GIBSON et at., 1997). A petrografia das amostras de Martin Vaz indica haver basanitos parcialmente alteradas, melanocráticas, textura afanítica, porosas, apresentando vesículas em torno de 1,0-5,0 milímetros. Apresenta fenocristais de piroxênio além de alguns fenocristais de olivina verde-oliva translúcido variando de 1,0-3,0 milímetros. A ilha principal apresenta também diques e necks fonolíticos apresentando matriz microlítica alterada, orientada, de cor verde apresentando minerais ripiformes de cor branca (feldspato alcalino) e outros de cor violácea (titanoaugita) além de pequenos opacos. Pequenos fenocristais de aegerina-augita fortemente pleocroica, alguns apresentando geminação simples, por vezes zonado, apresenta extinção variando de c ∧ α ou X = 23 a 33 (medida de 10 grãos). Biotita laranja amarronzada com textura poiquilítica (1,0 mm), minúsculos cristais euédricos de titanita (raros), além de cristais pseudohexagonais isotrópicos alterados de analcita e carbonatos. As amostras utilizadas neste trabalho de mestrado possuem valor mínimo de 33.91 % SiO2 (TRIM-01D) e máximo de 52,2 (MVA-01) variando de ultrabásicas a básicas. Através da análise dos óxidos SiO2 e MgO é possível distinguir dois grupos de rochas para Martin Vaz: um ultramáfico magnesiano (<42% SiO2 e >7% MgO) e um básico (>45% SiO2) e, para os montes submarinos, dois grupos: um ultramáfico magnesiano (>9% MgO <42% SiO2) e um básico (>45% SiO2 e com valores de MgO em torno de 4%). As análises de Ar-Ar para as quatro amostras de Martin Vaz apresentam idades para o derrame de basanito variando de 320366 Ka (MVA-10) à 623127 Ka (MVA-04). A única amostra datada representando do dique de fonólito é a MVA-05B e obteve idade de 64984 Ka, indicando ser contemporânea ao derrame basanítico.

Bundles of polytwins as meta-elastic domains in the thin polycrystalline simple multi-ferroic system PZT.

We used enhanced piezo-response force microscopy (E-PFM) to investigate both ferroelastic and ferroelectric nanodomains in thin films of the simple multi-ferroic system PbZr(0.3)Ti(0.7)O(3) (PZT). We show how the grains are organized into a new type of elastic domain bundles of the well-known periodic elastic twins. Here we present these bundle domains and discuss their stability and origin. Moreover, we show that they can arrange in such a way as to release strain in a more effective way than simple twinning. Finally, we show that these bundle domains can arrange to form the macroscopic ferroelectric domains that constitute the basis of ferroelectric-based memory devices.

Twin-free uniform epitaxial GaAs nanowires grown by a two-temperature process.

We demonstrate vertically aligned epitaxial GaAs nanowires of excellent crystallographic quality and optimal shape, grown by Au nanoparticle-catalyzed metalorganic chemical vapor deposition. This is achieved by a two-temperature growth procedure, consisting of a brief initial high-temperature growth step followed by prolonged growth at a lower temperature. The initial high-temperature step is essential for obtaining straight, vertically aligned epitaxial nanowires on the (111)B GaAs substrate. The lower temperature employed for subsequent growth imparts superior nanowire morphology and crystallographic quality by minimizing radial growth and eliminating twinning defects. Photoluminescence measurements confirm the excellent optical quality of these two-temperature grown nanowires. Two mechanisms are proposed to explain the success of this two-temperature growth process, one involving Au nanoparticle-GaAs interface conditions and the other involving melting-solidification temperature hysteresis of the Au-Ga nanoparticle alloy.

Nearly intrinsic exciton lifetimes in single twin-free GaAsAlGaAs core-shell nanowire heterostructures

CW and time-resolved photoluminescence measurements are used to investigate exciton recombination dynamics in GaAsAlGaAs heterostructure nanowires grown with a recently developed technique which minimizes twinning. A thin capping layer is deposited to eliminate the possibility of oxidation of the AlGaAs shell as a source of oxygen defects in the GaAs core. We observe exciton lifetimes of ∼1 ns, comparable to high quality two-dimensional double heterostructures. These GaAs nanowires allow one to observe state filling and many-body effects resulting from the increased carrier densities accessible with pulsed laser excitation. © 2008 American Institute of Physics.

Modelling and design of stress-induced martensite formation in metastable β Ti alloys

The temperature dependence of the stress-induced martensite (SIM) formation in a Ti-10V-2Fe-3Al (Ti-1023) alloy under compressive loading has been studied. At low temperatures, the stress level at which martensite starts to form increases linearly with the deformation temperature, while the stress at which the deformation switches to regular plastic deformation is roughly temperature independent. A thermostatistical model for dislocation evolution is employed to describe deformation twinning in martensite. Combined effects of twinning induced plasticity and solid solution strengthening are considered in terms of temperature variations. The SIM effect disappears on deformation at temperatures beyond ~ 233 ° C, which is close to the predicted Ms temperature of 240°C. The thermostatistical model predicts a transition from twinned martensite to pure slip at 250°C. By providing a model to predict the martensite formation, and by describing deformation twinning, the present work provides a number of tools that may be employed to conceive new titanium alloys combining improved strength and ductility. © 2013 Elsevier B.V.

Tensile and compressive mechanical behavior of twinned silicon carbide nanowires

Molecular dynamics simulations with the Tersoff potential were used to study the response of twinned SiC nanowires under tensile and compressive strain. The critical strain of the twinned nanowires can be enhanced by twin stacking faults, and their critical strains are larger than those of perfect nanowires with the same diameters. Under axial tensile strain, the bonds of the nanowires are stretched just before failure. The failure behavior is found to depend on the twin segment thickness and the diameter of the nanowires. An atomic chain is observed for thin nanowires with small twin segment thickness under tension strain. Under axial compressive strain, the collapse of twinned SiC nanowires exhibits two different failure modes, depending on the length and diameter of the nanowires, i.e., shell buckling for short nanowires and columnar buckling for longer nanowires.

Selective Synthesis of Single-Crystalline Rhombic Dodecahedral, Octahedral, and Cubic Gold Nanocrystals

This paper reports a versatile seed-mediated growth method for selectively synthesizing single-crystalline rhombic dodecahedral, octahedral, and cubic gold nanocrystals. In the seed-mediated growth method, cetylpyridinium chloride (CPC) and CPC-capped single-crystalline gold nanocrystals 41.3 nm in size are used as the surfactant and seeds, respectively. The CPC-capped gold seeds can avoid twinning during the growth process, which enables us to study the correlations between the growth conditions and the shapes of the gold nanocrystals. Surface-energy and kinetic considerations are taken into account to understand the formation mechanisms of the single-crystalline gold nanocrystals with varying shapes.

Template-free fabrication of hexagonal ZnO microprism with an interior space

Well-faceted hexagonal ZnO microprisms with regular interior space have been successfully prepared by a template-free hydrothermal synthetic route. The morphologies of the products depend on the experimental conditions such as the solvent, the concentration of ammonia aqueous solution, and the reaction temperature. Through manipulation of the aging time, the as-prepared ZnO can be controlled as a monodispersed hexagonal twinning solid or as hollow microprisms. Moreover, the evolution process of the hollow ZnO nanoarchitecture after reaction for 2, 6, 12, and 24 h has been investigated by field emission scanning electron microscopy (FE-SEM) and transmission electron microscopy (TEM). A possible growth mechanism has also been proposed and discussed. Furthermore, the photoluminescence (PL) measurement exhibits the unique emitting characteristic of hollow ZnO nanostructures.

大港油田复杂断块油气藏评价技术研究与应用

In this paper, in view of characteristics and appraisal difficulty of complicate fault block reservoir, based on the theory of reservoir geology, tectonic geology, sequence stratigraphy, logging geology and sedimentology, according to related multidisciplinary data, such as geological, geophysical, logging and testing, taking 3D seismic acquisition processing, delicate structure interpretation, fine reservoir prediction, reservoir engineering and logging evaluation as research methods, on the basis of comprehensive geologic research, innovative appraisal thinking of complicate fault block reservoir is probed into, reservoir evaluation technology series of different reservoirs were created, and appraisal circuit of complicate fault block reservoir was proposed. According to the research on main controlling factors of hydrocarbon and concentration and concealing laws of Dagang exploration area, family assemble patterns of complicate fault block reservoir were proposed firstly, such as twinning, juxtaposed, overlying, concurrent, connected in series and so on, so theory basement was provided to reservoir evaluation. Taking into account of difficulties of low signal to noise ratio and border identification of lithologic trap of complicate fault block reservoir, in the middle of seismic data acquisition, double focusing layout and optimization technology were adopted, geophysical model was set up, and best observation system was ascertained. During object processing, matching processing technology of removing time difference, phase difference, energy difference and frequency spectrum difference generated by different blazing and receiving elements was studied. Imagery, low amplitude structure and identification of subtle lithologic reservoir was interpreted reasonably. On the basis of characteristics and technical description analysis of structural, structural--lithologic and lithologic reservoirs, innovative appraisal thinking of these reservoirs was formatted, appraisal circuit of complicate fault block reservoir was proposed, and this method could be used in other similar oilfields

Origin of Ferroelastic domains in Free-Standing Single Crystal Ferroelectric Films

The origin of the unusual 90 degrees ferroelectric/ferroelastic domains, consistently observed in recent studies on mesoscale and nanoscale free-standing single crystals of BaTiO3 [Schilling , Phys. Rev. B 74, 024115 (2006); Schilling , Nano Lett. 7, 3787 (2007)], has been considered. A model has been developed which postulates that the domains form as a response to elastic stress induced by a surface layer which does not undergo the paraelectric-ferroelectric cubic-tetragonal phase transition. This model was found to accurately account for the changes in domain periodicity as a function of size that had been observed experimentally. The physical origin of the surface layer might readily be associated with patterning damage, seen in experiment; however, when all evidence of physical damage is removed from the BaTiO3 surfaces by thermal annealing, the domain configuration remains practically unchanged. This suggests a more intrinsic origin, such as the increased importance of surface tension at small dimensions. The effect of surface tension is also shown to be proportional to the difference in hardness between the surface and the interior of the ferroelectric. The present model for surface-tension induced twinning should also be relevant for finely grained or core-shell structured ceramics.

Atomistic studies of interactions between the dominant lattice dislocations and γ/γ-lamellar boundaries in lamellar γ-TiAl

This paper reports on atomistic simulations of the interactions between the dominant lattice dislocations in ?-TiAl (<1 0 1] superdislocations) with all three kinds of ?/?-lamellar boundaries in polysynthetically twinned (PST) TiAl. The purpose of this study is to clarify the early stage of lamellar boundary controlled plastic deformation in PST TiAl. The interatomic interactions in our simulations are described by a bond order potential for L10-TiAl which provides a proper quantum mechanical description of the bonding. We are interested in the dislocation core geometries that the lattice produces in proximity to lamellar boundaries and the way in which these cores are affected by the elastic and atomistic effects of dislocation-lamellar boundary interaction. We study the way in which the interfaces affect the activation of ordinary dislocation and superdislocation slip inside the ?-lamellae and transfer of plastic deformation across lamellar boundaries. We find three new phenomena in the atomic-scale plasticity of PST TiAl, particularly due to elastic and atomic mismatch associated with the 60° and 120° ?/?-interfaces: (i) two new roles of the ?/?-interfaces, i.e. decomposition of superdislocations within 120° and 60° interfaces and subsequent detachment of a single ordinary dislocation and (ii) blocking of ordinary dislocations by 60° and 120° interfaces resulting in the emission of a twinning dislocation.

Microscopic Origin of Channeled Flow in Lamellar Titanium Aluminide

We employ a quantum mechanical bond order potential in an atomistic simulation of channeled flow. We show that the original hypothesis that this is achieved by a cooperative deployment of slip and twinning is correct, first because a twin is able to “protect” a 60° ordinary dislocation from becoming sessile, and second because the two processes are found to be activated by Peierls stresses of similar magnitude. In addition we show an explicit demonstration of the lateral growth of a twin, again at a similar level of stress. Thus these simultaneous processes are shown to be capable of channeling deformation into the observed state of plane strain in so-called “A”-oriented mechanical testing of titanium aluminide superalloy.