The EUV provides most of the ionization that creates the high equivalent width (EW) broad and narrow emission lines (BELs, NELs) of quasars. Spectra of Hypermassive Schwarzschild black holes (HMBHs, $M_{BH} \geq 10^{10} M_{\odot}$) with $\alpha$-discs, decline rapidly in the EUV suggesting much lower EWs. Model spectra for black holes of mass $10^{6}-10^{12} M_{\odot}$ and accretion rates $0.03 \leq L_{bol}/L_{edd} \leq 1.0$ were input to the CLOUDY photoionization code. BELs become $\sim$100 times weaker in EW from $M_{BH} \sim 10^8 M_{\odot}$ to $M_{BH} \sim 10^{10} M_{\odot}$. The high ionization BELs (O VI 1034 $\overset{\circ}{\mathrm {A}}$, C IV 1549 $\overset{\circ}{\mathrm {A}}$, He II 1640 $\overset{\circ}{\mathrm {A}}$) decline in EW from ($M_{BH} \geq 10^6 M_{\odot}$, reproducing the Baldwin effect, but regain EW for $M_{BH} \geq 10^{10} M_{\odot}$). The low ionization lines (MgII 2798 $\overset{\circ}{\mathrm {A}}$, H$\beta$ 4861 $\overset{\circ}{\mathrm {A}}$ and H$\alpha$ 6563 $\overset{\circ}{\mathrm {A}}$) remain weak. Lines for maximally spinning HMBHs behave similarly. Line ratio diagrams for the BELs show that high OVI/H$\beta$ and low CIV/H$\alpha$ may pick out HMBH, although OVI is often hard to observe. In NEL BPT diagrams HMBHs lie among star-forming regions, except for highly spinning, high accretion rate HMBHs. In summary, the BELs expected from HMBHs would be hard to detect using the current optical facilities. From 100 to $10^{12} M_{\odot}$, the emission lines used to detect AGN only have high EW in the $10^6 - 10^9 M_{\odot}$ window, where most AGN are found. This selection effect may be distorting reported distributions of $M_{BH}$.

Rather general considerations of the string theory landscape imply a mild statistical draw towards large soft SUSY breaking terms tempered by the requirement of proper electroweak symmetry breaking where SUSY contributions to the weak scale are not too far from m(weak)~ 100 GeV. Such a picture leads to the prediction that m_h~ 125 GeV while most sparticles are beyond current LHC reach. Here we explore the possibility that the magnitude of the Peccei-Quinn (PQ) scale f_a is also set by string landscape considerations within the framework of a compelling SUSY axion model. First, we examine the case where the PQ symmetry arises as an accidental approximate global symmetry from a more fundamental gravity-safe Z(24)^R symmetry and where the SUSY mu parameter arises from a Kim-Nilles operator. The pull towards large soft terms then also pulls the PQ scale as large as possible. Unless this is tempered by rather severe (unknown) cosmological or anthropic bounds on the density of dark matter, then we would expect a far greater abundance of dark matter than is observed. This conclusion cannot be negated by adopting a tiny axion misalignment angle theta_i because WIMPs are also overproduced at large f_a. Hence, we conclude that setting the PQ scale via anthropics is highly unlikely. Instead, requiring soft SUSY breaking terms of order the gravity-mediation scale m_{3/2}~ 10-100 TeV places the mixed axion-neutralino dark matter abundance into the intermediate scale sweet zone where f_a~ 10^{11}-10^{12} GeV. We compare our analysis to the more general case of a generic SUSY DFSZ axion model with uniform selection on theta_i but leading to the measured dark matter abundance: this approach leads to a preference for f_a~ 10^{12} GeV.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

Particle physics models with Peccei-Quinn (PQ) symmetry breaking as a consequence of supersymmetry (SUSY) breaking are attractive in that they solve the strong CP problem with a SUSY DFSZ-like axion, link the SUSY breaking and PQ breaking intermediate mass scales and can resolve the SUSY mu problem with a naturalness-required weak scale mu term whilst soft SUSY breaking terms inhabit the multi-TeV regime as required by LHC sparticle mass limits and the Higgs mass measurement. On the negative ledger, models based on global symmetries suffer a generic gravity spoliation problem. We present a model based on the discrete R-symmetry Z_{24}^R-- which may emerge from compactification of 10-d Lorentzian spacetime in string theory-- where the mu term and dangerous proton decay and R-parity violating operators are either suppressed or forbidden while a gravity-safe PQ symmetry emerges as an accidental approximate global symmetry leading to a solution to the strong CP problem and a weak-scale/natural value for the mu term.

In spite of making a small contribution to total protoplanetary disk mass, dust affects the disk temperature by controlling absorption of starlight. As grains grow from their initial ISM-like size distribution, settling depletes the disk's upper layers of dust and decreases the optical depth, cooling the interior. Here we investigate the effect of collisional growth of dust grains and their dynamics on the thermal and optical profile of the disk, and explore the possibility that cooling induced by grain growth and settling could lead to gravitational instability. We develop a Monte Carlo dust collision model with a weighting technique and allow particles to collisionally evolve through sticking and fragmentation, along with vertical settling and turbulent mixing. We explore two disk models, the MMEN (minimum-mass extrasolar nebula), and a "heavy" disk with higher surface density than the MMEN, and perform simulations for both constant and spatially variable turbulence efficiency profile $\alpha(R,z)$. We then calculate mean wavelength-dependent opacities for the evolving disks and perform radiative transfer to calculate the temperature profile $T(R,z)$. Finally, we calculate the Toomre Q parameter, a measure of the disk's stability against self-gravity, for each disk model after it reaches a steady state dust-size distribution. We find that even weak turbulence can keep sub-micron sized particles stirred in the disk's upper layer, affecting its optical and thermal profiles, and the growth of large particles in the midplane can make a massive disk optically thick at millimeter wavelengths, making it difficult to calculate the surface density of dust available for planet formation in the inner disk. Also, for an initially massive disk, grain settling and growth can produce a drop in the Toomre Q parameter, driving the disk to $Q < 1.4$ and possibly triggering spiral instabilities.

The requirement of electroweak naturalness in supersymmetric (SUSY) models of particle physics necessitates light higgsinos not too far from the weak scale characterized by m(weak)~ m(W,Z,h)~100 GeV. On the other hand, LHC Higgs mass measurements and sparticle mass limits point to a SUSY breaking scale in the multi-TeV regime. Under such conditions, the lightest SUSY particle is expected to be a mainly higgsino-like neutralino with non-negligible gaugino components (required by naturalness). The computed thermal WIMP abundance in natural SUSY models is then found to be typically a factor 5-20 below its measured value. To gain concordance with observations, either an additional DM particle (the axion is a well-motivated possibility) must be present or additional non-thermal mechanisms must augment the neutralino abundance. We compare present direct and indirect WIMP detection limits to three natural SUSY models based on gravity-, anomaly- and mirage-mediation. We show that the case of natural higgsino-only dark matter where non-thermal production mechanisms augment its relic density, is essentially excluded by a combination of direct detection constraints from PandaX-II, LUX and Xenon-1t experiments, and by bounds from Fermi-LAT/MAGIC observations of gamma rays from dwarf spheroidal galaxies.