Commentaria in primam Diui thomae partem : in duos tomos diuisa ... /

Mode of access: Internet.

Miocene stable isotope record of IODP Site 320-U1335 and revised stable isotope record of IODP Site 321-U1337

The Miocene Climatic Optimum (MCO; ~16.9 to 14.7 Ma) provides an outstanding opportunity to investigate climate-carbon cycle dynamics during a geologically recent interval of global warmth. We present benthic stable oxygen (d18O) and carbon (d13C) isotope records (5-12 kyr time resolution) spanning the late early to middle Miocene interval (18 to 13 Ma) at Integrated Ocean Drilling Program (IODP) Site U1335 (eastern equatorial Pacific Ocean). The U1335 stable isotope series track the onset and development of the MCO as well as the transitional climatic phase culminating with global cooling and expansion of the East Antarctic ice-sheet at ~13.8 Ma. We integrate these new data with published stable isotope, geomagnetic polarity and X-ray fluorescence (XRF) scanner-derived carbonate records from IODP Sites U1335, U1336, U1337 and U1338 on a consistent, astronomically-tuned timescale. Benthic isotope and XRF scanner-derived CaCO3 records depict prominent 100 kyr variability with 400 kyr cyclicity additionally imprinted on d13C and CaCO3 records, pointing to a tight coupling between the marine carbon cycle and climate variations. Our inter-site comparison further indicates that the lysocline behaved in highly dynamic manner throughout the MCO, with >75% carbonate loss occurring at paleo-depths ranging from ~3.4 to ~4 km in the eastern equatorial Pacific Ocean. Carbonate dissolution maxima coincide with warm phases (d18O minima) and d13C decreases, implying that climate-carbon cycle feedbacks fundamentally differed from the late Pleistocene glacial-interglacial pattern, where dissolution maxima correspond to d13C maxima and d18O minima. Carbonate dissolution cycles during the MCO were, thus, more similar to Paleogene hyperthermal patterns.

Changes in organic matter production and accumulation as a mechanism for isotopic evolution in the Mesoproterozoic ocean

Mesoproterozoic marine successions worldwide record a shift in average delta(13)C values from 0 to +3.5parts per thousand, with the latter value evident in successions younger than 1250 Ma. New carbon isotope data from the similar to 1300 to 1270 Ma Dismal Lakes Group, Arctic Canada, provide further insight into this fundamental transition. Data reveal that the shift to higher VC values was gradual and marked by occasional excursions to values less than 0 parts per thousand. When compared to records from older and younger marine successions, it is evident that the difference between isotopic minima and maxima increased with time, indicating that the marine system evolved to become isotopically more variable. We interpret these patterns to record an increase in the crustal inventory of organic carbon, reflecting eukaryotic diversification and a change in the locus of organic carbon burial to include anoxic deep marine sites where preservation potential was high. We speculate that the release of O-2 to Earth's surface environments associated with increased organic carbon storage induced irreversible changes in the Mesoproterozoic biosphere, presaging the more extreme environmental and evolutionary developments of the Neoproterozoic.

Genetic algorithms in statistical tolerancing

An important aspect in manufacturing design is the distribution of geometrical tolerances so that an assembly functions with given probability, while minimising the manufacturing cost. This requires a complex search over a multidimensional domain, much of which leads to infeasible solutions and which can have many local minima. As well, Monte-Carlo methods are often required to determine the probability that the assembly functions as designed. This paper describes a genetic algorithm for carrying out this search and successfully applies it to two specific mechanical designs, enabling comparisons of a new statistical tolerancing design method with existing methods. (C) 2003 Elsevier Ltd. All rights reserved.

Partial regularity of weak solutions of the liquid crystal equilibrium system

For a parameter, we consider the modified relaxed energy of the liquid crystal system. Each minimizer of the modified relaxed energy is a weak solution to the liquid crystal equilibrium system. We prove the partial regularity of minimizers of the modified relaxed energy. We also prove the existence of infinitely many weak solutions for the special boundary value x.

High-frequency winter cooling and reef coral mortality during the Holocene climatic optimum

A detailed ecological, micro-structural and skeletal Sr/Ca study of a 3.42 m thick Goniopora reef profile from an emerged Holocene reef terrace at the northern South China Sea reveals at least nine abrupt massive Goniopora stress and mortality events occurred in winter during the 7.0-7.5 thousand calendar years before present (cal. ka BP) (within the Holocene climatic optimum). Whilst calculated Sr/Ca-SST (sea surface temperature) maxima during this period are comparable to those in the 1990s, Sr/Ca-SST minima are significantly lower, probably due to stronger winter monsoons. Such generally cooler winters, superimposed by further exceptional winter cooling on inter-annual to decadal scales, may have caused stress and mortality of the corals about every 50 years. Sea level rose by similar to 3.42 m during this period, with present sea-level reached at similar to 7.3 ka BP and a sea-level highstand of at least similar to 1.8 m occurred at similar to 7.0 ka. The results show that it took about 20-25 years for a killed Goniopora coral reef to recover. (C) 2004 Elsevier B.V. All rights reserved.

Behaviour of the energy gap in a model of Josephson coupled Bose-Einstein condensates

In this work we investigate the energy gap between the ground state and the first excited state in a model of two single-mode Bose-Einstein condensates coupled via Josephson tunnelling. The ene:rgy gap is never zero when the tunnelling interaction is non-zero. The gap exhibits no local minimum below a threshold coupling which separates a delocalized phase from a self-trapping phase that occurs in the absence of the external potential. Above this threshold point one minimum occurs close to the Josephson regime, and a set of minima and maxima appear in the Fock regime. Expressions for the position of these minima and maxima are obtained. The connection between these minima and maxima and the dynamics for the expectation value of the relative number of particles is analysed in detail. We find that the dynamics of the system changes as the coupling crosses these points.

Anomalous Excitation Spectra of Frustrated Quantum Antiferromagnets

We use series expansions to study the excitation spectra of spin-1/2 antiferromagnets on anisotropic triangular lattices. For the isotropic triangular lattice model (TLM), the high-energy spectra show several anomalous features that differ strongly from linear spin-wave theory (LSWT). Even in the Neel phase, the deviations from LSWT increase sharply with frustration, leading to rotonlike minima at special wave vectors. We argue that these results can be interpreted naturally in a spinon language and provide an explanation for the previously observed anomalous finite-temperature properties of the TLM. In the coupled-chains limit, quantum renormalizations strongly enhance the one-dimensionality of the spectra, in agreement with experiments on Cs2CuCl4.

Orbital currents and charge density waves in a generalized Hubbard ladder

We study a generalized Hubbard model on the two-leg ladder at zero temperature, focusing on a parameter region with staggered flux (SF)/d-density wave (DDW) order. To guide our numerical calculations, we first investigate the location of a SF/DDW phase in the phase diagram of the half-filled weakly interacting ladder using a perturbative renormalization group (RG) and bosonization approach. For hole doping 6 away from half-filling, finite-system density-matrix renormalizationgroup (DMRG) calculations are used to study ladders with up to 200 rungs for intermediate-strength interactions. In the doped SF/DDW phase, the staggered rung current and the rung electron density both show periodic spatial oscillations, with characteristic wavelengths 2/delta and 1/delta, respectively, corresponding to ordering wavevectors 2k(F) and 4k(F) for the currents and densities, where 2k(F) = pi(1 - delta). The density minima are located at the anti-phase domain walls of the staggered current. For sufficiently large dopings, SF/DDW order is suppressed. The rung density modulation also exists in neighboring phases where currents decay exponentially. We show that most of the DMRG results can be qualitatively understood from weak-coupling RG/bosonization arguments. However, while these arguments seem to suggest a crossover from non-decaying correlations to power-law decay at a length scale of order 1/delta, the DMRG results are consistent with a true long-range order scenario for the currents and densities. (c) 2005 Elsevier Inc. All rights reserved.

Excitation spectra of the spin-1/2 triangular-lattice Heisenberg antiferromagnet

We use series expansion methods to calculate the dispersion relation of the one-magnon excitations for the spin-(1)/(2) triangular-lattice nearest-neighbor Heisenberg antiferromagnet above a three-sublattice ordered ground state. Several striking features are observed compared to the classical (large-S) spin-wave spectra. Whereas, at low energies the dispersion is only weakly renormalized by quantum fluctuations, significant anomalies are observed at high energies. In particular, we find rotonlike minima at special wave vectors and strong downward renormalization in large parts of the Brillouin zone, leading to very flat or dispersionless modes. We present detailed comparison of our calculated excitation energies in the Brillouin zone with the spin-wave dispersion to order 1/S calculated recently by Starykh, Chubukov, and Abanov [Phys. Rev. B74, 180403(R) (2006)]. We find many common features but also some quantitative and qualitative differences. We show that at temperatures as low as 0.1J the thermally excited rotons make a significant contribution to the entropy. Consequently, unlike for the square lattice model, a nonlinear sigma model description of the finite-temperature properties is only applicable at temperatures < 0.1J. Finally, we review recent NMR measurements on the organic compound kappa-(BEDT-TTF)(2)Cu-2(CN)(3). We argue that these are inconsistent with long-range order and a description of the low-energy excitations in terms of interacting magnons, and that therefore a Heisenberg model with only nearest-neighbor exchange does not offer an adequate description of this material.

Description and training of neural network dynamics

Attractor properties of a popular discrete-time neural network model are illustrated through numerical simulations. The most complex dynamics is found to occur within particular ranges of parameters controlling the symmetry and magnitude of the weight matrix. A small network model is observed to produce fixed points, limit cycles, mode-locking, the Ruelle-Takens route to chaos, and the period-doubling route to chaos. Training algorithms for tuning this dynamical behaviour are discussed. Training can be an easy or difficult task, depending whether the problem requires the use of temporal information distributed over long time intervals. Such problems require training algorithms which can handle hidden nodes. The most prominent of these algorithms, back propagation through time, solves the temporal credit assignment problem in a way which can work only if the relevant information is distributed locally in time. The Moving Targets algorithm works for the more general case, but is computationally intensive, and prone to local minima.

The dynamics of a genetic algorithm under stabilizing selection

A formalism recently introduced by Prugel-Bennett and Shapiro uses the methods of statistical mechanics to model the dynamics of genetic algorithms. To be of more general interest than the test cases they consider. In this paper, the technique is applied to the subset sum problem, which is a combinatorial optimization problem with a strongly non-linear energy (fitness) function and many local minima under single spin flip dynamics. It is a problem which exhibits an interesting dynamics, reminiscent of stabilizing selection in population biology. The dynamics are solved under certain simplifying assumptions and are reduced to a set of difference equations for a small number of relevant quantities. The quantities used are the population's cumulants, which describe its shape, and the mean correlation within the population, which measures the microscopic similarity of population members. Including the mean correlation allows a better description of the population than the cumulants alone would provide and represents a new and important extension of the technique. The formalism includes finite population effects and describes problems of realistic size. The theory is shown to agree closely to simulations of a real genetic algorithm and the mean best energy is accurately predicted.

Using generative models for handwritten digit recognition

We describe a method of recognizing handwritten digits by fitting generative models that are built from deformable B-splines with Gaussian ``ink generators'' spaced along the length of the spline. The splines are adjusted using a novel elastic matching procedure based on the Expectation Maximization (EM) algorithm that maximizes the likelihood of the model generating the data. This approach has many advantages. (1) After identifying the model most likely to have generated the data, the system not only produces a classification of the digit but also a rich description of the instantiation parameters which can yield information such as the writing style. (2) During the process of explaining the image, generative models can perform recognition driven segmentation. (3) The method involves a relatively small number of parameters and hence training is relatively easy and fast. (4) Unlike many other recognition schemes it does not rely on some form of pre-normalization of input images, but can handle arbitrary scalings, translations and a limited degree of image rotation. We have demonstrated our method of fitting models to images does not get trapped in poor local minima. The main disadvantage of the method is it requires much more computation than more standard OCR techniques.

Topographic mappings and feed-forward neural networks

This thesis is a study of the generation of topographic mappings - dimension reducing transformations of data that preserve some element of geometric structure - with feed-forward neural networks. As an alternative to established methods, a transformational variant of Sammon's method is proposed, where the projection is effected by a radial basis function neural network. This approach is related to the statistical field of multidimensional scaling, and from that the concept of a 'subjective metric' is defined, which permits the exploitation of additional prior knowledge concerning the data in the mapping process. This then enables the generation of more appropriate feature spaces for the purposes of enhanced visualisation or subsequent classification. A comparison with established methods for feature extraction is given for data taken from the 1992 Research Assessment Exercise for higher educational institutions in the United Kingdom. This is a difficult high-dimensional dataset, and illustrates well the benefit of the new topographic technique. A generalisation of the proposed model is considered for implementation of the classical multidimensional scaling (¸mds}) routine. This is related to Oja's principal subspace neural network, whose learning rule is shown to descend the error surface of the proposed ¸mds model. Some of the technical issues concerning the design and training of topographic neural networks are investigated. It is shown that neural network models can be less sensitive to entrapment in the sub-optimal global minima that badly affect the standard Sammon algorithm, and tend to exhibit good generalisation as a result of implicit weight decay in the training process. It is further argued that for ideal structure retention, the network transformation should be perfectly smooth for all inter-data directions in input space. Finally, there is a critique of optimisation techniques for topographic mappings, and a new training algorithm is proposed. A convergence proof is given, and the method is shown to produce lower-error mappings more rapidly than previous algorithms.

Strain distributionand electronic property modifications in Si/Ge axial nanowires hetrostructures

Molecular dynamics simulations were carried out for Si/Ge axial nanowire heterostructures using modified effective atom method (MEAM) potentials. A Si–Ge MEAM interatomic cross potential was developed based on available experimental data and was used for these studies. The atomic distortions and strain distributions near the Si/Ge interfaces are predicted for nanowires with their axes oriented along the [111] direction. The cases of 10 and 25 nm diameter Si/Ge biwires and of 25 nm diameter Si/Ge/Si axial heterostructures with the Ge disk 1 nm thick were studied. Substantial distortions in the height of the atoms adjacent to the interface were found for the biwires but not for the Ge disks. Strains as high as 3.5% were found for the Ge disk and values of 2%–2.5% were found at the Si and Ge interfacial layers in the biwires. Deformation potential theory was used to estimate the influence of the strains on the band gap, and reductions in band gap to as small as 40% of bulk values are predicted for the Ge disks. The localized regions of increased strain and resulting energy minima were also found within the Si/Ge biwire interfaces with the larger effects on the Ge side of the interface. The regions of strain maxima near and within the interfaces are anticipated to be useful for tailoring band gaps and producing quantum confinement of carriers. These results suggest that nanowire heterostructures provide greater design flexibility in band structure modification than is possible with planar layer growth.