In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

Self-organization properties of sustained magnetized plasma are applied to selected solar data to understand solar magnetic fields. Torsional oscillations are speed-up and slow-down bands of the azimuthal flow that correlate with the solar cycle, and they imply the existence of a symmetric solar dynamo with a measured polar flux of 3x10^14 Wb. It is shown that the solar dynamo is thin (~0.1 Mm gradient scale size) and powerful (~10^23 W). These properties are found from the amplitude of the torsional oscillations and the relationship of their velocity contours to solar magnetograms supports the result. The dynamo has enough power to heat the chromosphere and to power the corona and the solar wind. The dynamo also causes a rigid rotation of the heliosphere out to at least the corona and the relationship of the rotation of the corona to solar magnetograms supports this result as well. The thin solar dynamo sustains a thin stable minimum energy state that seems to be covering most of the solar surface just below the photosphere. The magnetic field lines of the minimum energy state should be parallel to the solar surface and rotate with distance from the surface with 2{\pi} radians of rotation in ~1 Mm Resistive diffusion helps to push the magnetic fields to the surface and the global magnetic structure (GMS) seems to lose {\pi} radians every 11 years, causing the observed 180 degree flipping of the solar magnetic field. The thin sheets of magnetized plasma in solar prominences may be the lost thin sheets of the GMS. For completeness, the formation of sunspots, CMEs and flares is discussed.

The solar dynamo and the solar Global internal Magnetic Structure (GMS) appear to be a thin ($\sim$2 Mm thick) structure near ($\sim$1 Mm below) the solar surface. Evidence for these properties are found from the amplitude of the torsional oscillations and in their velocity contours relationship to solar magnetogram; the power to the chromosphere; power to the corona and the solar wind; the current in the helio-current-sheet measured at the radius of the orbit of Earth; the calculated size ($\sim$1 Mm) of the expanding polar flux when it enters the photosphere; and from the observation that solar magnetic activity is generated near the surface. A thin stable minimum energy state seems to be covering most of the solar surface just below the photosphere. The magnetic field lines should be parallel to the solar surface and rotate with distance from the surface for 2$\pi$ radians in $\sim$2 Mm. Resistive diffusion helps to push the magnetic fields to the surface and the GMS seems to lose $\pi$ radians every 11 years, causing the observed 180$^\circ$ flipping of the solar magnetic fields including the flipping of the polar flux. Further evidence for this GMS and its loss is that solar prominences are made of thin sheets of magnetized plasma, which are, likely, remnants of the lost thin sheet of the GMS. The loss process is consistent with the butterfly pattern of the sunspots and with the differences observed between solar maximum and solar minimum in the corona. The solar dynamo drives current parallel to the polar flux, which, in turn, sustains the GMS using cross-field current drive. For completeness, the formation of sunspots, CMEs and flares is discussed.

We study dynamical mass measurements of galaxy clusters contaminated by interlopers and show that a modern machine learning (ML) algorithm can predict masses by better than a factor of two compared to a standard scaling relation approach. We create two mock catalogs from Multidark's publicly available $N$-body MDPL1 simulation, one with perfect galaxy cluster membership information and the other where a simple cylindrical cut around the cluster center allows interlopers to contaminate the clusters. In the standard approach, we use a power-law scaling relation to infer cluster mass from galaxy line-of-sight (LOS) velocity dispersion. Assuming perfect membership knowledge, this unrealistic case produces a wide fractional mass error distribution, with a width of $\Delta\epsilon\approx0.87$. Interlopers introduce additional scatter, significantly widening the error distribution further ($\Delta\epsilon\approx2.13$). We employ the support distribution machine (SDM) class of algorithms to learn from distributions of data to predict single values. Applied to distributions of galaxy observables such as LOS velocity and projected distance from the cluster center, SDM yields better than a factor-of-two improvement ($\Delta\epsilon\approx0.67$) for the contaminated case. Remarkably, SDM applied to contaminated clusters is better able to recover masses than even the scaling relation approach applied to uncontaminated clusters. We show that the SDM method more accurately reproduces the cluster mass function, making it a valuable tool for employing cluster observations to evaluate cosmological models.

We present a modern machine learning approach for cluster dynamical mass measurements that is a factor of two improvement over using a conventional scaling relation. Different methods are tested against a mock cluster catalog constructed using halos with mass >= 10^14 Msolar/h from Multidark's publicly-available N-body MDPL halo catalog. In the conventional method, we use a standard M(sigma_v) power law scaling relation to infer cluster mass, M, from line-of-sight (LOS) galaxy velocity dispersion, sigma_v. The resulting fractional mass error distribution is broad, with width=0.87 (68% scatter), and has extended high-error tails. The standard scaling relation can be simply enhanced by including higher-order moments of the LOS velocity distribution. Applying the kurtosis as a correction term to log(sigma_v) reduces the width of the error distribution to 0.74 (16% improvement). Machine learning can be used to take full advantage of all the information in the velocity distribution. We employ the Support Distribution Machines (SDMs) algorithm that learns from distributions of data to predict single values. SDMs trained and tested on the distribution of LOS velocities yield width=0.46 (47% improvement). Furthermore, the problematic tails of the mass error distribution are effectively eliminated. Decreasing cluster mass errors will improve measurements of the growth of structure and lead to tighter constraints on cosmological parameters.