作者:Muhammad A. Latif, Daniel J. Whalen, Sadegh Khochfar
How quasars powered by supermassive black holes (SMBHs) formed less than a billion years after the Big Bang is still one of the outstanding problems in astrophysics 20 years after their discovery$^{1-4}$. Cosmological simulations suggest that rare cold flows converging on primordial haloes in low-shear environments could have created these quasars if they were 10$^4$ - 10$^5$ M$_{\odot}$ at birth but could not resolve their formation$^{5-8}$. Semianalytical studies of the progenitor halo of a primordial quasar found that it favours the formation of such seeds but could not verify if one actually appeared$^9$. Here we show that a halo at the rare convergence of strong, cold accretion flows creates massive BH seeds without the need for UV backgrounds, supersonic streaming motions, or even atomic cooling. Cold flows drive violent, supersonic turbulence in the halo that prevents star formation until it reaches a mass that triggers sudden, catastrophic baryon collapse that forms 31,000 and 40,000 M$_{\odot}$ stars. This simple, robust process ensures that haloes capable of forming quasars by z $>$ 6 produce massive seeds. The first quasars were thus a natural consequence of structure
作者:Ana López-Sepulcre, Mihkel Kama, Cecilia Ceccarelli
As part of the Herschel Guaranteed Time Key Programme CHESS, we present the discovery of a diffuse gas component in the foreground of the intermediate-mass protostar OMC-2 FIR 4, located in the Orion A region. Making use of the full HIFI spectrum of OMC-2 FIR 4 obtained in CHESS, we detected several ground-state lines from OH+, H2O+, HF, and CH+, all of them seen in absorption against the dust continuum emission of the protostar's envelope. The lines peak at a velocity of 9 km/s, which is blue-shifted by 2 km/s with respect to the systemic velocity of OMC-2 FIR 4 (V = 11.4 km/s). We derived column densities for each species, as well as an upper limit to the column density of the undetected H3O+. In order to model and characterise the foreground cloud, we used the Meudon PDR code to run a homogeneous grid of models that spans a reasonable range of densities, visual extinctions, cosmic ray ionisation rates and far-ultraviolet (FUV) radiation fields. The results of our modelling indicate that the foreground cloud is composed of predominantly neutral diffuse gas (n(H) = 100 cm-3) and is heavily irradiated by an external source of FUV that most likely arises from the nearby Trapezium OB
作者:Daniel Rahner, Eric W. Pellegrini, Simon C. O. Glover
The 30 Doradus Nebula in the Large Magellanic Cloud (LMC) contains the massive starburst cluster NGC 2070 with a massive and probably younger stellar sub clump at its center: R136. It is not clear how such a massive inner cluster could form several million years after the older stars in NGC 2070, given that stellar feedback is usually thought to expel gas and inhibit further star formation. Using the recently developed 1D feedback scheme WARPFIELD to scan a large range of cloud and cluster properties, we show that an age offset of several million years between the stellar populations is in fact to be expected given the interplay between feedback and gravity in a giant molecular cloud (GMC) with a density $\gtrsim 500$ cm$^{-3}$ due to re-accretion of gas onto the older stellar population. Neither capture of field stars nor gas retention inside the cluster have to be invoked in order to explain the observed age offset in NGC 2070 as well as the structure of the interstellar medium (ISM) around it.
The Traveling Salesman Problem (TSP) in the $d$-dimensional Euclidean space is among the oldest and most famous NP-hard optimization problems. In breakthrough works, Arora [J. ACM 1998] and Mitchell [SICOMP 1999] gave the first polynomial time approximation schemes. To improve the running time, Rao and Smith [STOC 1998] gave a randomized $(1/\varepsilon)^{O(1/\varepsilon^{d-1})}\cdot n\log n$ time approximation scheme. Bartal and Gottlieb [FOCS 2013] gave a randomized approximation scheme in $2^{(1/\varepsilon)^{O(d)}} n$ time, which is linear in $n$. Recently, Kisfaludi-Bak, Nederlof, and Węgrzycki [FOCS 2021] gave a randomized approximation scheme in $2^{O(1/\varepsilon^{d-1})} n \log n$ time, achieving a Gap-ETH tight dependence on $\varepsilon$. It is raised as a challenging open question by Kisfaludi-Bak, Nederlof, and Węgrzycki [FOCS 2021] whether a running time of $2^{O(1/\varepsilon^{d-1})}n$ is achievable. We answer their question positively by giving a randomized $2^{O(1/\varepsilon^{d-1})} n$ time approximation scheme for Euclidean TSP.
Randomized iterative methods, such as the Kaczmarz method and its variants, have gained growing attention due to their simplicity and efficiency in solving large-scale linear systems. Meanwhile, absolute value equations (AVE) have attracted increasing interest due to their connection with the linear complementarity problem. In this paper, we investigate the application of randomized iterative methods to generalized AVE (GAVE). Our approach differs from most existing works in that we tackle GAVE with non-square coefficient matrices. We establish more comprehensive sufficient and necessary conditions for characterizing the solvability of GAVE and propose precise error bound conditions. Furthermore, we introduce a flexible and efficient randomized iterative algorithmic framework for solving GAVE, which employs randomized sketching matrices drawn from user-specified distributions. This framework is capable of encompassing many well-known methods, including the Picard iteration method and the randomized Kaczmarz method. Leveraging our findings on solvability and error bounds, we establish both almost sure convergence and linear convergence rates for this versatile algorithmic framework.
作者:Jeremy S. C. Clark, Piotr Kulig, Konrad Podsiadlo
Bernstein fits implemented into R allow another route for Kruskal-Wallis power-study tool development. Monte-Carlo Kruskal-Wallis power studies were compared with measured power, with Monte-Carlo ANOVA equivalent and with an analytical method, with or without normalization, using four simulated runs each with 60-100 populations (each population with N=30000 from a set of Pearson-type ranges): random selection gave 6300 samples analysed for predictive power. Three medical-study datasets (Dialysis/systolic blood pressure; Diabetes/sleep-hours; Marital-status/high-density-lipoprotein cholesterol) were also analysed. In three from four simulated runs (run_one, run_one_relaxed, and run_three) with Pearson types pooled, Monte-Carlo Kruskal-Wallis gave predicted sample sizes significantly slightly lower than measured but more accurate than with ANOVA methods; the latter gave high sample-size predictions. Populations (run_one_relaxed) with ANOVA assumptions invalid gave Kruskal-Wallis predictions similar to those measured. In two from three medical studies, Kruskal-Wallis predictions (Dialysis: similar predictions; Marital: higher than measured) were more accurate than ANOVA (both higher t