GRMHD simulations of accretion onto Sgr A*: How important are radiative losses?

We present general relativistic magnetohydrodynamic (GRMHD) numerical simulations of the accretion flow around the supermassive black hole in the Galactic centre, Sagittarius A* (Sgr A*). The simulations include for the first time radiative cooling processes (synchrotron, bremsstrahlung, and inverse Compton) self-consistently in the dynamics, allowing us to test the common simplification of ignoring all cooling losses in the modeling of Sgr A*. We confirm that for Sgr A*, neglecting the cooling losses is a reasonable approximation if the Galactic centre is accreting below ~10^{-8} Msun/yr i.e. Mdot < 10^{-7} Mdot_Edd. But above this limit, we show that radiative losses should be taken into account as significant differences appear in the dynamics and the resulting spectra when comparing simulations with and without cooling. This limit implies that most nearby low-luminosity active galactic nuclei are in the regime where cooling should be taken into account. We further make a parameter study of axisymmetric gas accretion around the supermassive black hole at the Galactic centre. This approach allows us to investigate the physics of gas accretion in general, while confronting our results with the well studied and observed source, Sgr A*, as a test case. We confirm that the nature of the accretion flow and outflow is strongly dependent on the initial geometry of the magnetic field. For example, we find it difficult, even with very high spins, to generate powerful outflows from discs threaded with multiple, separate poloidal field loops.Comment: Resubmitted to MNRAS, including modifications in response to referee report. 13 pages, 15 figure

Measurement and analysis of critical crack tip processes during fatigue crack growth

The mechanics of fatigue crack growth under constant-amplitudes and variable-amplitude loading were examined. Critical loading histories involving relatively simple overload and overload/underload cycles were studied to provide a basic understanding of the underlying physical processes controlling crack growth. The material used for this study was 7091-T7E69, a powder metallurgy aluminum alloy. Local crack-tip parameters were measured at various times before, during, and after the overloads, these include crack-tip opening loads and displacements, and crack-tip strain fields. The latter were useed, in combination with the materials cyclic and monotonic stress-strain properties, to compute crack-tip residual stresses. The experimental results are also compared with analytical predictions obtained using the FAST-2 computer code. The sensitivity of the analytical model to constant-amplitude fatigue crack growth rate properties and to through-thickness constrain are studied

Radiative Models of Sagittarius A* and M87 from Relativistic MHD Simulations

Ongoing millimeter VLBI observations with the Event Horizon Telescope allow unprecedented study of the innermost portion of black hole accretion flows. Interpreting the observations requires relativistic, time-dependent physical modeling. We discuss the comparison of radiative transfer calculations from general relativistic MHD simulations of Sagittarius A* and M87 with current and future mm-VLBI observations. This comparison allows estimates of the viewing geometry and physical conditions of the Sgr A* accretion flow. The viewing geometry for M87 is already constrained from observations of its large-scale jet, but, unlike Sgr A*, there is no consensus for its millimeter emission geometry or electron population. Despite this uncertainty, as long as the emission region is compact, robust predictions for the size of its jet launching region can be made. For both sources, the black hole shadow may be detected with future observations including ALMA and/or the LMT, which would constitute the first direct evidence for a black hole event horizon.Comment: 8 pages, 2 figures, submitted to the proceedings of AHAR 2011: The Central Kiloparse

Analytical solutions of bound timelike geodesic orbits in Kerr spacetime

We derive the analytical solutions of the bound timelike geodesic orbits in Kerr spacetime. The analytical solutions are expressed in terms of the elliptic integrals using Mino time $\lambda$ as the independent variable. Mino time decouples the radial and polar motion of a particle and hence leads to forms more useful to estimate three fundamental frequencies, radial, polar and azimuthal motion, for the bound timelike geodesics in Kerr spacetime. This paper gives the first derivation of the analytical expressions of the fundamental frequencies. This paper also gives the first derivation of the analytical expressions of all coordinates for the bound timelike geodesics using Mino time. These analytical expressions should be useful not only to investigate physical properties of Kerr geodesics but more importantly to applications related to the estimation of gravitational waves from the extreme mass ratio inspirals.Comment: A typo in the first expression in equation 21 was fixe

Certification of Compiler Optimizations using Kleene Algebra with Tests

We use Kleene algebra with tests to verify a wide assortment of common compiler optimizations, including dead code elimination, common subexpression elimination, copy propagation, loop hoisting, induction variable elimination, instruction scheduling, algebraic simplification, loop unrolling, elimination of redundant instructions, array bounds check elimination, and introduction of sentinels. In each of these cases, we give a formal equational proof of the correctness of the optimizing transformation

Automatic Labeling of Special Diagnostic Mammography Views from Images and DICOM Headers.

Applying state-of-the-art machine learning techniques to medical images requires a thorough selection and normalization of input data. One of such steps in digital mammography screening for breast cancer is the labeling and removal of special diagnostic views, in which diagnostic tools or magnification are applied to assist in assessment of suspicious initial findings. As a common task in medical informatics is prediction of disease and its stage, these special diagnostic views, which are only enriched among the cohort of diseased cases, will bias machine learning disease predictions. In order to automate this process, here, we develop a machine learning pipeline that utilizes both DICOM headers and images to predict such views in an automatic manner, allowing for their removal and the generation of unbiased datasets. We achieve AUC of 99.72% in predicting special mammogram views when combining both types of models. Finally, we apply these models to clean up a dataset of about 772,000 images with expected sensitivity of 99.0%. The pipeline presented in this paper can be applied to other datasets to obtain high-quality image sets suitable to train algorithms for disease detection

Resonant Energy Exchange between Atoms in Dispersing and Absorbing Surroundings

Within the framework of quantization of the macroscopic electromagnetic field, a master equation describing both the resonant dipole-dipole interaction (RDDI) and the resonant atom-field interaction (RAFI) in the presence of dispersing and absorbing macroscopic bodies is derived, with the relevant couplings being expressed in terms of the surroundings-assisted Green tensor. It is shown that under certain conditions the RDDI can be regarded as being governed by an effective Hamiltonian. The theory, which applies to both weak and strong atom-field coupling, is used to study the resonant energy exchange between two (two-level) atoms sharing initially a single excitation. In particular, it is shown that in the regime of weak atom-field coupling there is a time window, where the energy transfer follows a transfer-rate law of the type obtained by ordinary second-order perturbation theory. Finally, the spectrum of the light emitted during the energy transfer is studied and the line splittings are discussed.Comment: 9 pages, 5 figs, Proceedings of ICQO'2002, Raubichi, to appear in Optics and Spectroscop

The Post-Pericenter Evolution of the Galactic Center Source G2

In early 2014 the fast-moving near-infrared source G2 reached its closest approach to the supermassive black hole Sgr A* in the Galactic Center. We report on the evolution of the ionized gaseous component and the dusty component of G2 immediately after this event, revealed by new observations obtained in 2015 and 2016 with the SINFONI integral field spectrograph and the NACO imager at the ESO VLT. The spatially resolved dynamics of the Br$\gamma$ line emission can be accounted for by the ballistic motion and tidal shearing of a test-particle cloud that has followed a highly eccentric Keplerian orbit around the black hole for the last 12 years. The non-detection of a drag force or any strong hydrodynamic interaction with the hot gas in the inner accretion zone limits the ambient density to less than a few 10$^3$ cm$^{-3}$ at the distance of closest approach (1500 $R_s$), assuming G2 is a spherical cloud moving through a stationary and homogeneous atmosphere. The dust continuum emission is unresolved in L'-band, but stays consistent with the location of the Br$\gamma$ emission. The total luminosity of the Br$\gamma$ and L' emission has remained constant to within the measurement uncertainty. The nature and origin of G2 are likely related to that of the precursor source G1, since their orbital evolution is similar, though not identical. Both object are also likely related to a trailing tail structure, which is continuously connected to G2 over a large range in position and radial velocity.Comment: 17 pages, 12 figures; accepted for publication in Ap