The growth, lifetime, number density, and size of water droplets in warm atmospheric clouds determine the evolution, lifetime and light transmission properties of those clouds. These small-scale cloud properties, in addition to precipitation initiation, have strong implications for the Earth's energy budget since warm clouds cover large geographic areas. Spatio-temporal correlations on the millimetre scale and smaller may or may not affect these properties of clouds. To date, the pioneering measurements of such correlations in marine stratocumulus clouds have relied on averaging over holographically reconstructed volumes spanning at least ten kilometres. These have revealed weak but widespread spatial clustering of cloud droplets. Here we present results of strong localised clustering on scales of half a metre or less from holographic measurements collected with the Max Planck CloudKite in shallow cumulus clouds in the mid-Atlantic trade wind region near Barbados, with a spatial separation of only 12~cm between measurement volumes. This observation challenges the foundations of our understanding of cloud microphysics at the droplet scale, with implications for cloud modelling in weather and climate prediction.

We present a machine learning system that can quantify fine art paintings with a set of visual elements and principles of art. This formal analysis is fundamental for understanding art, but developing such a system is challenging. Paintings have high visual complexities, but it is also difficult to collect enough training data with direct labels. To resolve these practical limitations, we introduce a novel mechanism, called proxy learning, which learns visual concepts in paintings though their general relation to styles. This framework does not require any visual annotation, but only uses style labels and a general relationship between visual concepts and style. In this paper, we propose a novel proxy model and reformulate four pre-existing methods in the context of proxy learning. Through quantitative and qualitative comparison, we evaluate these methods and compare their effectiveness in quantifying the artistic visual concepts, where the general relationship is estimated by language models; GloVe or BERT. The language modeling is a practical and scalable solution requiring no labeling, but it is inevitably imperfect. We demonstrate how the new proxy model is robust to the imperfection, while the other models are sensitively affected by it.

We present a machine learning system that can quantify fine art paintings with a set of visual elements and principles of art. This formal analysis is fundamental for understanding art, but developing such a system is challenging. Paintings have high visual complexities, but it is also difficult to collect enough training data with direct labels. To resolve these practical limitations, we introduce a novel mechanism, called proxy learning, which learns visual concepts in paintings though their general relation to styles. This framework does not require any visual annotation, but only uses style labels and a general relationship between visual concepts and style. In this paper, we propose a novel proxy model and reformulate four pre-existing methods in the context of proxy learning. Through quantitative and qualitative comparison, we evaluate these methods and compare their effectiveness in quantifying the artistic visual concepts, where the general relationship is estimated by language models; GloVe or BERT. The language modeling is a practical and scalable solution requiring no labeling, but it is inevitably imperfect. We demonstrate how the new proxy model is robust to the imperfection, while the other models are sensitively affected by it.

We show that discrete $W_m$ lattices are bi-Hamiltonian, using geometric realizations of discretizations of the Adler-Gel'fand-Dikii flows as local evolutions of arc length-parametrized polygons in centro-affine space. We prove the compatibility of two known Hamiltonian structure defined on the space of geometric invariants by lifting them to a pair of pre-symplectic forms on the space of arc length parametrized polygons. The simplicity of the expressions of the pre-symplectic forms makes the proof of compatibility straightforward. We also study their kernels and possible integrable systems associated to the pair.

We present new high-resolution ($\sim$0\farcs09) $H$-band imaging observations of the circumstellar disk around the T Tauri star SU Aur. Our observations with Subaru-HiCIAO have revealed the presence of scattered light as close as 0\farcs15 ($\sim$20 AU) to the star. Within our image, we identify bright emission associated with a disk with a minimum radius of $\sim$90 AU, an inclination of $\sim$35\degr from the plane of the sky, and an approximate P.A. of 15\degr for the major axis. We find a brightness asymmetry between the northern and southern sides of the disk due to a non-axisymmetric disk structure. We also identify a pair of asymmetric tail structures extending east and west from the disk. The western tail extends at least 2\farcs5 (350 AU) from the star, and is probably associated with a reflection nebula previously observed at optical and near-IR wavelengths. The eastern tail extends at least 1\arcsec (140 AU) at the present signal-to-noise. These tails are likely due to an encounter with an unseen brown dwarf, but our results do not exclude the explanation that these tails are outflow cavities or jets.

We compare X-ray emission from several general relativistic, multi-frequency, radiation magnetohydrodynamic simulations of thin black hole accretion disks with different accretion rates and spins. The simulations were performed using the M1 closure scheme, resolved with twelve frequency (energy) bins logarithmically spaced from $5 \times 10^{-3}$ to $5 \times 10^3$ keV. We apply a general relativistic Monte Carlo transport code to post-process the simulation data with greater fidelity in frequency resolution and Compton scattering treatment. Despite the relatively few energy bins and Kompaneets approximation to Compton scattering utilized in the M1 method, we find generally good agreement between the methods. Both produce prominent thermal profiles with peaks around 2 - 2.5 keV, where agreement is particularly strong and representative of the soft state. Both also find weaker (lower luminosity) thermally sourced emission extending out to 100 keV due to the hotter innermost regions of the disks. Inverse Compton scattering becomes increasingly effective at hardening spectral outputs with increasing black hole spin, and becomes the dominant mechanism for photons that escape with energies between 10 to several hundred keV. At very high rates of spin the radiation flux in this upscattered component becomes comparable to the thermal flux, a phenomenon typically associated with intermediate states. Beyond $10^4$ keV, we observe faint, free-free emission from hot, optically thin coronal regions developing near the horizon, common to both spinning and nonspinning black holes.

We present high resolution H-band polarized intensity (PI; FWHM = 0."1: 14 AU) and L'-band imaging data (FWHM = 0."11: 15 AU) of the circumstellar disk around the weak-lined T Tauri star PDS 70 in Centaurus at a radial distance of 28 AU (0."2) up to 210 AU (1."5). In both images, a giant inner gap is clearly resolved for the first time, and the radius of the gap is ~70 AU. Our data show that the geometric center of the disk shifts by ~6 AU toward the minor axis. We confirm that the brown dwarf companion candidate to the north of PDS 70 is a background star based on its proper motion. As a result of SED fitting by Monte Carlo radiative transfer modeling, we infer the existence of an optically thick inner disk at a few AU. Combining our observations and modeling, we classify the disk of PDS 70 as a pre-transitional disk. Furthermore, based on the analysis of L'-band imaging data, we put an upper limit mass of companions at ~30 to ~50MJ within the gap. Taking account of the presence of the large and sharp gap, we suggest that the gap could be formed by dynamical interactions of sub-stellar companions or multiple unseen giant planets in the gap.

Much of the baryonic matter in the Universe, including the most active and luminous sources, are best studied in the X-ray band. Key advances in X-ray optics and detectors have paved the way for the Advanced X-ray Imaging Satellite (AXIS), a Probe-class mission that is a major improvement over Chandra, which has generated a steady stream of important discoveries for the past 2 decades. AXIS can be launched in the late 2020s and will transform our understanding in several major areas of astrophysics, including the growth and fueling of supermassive black holes, galaxy formation and evolution, the microphysics of cosmic plasmas, the time-variable universe, and a wide variety of cutting-edge studies. Relative to Chandra, the AXIS PSF is nearly twice as sharp on-axis; its field of view for subarcsecond imaging 70 times larger by area; its effective area at 1 keV is 10 times larger. The low-Earth orbit ensures a low and stable detector background, resulting in 50 times greater sensitivity than Chandra for extended sources. AXIS has a rapid repointing response with operations similar to Swift, but is 100 times more sensitive for time-domain science. These capabilities open up a vast discovery space and complement the next generation of astronomical observatories. A high-spectral-resolution mission (Athena) operating at the same time as a high-angular-resolution mission (AXIS) greatly increases the range of scientific discovery. AXIS will use lightweight X-ray optics made of thin single-crystal silicon mirrors developed at NASA Goddard. The detector array builds on a long legacy of X-ray CCD and provides improved photon localization, much faster readout time, and broader energy band. The estimated mission costs are consistent with the $1B Probe mission cost guideline.

This paper presents a continuation of our efforts to numerically study accretion disks that are misaligned (tilted) with respect to the rotation axis of a Kerr black hole. Here we present results of a global numerical simulation which fully incorporates the effects of the black hole spacetime as well as magnetorotational turbulence that is the primary source of angular momentum transport in the flow. This simulation shows dramatic differences from comparable simulations of untilted disks. Accretion onto the hole occurs predominantly through two opposing plunging streams that start from high latitudes with respect to both the black-hole and disk midplanes. This is due to the aspherical nature of the gravitational spacetime around the rotating black hole. These plunging streams start from a larger radius than would be expected for an untilted disk. In this regard the tilted black hole effectively acts like an untilted black hole of lesser spin. Throughout the duration of the simulation, the main body of the disk remains tilted with respect to the symmetry plane of the black hole; thus there is no indication of a Bardeen-Petterson effect in the disk at large. The torque of the black hole instead principally causes a global precession of the main disk body. In this simulation the precession has a frequency of $3 (M_\odot/M)$ Hz, a value consistent with many observed low-frequency quasi-periodic oscillations. However, this value is strongly dependent on the size of the disk, so this frequency may be expected to vary over a large range.

The problem of interpreting or aggregating multiple rankings is common to many real-world applications. Perhaps the simplest and most common approach is a weighted rank aggregation, wherein a (convex) weight is applied to each input ranking and then ordered. This paper describes a new tool for visualizing and displaying ranking information for the weighted rank aggregation method. Traditionally, the aim of rank aggregation is to summarize the information from the input rankings and provide one final ranking that hopefully represents a more accurate or truthful result than any one input ranking. While such an aggregated ranking is, and clearly has been, useful to many applications, it also obscures information. In this paper, we show the wealth of information that is available for the weighted rank aggregation problem due to its structure. We apply weight set decomposition to the set of convex multipliers, study the properties useful for understanding this decomposition, and visualize the indifference regions. This methodology reveals information--that is otherwise collapsed by the aggregated ranking--into a useful, interpretable, and intuitive decision support tool. Included are multiple illustrative examples, along with heuristic and exact algorithms for computing the weight set decomposition.

Matter flows in the central regions of quasars during their active phases are probably responsible for the properties of the super-massive black holes and that of the bulges of host galaxies. To understand how this mechanism works, we need to characterize the geometry and the physical state of the accreting matter at cosmological redshifts. The few high quality X-ray spectra of distant QSO have been collected by adding sparse pointings of single objects obtained during X-ray monitoring campaigns. This could have introduced spurious spectral features due to source variability. Here we present a single epoch, high-quality X-ray spectrum of the z=3.62 quasar B1422+231 whose flux is enhanced by gravitationally lensing (F$_{2-10 keV}\sim$10$^{-12}$erg s$^{-1}$ cm$^{-2}$). The X-ray spectrum of B1422+231 is found to be very similar to the one of a typical nearby Seyfert galaxy. Neutral absorption is detected (N$_{H}\sim$5$\times$10$^{21}$ cm$^{-2}$ at the redshift of the source) while a strong absorption edge is measured at E$\sim$7.5 keV with an optical depth of $\tau\sim$0.14. We also find hints of the FeK$\alpha$ line in emission at E$\sim$6.4 keV line (EW$\lesssim$70 eV) and a hump is detected in the E$\sim$15-20 keV energy band (rest-frame) suggesting the presence of a reflection component. In this scenario, the primary emission of B1422+231 is most probably dominated by the thermal Comptonization of UV seed photons in a corona with kT$\sim$40 keV and the reflection component has a relative direct-to-reflect normalization r$\sim$1. These findings confirm that gravitational lensing is effective to obtain good quality X-ray spectral information of quasar at high-z, moreover they support the idea that the same general picture characterizing active galactic nuclei in the nearby Universe is valid also at high redshift.

Dust trapping accelerates the coagulation of dust particles, and thus it represents an initial step toward the formation of planetesimals. We report $H$-band (1.6 um) linear polarimetric observations and 0.87 mm interferometric continuum observations toward a transitional disk around LkHa 330. As results, a pair of spiral arms were detected in the $H$-band emission and an asymmetric (potentially arm-like) structure was detected in the 0.87 mm continuum emission. We discuss the origin of the spiral arm and the asymmetric structure, and suggest that a massive unseen planet is the most plausible explanation. The possibility of dust trapping and grain growth causing the asymmetric structure was also investigated through the opacity index (beta) by plotting the observed SED slope between 0.87 mm from our SMA observation and 1.3 mm from literature. The results imply that grains are indistinguishable from ISM-like dust in the east side ($beta = 2.0 pm 0.5$), but much smaller in the west side $beta = 0.7^{+0.5}_{-0.4}$, indicating differential dust size distribution between the two sides of the disk. Combining the results of near-infrared and submillimeter observations, we conjecture that the spiral arms exist at the upper surface and an asymmetric structure resides in the disk interior. Future observations at centimeter wavelengths and differential polarization imaging in other bands (Y to K) with extreme AO imagers are required to understand how large dust grains form and to further explore the dust distribution in the disk.

We consider a nonlinear (Berger or Von Karman) clamped plate model with a {\em piston-theoretic} right hand side---which include non-dissipative, non-conservative lower order terms. The model arises in aeroelasticity when a panel is immersed in a high velocity linear potential flow; in this case the effect of the flow can be captured by a dynamic pressure term written in terms of the material derivative of the plate's displacement. The effect of fully-supported internal damping is studied for both Berger and von Karman dynamics. The non-dissipative nature of the dynamics preclude the use of strong tools such as backward-in-time smallness of velocities and finiteness of the dissipation integral. Modern quasi-stability techniques are utilized to show the existence of compact global attractors and generalized fractal exponential attractors. Specific results depending on the size of the damping parameter and the nonlinearity in force. For the Berger plate, in the presence of large damping, the existence of a proper global attractor (whose fractal dimension is finite in the state space) is shown via a decomposition of the nonlinear dynamics. This leads to the construction of a compact set upon which quasi-stability theory can be implemented. Numerical investigations for appropriate 1-D models are presented which explore and support the abstract results presented herein.

Starshade is one of the technologies that will enable the observation and characterization of small planets around nearby stars through direct imaging. The Starshade Exoplanetary Data Challenge (SEDC) was designed to validate starshade-imaging's noise budget and evaluate the capabilities of image-processing techniques, by inviting community participating teams to analyze >1000 simulated images of hypothetical exoplanetary systems observed through a starshade. Because the starshade would suppress the starlight so well, the dominant noise source and the main challenge for the planet detection becomes the exozodiacal disks and their structures. In this paper, we summarize the techniques used by the participating teams and compare their findings with the truth. With an independent component analysis to remove the background, about 70% of the inner planets (close to the inner working angle) have been detected and ~40% of the outer planet (fainter than the inner counterparts) have been identified. Planet detection becomes more difficult in the cases of higher disk inclination, as the false negative and false positive counts increase. Interestingly, we found little difference in the planet detection ability between 1e-10 and 1e-9 instrument contrast, confirming that the dominant limitations are from the astrophysical background and not due to the performance of the starshade. Finally, we find that a non-parametric background calibration scheme, such as the independent component analysis reported here, results in a mean residual of 10% the background brightness. This background estimation error leads to substantial false positives and negatives and systematic bias in the planet flux estimation, and should be included in the estimation of the planet detection signal-to-noise ratio for imaging using a starshade and also a coronagraph that delivers exozodi-limited imaging.

The Fermi GBM Catalog has been recently published. Previous classification analyses of the BATSE, RHESSI, BeppoSAX, and Swift databases found three types of gamma-ray bursts. Now we analyzed the GBM catalog to classify the GRBs. PCA and Multiclustering analysis revealed three groups. Validation of these groups, in terms of the observed variables, shows that one of the groups coincides with the short GRBs. The other two groups split the long class into a bright and dim part, as defined by the peak flux. Additional analysis is needed to determine whether this splitting is only a mathematical byproduct of the analysis or has some real physical meaning.

The immense tidal forces of massive black holes can rip apart stars that come too close to them. As the resulting stellar debris spirals inwards, it heats up and emits x-rays when near the black hole. Here, we report the discovery of an exceptionally stable 131-second x-ray quasi-periodicity from a black hole after it disrupted a star. Using a black hole mass indicated from host galaxy scaling relations implies that, (1) this periodicity originates from very close to the black hole's event horizon, and (2) the black hole is rapidly spinning. Our findings suggest that other disruption events with similar highly sensitive observations likely also exhibit quasi-periodicities that encode information about the fundamental properties of their black holes.

We present a new Subaru/HiCIAO high-contrast H-band polarized intensity (PI) image of a nearby transitional disk associated with TW Hydrae. The scattered light from the disk was detected from 0.2" to 1.5" (11 - 81 AU) and the PI image shows a clear axisymmetric depression in polarized intensity at ~ 0.4" (~ 20 AU) from the central star, similar to the ~ 80 AU gap previously reported from HST images. Azimuthal polarized intensity profile also shows the disk beyond 0.2" is almost axisymmetric. We discuss two possible scenarios explaining the origin of the polarized intensity depression: 1) a gap structure may exist at ~ 20 AU from the central star because of shallow slope seen in the polarized intensity profile, and 2) grain growth may be occurring in the inner region of the disk. Multi-band observations at NIR and millimeter/sub-millimeter wavelengths play a complementary role in investigating dust opacity and may help reveal the origin of the gap more precisely.

We analyze the two brightest Chandra X-ray flares detected from Sagittarius A*, with peak luminosities more than 600 x and 245 x greater than the quiescent X-ray emission. The brightest flare has a distinctive double-peaked morphology --- it lasts 5.7 ksec ($\sim 2$ hours), with a rapid rise time of 1500 sec and a decay time of 2500 sec. The second flare lasts 3.4 ksec, with rise and decay times of 1700 sec and 1400 sec. These luminous flares are significantly harder than quiescence: the first has a power law spectral index $\Gamma = 2.06\pm 0.14$ and the second has $\Gamma = 2.03\pm 0.27$, compared to $\Gamma = 3.0\pm0.2$ for the quiescent accretion flow. These spectral indices (as well as the flare hardness ratios) are consistent with previously-detected Sgr A* flares, suggesting that bright and faint flares arise from similar physical processes. Leveraging the brightest flare's long duration and high signal-to-noise, we search for intraflare variability and detect excess X-ray power at a frequency of $\nu \approx 3$ mHz, but show that it is an instrumental artifact and not of astrophysical origin. We find no other evidence (at the 95% confidence level) for periodic or quasi-periodic variability in either flares' time series. We also search for non-periodic excess power but do not find compelling evidence in the power spectrum. Bright flares like these remain our most promising avenue for identifying Sgr A*'s short timescale variability in the X-ray, which may probe the characteristic size scale for the X-ray emission region.

We present a statistical analysis of the X-ray flux distribution of Sgr A* from the Chandra X-ray Observatory's 3 Ms Sgr A* X-ray Visionary Project (XVP) in 2012. Our analysis indicates that the observed X-ray flux distribution can be decomposed into a steady quiescent component, represented by a Poisson process with rate $Q=(5.24\pm0.08)\times10^{-3}$ cts s$^{-1},$ and a variable component, represented by a power law process ($dN/dF\propto F^{-\xi},$ $\xi=1.92_{-0.02}^{+0.03}$). This slope matches our recently-reported distribution of flare luminosities. The variability may also be described by a log-normal process with a median unabsorbed 2-8 keV flux of $1.8^{+0.9}_{-0.6}\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$ and a shape parameter $\sigma=2.4\pm0.2,$ but the power law provides a superior description of the data. In this decomposition of the flux distribution, all of the intrinsic X-ray variability of Sgr A* (spanning at least three orders of magnitude in flux) can be attributed to flaring activity, likely in the inner accretion flow. We confirm that at the faint end, the variable component contributes ~10% of the apparent quiescent flux, as previously indicated by our statistical analysis of X-ray flares in these Chandra observations. Our flux distribution provides a new and important observational constraint on theoretical models of Sgr A*, and we use simple radiation models to explore the extent to which a statistical comparison of the X-ray and infrared can provide insights into the physics of the X-ray emission mechanism.