Update on the correlation of the highest energy cosmic rays with nearby extragalactic matter

Data collected by the Pierre Auger Observatory through 31 August 2007 showed evidence for anisotropy in the arrival directions of cosmic rays above the Greisen-Zatsepin-Kuz`min energy threshold, 6 x 10(19) eV. The anisotropy was measured by the fraction of arrival directions that are less than 3.1 degrees from the position of an active galactic nucleus within 75 Mpc (using the Veron-Cetty and Veron 12th catalog). An updated measurement of this fraction is reported here using the arrival directions of cosmic rays recorded above the same energy threshold through 31 December 2009. The number of arrival directions has increased from 27 to 69, allowing a more precise measurement. The correlating fraction is (38(-6)(+7))%, compared with 21% expected for isotropic cosmic rays. This is down from the early estimate of (69-(+11)(13))%. The enlarged set of arrival directions is examined also in relation to other populations of nearby extragalactic objects: galaxies in the 2 Microns All Sky Survey and active galactic nuclei detected in hard X-rays by the Swift Burst Alert Telescope. A celestial region around the position of the radiogalaxy Cen A has the largest excess of arrival directions relative to isotropic expectations. The 2-point autocorrelation function is shown for the enlarged set of arrival directions and compared to the isotropic expectation. (C) 2010 Elsevier B.V. All rights reserved.

The fluorescence detector of the Pierre Auger Observatory

The Pierre Auger Observatory is a hybrid detector for ultra-high energy cosmic rays. It combines a surface array to measure secondary particles at ground level together with a fluorescence detector to measure the development of air showers in the atmosphere above the array. The fluorescence detector comprises 24 large telescopes specialized for measuring the nitrogen fluorescence caused by charged particles of cosmic ray air showers. In this paper we describe the components of the fluorescence detector including its optical system, the design of the camera, the electronics, and the systems for relative and absolute calibration. We also discuss the operation and the monitoring of the detector. Finally, we evaluate the detector performance and precision of shower reconstructions. (C) 2010 Elsevier B.V All rights reserved.

Measurement of the energy spectrum of cosmic rays above 10(18) eV using the Pierre Auger Observatory

We report a measurement of the flux of cosmic rays with unprecedented precision and Statistics using the Pierre Auger Observatory Based on fluorescence observations in coincidence with at least one Surface detector we derive a spectrum for energies above 10(18) eV We also update the previously published energy spectrum obtained with the surface detector array The two spectra are combined addressing the systematic uncertainties and, in particular. the influence of the energy resolution on the spectral shape The spectrum can be described by a broken power law E(-gamma) with index gamma = 3 3 below the ankle which is measured at log(10)(E(ankle)/eV) = 18 6 Above the ankle the spectrum is described by a power law with index 2 6 followed by a flux suppression, above about log(10)(E/eV) = 19 5, detected with high statistical significance (C) 2010 Elsevier B V All rights reserved

A study of the effect of molecular and aerosol conditions in the atmosphere on air fluorescence measurements at the Pierre Auger Observatory

The air fluorescence detector of the Pierre Auger Observatory is designed to perforin calorimetric measurements of extensive air showers created by Cosmic rays of above 10(18) eV. To correct these measurements for the effects introduced by atmospheric fluctuations, the Observatory contains a group Of monitoring instruments to record atmospheric conditions across the detector site, ail area exceeding 3000 km(2). The atmospheric data are used extensively in the reconstruction of air showers, and are particularly important for the correct determination of shower energies and the depths of shower maxima. This paper contains a summary of the molecular and aerosol conditions measured at the Pierre Auger Observatory since the start of regular operations in 2004, and includes a discussion of the impact of these measurements oil air shower reconstructions. Between 10(18) and 10(20) eV, the systematic Uncertainties due to all atmospheric effects increase from 4% to 8% in measurements of shower energy, and 4 g cm(-2) to 8 g cm(-2) in measurements of the shower maximum. (C) 2010 Elsevier B.V. All rights reserved.

Effects of electric fields on ultracold Rydberg atom interactions

The behaviour of interacting ultracold Rydberg atoms in both constant electric fields and laser fields is important for designing experiments and constructing realistic models of them. In this paper, we briefly review our prior work and present new results on how electric fields affect interacting ultracold Rydberg atoms. Specifically, we address the topics of constant background electric fields on Rydberg atom pair excitation and laser-induced Stark shifts on pair excitation.

Precambrian history of the Zona Transversal of the Borborema Province, NE Brazil: Insights from Sm-Nd and U-Pb geochronology

The Borborema Province has three major subprovinces. The northern subprovince lies north of the Patos shear zone and is comprised of Paleoproterozoic cratonic basement with Archean nuclei, plus overlying Neoproterozoic supracrustal rocks and Brasiliano plutonic rocks. The central subprovince occurs between the Patos and Pernambuco shear zones and is mainly comprised of the Zona Transversal. The southern subprovince occurs between the Pernamabuco shear zone and the Sao Francisco craton and is comprised of a tectonic collage of various blocks, terranes, or domains ranging in age from Archean to Neoproterozoic. This report focuses on the Zona Transversal, especially on Brasiliano rocks for which we have the most new information. Paleoproterozoic gneisses with ages of 2.0-2.2 Ga occur discontinuously throughout the Zona Transversal. The Cariris Velhos suite consists of metavolcanic, metasedimentary, and metaplutonic rocks yielding U-Pb zircon ages of 995-960 Ma. This suite is mainly confined to a 100 km wide belt that extends for more than 700 km within the Alto Pajeu terrane. Sm-Nd model ages in metaigneous rocks cluster about 1.3-1.6 Ga, indicating that older crust was involved in genesis of their magmas. Brasiliano supracrustal rocks dominate the Pianco-Alto Brigida terrane, and they probably also constitute significant parts of the Alto Pajeu and Rio Capibaribe terranes. They are only slightly older than early stages of Brasiliano plutonism, with detrital zircon ages at least as young as 620 Ma; most T(DM) ages range from 1.2 to 1.6 Ga. Brasiliano plutons range from ca. 640 to 540 Ma, and their T(DM) ages range from 1.2 to 2.5 Ga. Previous workers have shown significant correlations among U-Pb ages, Sm-Nd model ages, petrology, and geochemistry, and we are able to reinforce and extend these correlations. Stage I plutons formed 640 -610 Ma and have T(DM) ages less than 1.5 Ga. Stage 11 (610-590 Ma) contains few plutons, but coincides with the peak of compressional deformation, metamorphism, and formation of migmatites. Stage III plutons (590 to ca. 575 Ma) have older T(DM) ages (ca. 1.8-2.0 Ga), as do Stage IV plutons (575 to ca. 550 Ma; T(DM) from 1.9 to 2.4 Ga). Stage III plutons formed during the transition from compressional to transcurrent deformation, while Stage IV plutons are mainly post-tectonic. Stage V plutons (550-530 Ma) are commonly undeformed (except along younger shear zones) and have A-type geochemistry. The five stages have distinct geochemical properties, which suggest that the tectonic settings evolved from early, arc-related magma-genesis (Stage I) to within-plate magma-genesis (Stage V), with perhaps some intermediate phases of extensional environments. (C) 2011 Elsevier Ltd. All rights reserved.

The Cariris Velhos tectonic event in Northeast Brazil

The Borborema Province in northeastern South America is a typical Brasiliano-Pan-African branching system of Neoproterozoic orogens that forms part of the Western Gondwana assembly. The province is positioned between the Sao Luis-West Africa craton to the north and the Sao Francisco (Congo-Kasai) craton to the south. For this province the main characteristics are (a) its subdivision into five major tectonic domains, bounded mostly by long shear zones, as follows: Medio Coreau, Ceara Central, Rio Grande do Norte, Transversal, and Southern; (b) the alternation of supracrustal belts with reworked basement inliers (Archean nuclei + Paleoproterozoic belts); and (c) the diversity of granitic plutonism, from Neoproterozoic to Early Cambrian ages, that affect supracrustal rocks as well as basement inliers. Recently, orogenic rock assemblages of early Tonian (1000-920 Ma) orogenic evolution have been recognized, which are restricted to the Transversal and Southern domains of the Province. Within the Transversal Zone, the Alto Pajeu terrane locally includes some remnants of oceanic crust along with island arc and continental arc rock assemblages, but the dominant supracrustal rocks are mature and immature pelitic metasedimentary and metavolcaniclastic rocks. Contiguous and parallel to the Alto Pajeu terrane, the Riacho Gravata subterrane consists mainly of low-grade metamorphic successions of metarhythmites, some of which are clearly turbiditic in origin, metaconglomerates, and sporadic marbles, along with interbedded metarhyolitic and metadacitic volcanic or metavolcaniclastic rocks. Both terrane and subterrane are cut by syn-contractional intrusive sheets of dominantly peraluminous high-K calc-alkaline, granititic to granodioritic metaplutonic rocks. The geochemical patterns of both supracrustal and intrusive rocks show similarities with associations of mature continental arc volcano-sedimentary sequences, but some subordinate intra-plate characteristics are also found. In both the Alto Pajeu and Riacho Gravata terranes, TIMS and SHRIMP U-Pb isotopic data from zircons from both metavolcanic and metaplutonic rocks yield ages between 1.0 and 0.92 Ga, which define the time span for an event of orogenic character, the Cariris Velhos event. Less extensive occurrences of rocks of Cariris Velhos age are recognized mainly in the southernmost domains of the Province, as for example in the Polo Redondo-Maranco terrane, where arc-affinity migmatite-granitic and meta-volcano-sedimentary rocks show U-Pb ages (SHRIMP data) around 0.98-0.97 Ga. For all these domains, Sm-Nd data exhibit Tom model ages between 1.9 and 1.1 Ga with corresponding slightly negative to slightly positive epsilon(Nd)(t) values. These domains, along with the Borborema Province as a whole, were significantly affected by tectonic and magmatic events of the Brasiliano Cycle (0.7-0.5 Ga), so that it is possible that there are some other early Tonian rock assemblages which were completely masked and hidden by these later Brasiliano events. Cariris Velhos processes are younger than the majority of orogenic systems at the end of Mesoproterozoic Era and beginning of Neoproterozoic throughout the world, e.g. Irumide belt, Kibaride belt and Namaqua-Natal belt, and considerably younger than those of the youngest orogenic process (Ottawan) in the Grenvillian System. Therefore, they were probably not associated with the proposed assembly of Rodinia. We suggest, instead, that Cariris Velhos magmatism and tectonism could have been related to a continental margin magmatic arc, with possible back-arc associations, and that this margin may have been a short-lived (<100 m.y.) leading edge of the newly assembled Rodinia supercontinent. (C) 2009 Elsevier Ltd. All rights reserved.

Closing the Clymene ocean and bending a Brasiliano belt: Evidence for the Cambrian formation of Gondwana, southeast Amazon craton

We report new paleomagnetic and geochronological data from Ediacaran rift-drift carbonates in the Paraguai belt at the southern end of the suture zone between the Amazon craton and the Sao Francisco and Rio de Plata cratons, South America. Early thrusting resulted in remagnetization ca. 528 +/- 36 Ma or later; the mean age is established by (40)Ar/(39)Ar encapsulation dating of mixed authigenic and detrital illite from remagnetized carbonates from the unmetamorphosed fold-thrust belt. This remagnetization overlaps with a 525 Ma Gondwana reference pole. Metamorphic illite from the slate belt yields (40)Ar/(39)Ar ages of 496-484 Ma, the timing of peak regional metamorphism. Oroclinal bending of the Paraguai belt was caused by a 90 degrees clockwise rotation of the east-west limb after ca. 528 Ma, probably reflecting the irregular margin of the southeast Amazon craton. The age of the Paraguai belt overlaps with that of the Pampean orogeny farther south along the western margin of the Rio de Plata craton, suggesting a coeval closure for the Clymene ocean separating the Amazon craton from the Sao Francisco and Rio de Plata cratons.

Timing and structure of the 8.2 kyr BP event inferred from delta(18)O records of stalagmites from China, Oman, and Brazil

Oxygen isotope records of stalagmites from China and Oman reveal a weak summer monsoon event, with a double-plunging structure, that started 8.21 +/- 0.02 kyr B. P. An identical but antiphased pattern is also evident in two stalagmite records from eastern Brazil, indicating that the South American Summer Monsoon was intensified during the 8.2 kyr B. P. event. These records demonstrate that the event was of global extent and synchronous within dating errors of <50 years. In comparison with recent model simulations, it is plausible that the 8.2 kyr B. P. event can be tied in changes of the Atlantic Meridional Overturning Circulation triggered by a glacial lake draining event. This, in turn, affected North Atlantic climate and latitudinal position of the Intertropical Convergence Zone, resulting in the observed low-latitude monsoonal precipitation patterns.

Global perturbation of the marine Ca isotopic composition in the aftermath of the Marinoan global glaciation

Ca isotopic compositions of Marinoan post-glacial carbonate successions in Brazil and NW Canada were measured Both basal dolostones display delta(44/40)Ca values between 1 and 0 7 parts per thousand overlying limestones show a negative Ca isotope excursion to values around 0 1 parts per thousand and delta(44/40)Ca values rapidly increase up-section to near 2 0 parts per thousand In the Brazilian successions those high delta(44/40)Ca values rapidly decrease and stabilize to values between 0 6 and 0 9 parts per thousand These Ca isotope secular variation trends are unlike those of Sturtian post-glacial carbonate successions but similar to those of Marinoan post-glacial carbonate successions in Namibia suggesting that the perturbation of the marine Ca cycle was global This recommends Ca isotope stratigraphy as a tool to correlate Neoproterozoic post-glacial carbonate successions worldwide While the lowermost and uppermost strata have delta(44/40)Ca values typical of Phanerozoic carbonates the extremes 0 1 and 2 0 parts per thousand have not been thus far reported for other marine carbonates These extreme values suggest a short-lived non-actualistic perturbation in the marine Ca cycle Simple box modelling of the Marinoan post-glacial marine Ca cycle can reproduce the extreme values only by postulating a two-step process with Ca input initially exceeding Ca removal trough carbonate precipitation followed by precipitation overtaking a decreased Ca Input (C) 2010 Elsevier B V All rights reserved

Global minimization using an Augmented Lagrangian method with variable lower-level constraints

A novel global optimization method based on an Augmented Lagrangian framework is introduced for continuous constrained nonlinear optimization problems. At each outer iteration k the method requires the epsilon(k)-global minimization of the Augmented Lagrangian with simple constraints, where epsilon(k) -> epsilon. Global convergence to an epsilon-global minimizer of the original problem is proved. The subproblems are solved using the alpha BB method. Numerical experiments are presented.

Intrinsically Multivariate Predictive Genes

Canalizing genes possess such broad regulatory power, and their action sweeps across a such a wide swath of processes that the full set of affected genes are not highly correlated under normal conditions. When not active, the controlling gene will not be predictable to any significant degree by its subject genes, either alone or in groups, since their behavior will be highly varied relative to the inactive controlling gene. When the controlling gene is active, its behavior is not well predicted by any one of its targets, but can be very well predicted by groups of genes under its control. To investigate this question, we introduce in this paper the concept of intrinsically multivariate predictive (IMP) genes, and present a mathematical study of IMP in the context of binary genes with respect to the coefficient of determination (CoD), which measures the predictive power of a set of genes with respect to a target gene. A set of predictor genes is said to be IMP for a target gene if all properly contained subsets of the predictor set are bad predictors of the target but the full predictor set predicts the target with great accuracy. We show that logic of prediction, predictive power, covariance between predictors, and the entropy of the joint probability distribution of the predictors jointly affect the appearance of IMP genes. In particular, we show that high-predictive power, small covariance among predictors, a large entropy of the joint probability distribution of predictors, and certain logics, such as XOR in the 2-predictor case, are factors that favor the appearance of IMP. The IMP concept is applied to characterize the behavior of the gene DUSP1, which exhibits control over a central, process-integrating signaling pathway, thereby providing preliminary evidence that IMP can be used as a criterion for discovery of canalizing genes.

Bacterial Community Composition in Brazilian Anthrosols and Adjacent Soils Characterized Using Culturing and Molecular Identification

Microbial community composition was examined in two soil types, Anthrosols and adjacent soils, sampled from three locations in the Brazilian Amazon. The Anthrosols, also known as Amazonian dark earths, are highly fertile soils that are a legacy of pre-Columbian settlement. Both Anthrosols and adjacent soils are derived from the same parent material and subject to the same environmental conditions, including rainfall and temperature; however, the Anthrosols contain high levels of charcoal-like black carbon from which they derive their dark color. The Anthrosols typically have higher cation exchange capacity, higher pH, and higher phosphorus and calcium contents. We used culture media prepared from soil extracts to isolate bacteria unique to the two soil types and then sequenced their 16S rRNA genes to determine their phylogenetic placement. Higher numbers of culturable bacteria, by over two orders of magnitude at the deepest sampling depths, were counted in the Anthrosols. Sequences of bacteria isolated on soil extract media yielded five possible new bacterial families. Also, a higher number of families in the bacteria were represented by isolates from the deeper soil depths in the Anthrosols. Higher bacterial populations and a greater diversity of isolates were found in all of the Anthrosols, to a depth of up to 1 m, compared to adjacent soils located within 50-500 m of their associated Anthrosols. Compared to standard culture media, soil extract media revealed diverse soil microbial populations adapted to the unique biochemistry and physiological ecology of these Anthrosols.

Efficient solutions to robust, semi-implicit discretizations of the immersed boundary method

The immersed boundary method is a versatile tool for the investigation of flow-structure interaction. In a large number of applications, the immersed boundaries or structures are very stiff and strong tangential forces on these interfaces induce a well-known, severe time-step restriction for explicit discretizations. This excessive stability constraint can be removed with fully implicit or suitable semi-implicit schemes but at a seemingly prohibitive computational cost. While economical alternatives have been proposed recently for some special cases, there is a practical need for a computationally efficient approach that can be applied more broadly. In this context, we revisit a robust semi-implicit discretization introduced by Peskin in the late 1970s which has received renewed attention recently. This discretization, in which the spreading and interpolation operators are lagged. leads to a linear system of equations for the inter-face configuration at the future time, when the interfacial force is linear. However, this linear system is large and dense and thus it is challenging to streamline its solution. Moreover, while the same linear system or one of similar structure could potentially be used in Newton-type iterations, nonlinear and highly stiff immersed structures pose additional challenges to iterative methods. In this work, we address these problems and propose cost-effective computational strategies for solving Peskin`s lagged-operators type of discretization. We do this by first constructing a sufficiently accurate approximation to the system`s matrix and we obtain a rigorous estimate for this approximation. This matrix is expeditiously computed by using a combination of pre-calculated values and interpolation. The availability of a matrix allows for more efficient matrix-vector products and facilitates the design of effective iterative schemes. We propose efficient iterative approaches to deal with both linear and nonlinear interfacial forces and simple or complex immersed structures with tethered or untethered points. One of these iterative approaches employs a splitting in which we first solve a linear problem for the interfacial force and then we use a nonlinear iteration to find the interface configuration corresponding to this force. We demonstrate that the proposed approach is several orders of magnitude more efficient than the standard explicit method. In addition to considering the standard elliptical drop test case, we show both the robustness and efficacy of the proposed methodology with a 2D model of a heart valve. (C) 2009 Elsevier Inc. All rights reserved.

A Robust, Fully Adaptive Hybrid Level-Set/Front-Tracking Method for Two-Phase Flows with an Accurate Surface Tension Computation

We present a variable time step, fully adaptive in space, hybrid method for the accurate simulation of incompressible two-phase flows in the presence of surface tension in two dimensions. The method is based on the hybrid level set/front-tracking approach proposed in [H. D. Ceniceros and A. M. Roma, J. Comput. Phys., 205, 391400, 2005]. Geometric, interfacial quantities are computed from front-tracking via the immersed-boundary setting while the signed distance (level set) function, which is evaluated fast and to machine precision, is used as a fluid indicator. The surface tension force is obtained by employing the mixed Eulerian/Lagrangian representation introduced in [S. Shin, S. I. Abdel-Khalik, V. Daru and D. Juric, J. Comput. Phys., 203, 493-516, 2005] whose success for greatly reducing parasitic currents has been demonstrated. The use of our accurate fluid indicator together with effective Lagrangian marker control enhance this parasitic current reduction by several orders of magnitude. To resolve accurately and efficiently sharp gradients and salient flow features we employ dynamic, adaptive mesh refinements. This spatial adaption is used in concert with a dynamic control of the distribution of the Lagrangian nodes along the fluid interface and a variable time step, linearly implicit time integration scheme. We present numerical examples designed to test the capabilities and performance of the proposed approach as well as three applications: the long-time evolution of a fluid interface undergoing Rayleigh-Taylor instability, an example of bubble ascending dynamics, and a drop impacting on a free interface whose dynamics we compare with both existing numerical and experimental data.