We examine the dynamical properties of interacting galaxies and the properties of shocked gas produced as a result of the interaction. We observed 22 galaxy mergers using the SparsePak IFU at Kitt Peak National Observatory (KPNO). The goal of the observations was to obtain the \ha\ velocity maps over the entire luminous parts of the galaxies including the faint tidal tails and to find extended shocks and outflows. Our sample consists of major and minor galaxy mergers with mass ratios $1<\mu<8$. We fit multiple kinematic components to the \ha\ and \nii\ emission lines, develop an MCMC code to robustly estimate the error of fit parameters, and use the F-test to determine the best number of kinematic components for each fiber. We use \nii/\ha\ and velocity dispersion of components to separate star-forming (HII) regions from shocks. We use the kinematics of the \ha\ emission from HII regions and an automated modeling method to put the first-ever constraints on the encounter parameters of one of the observed systems. Besides, we roughly estimate the fraction of shocked \ha\ emission , $\text{f}_\text{shocked}$, without taking extinction into account and examine the spatial distribution of shocks. We find that close galaxy pairs have, on average, a higher shock fraction than wide pairs, and coalesced mergers have the highest average $\text{f}_\text{shocked}$. In addition, galaxy pairs with more equal mass ratio tend to have a higher $\text{f}_\text{shocked}$. Combining the dynamical models from the literature and this work, we inspect trends between $\text{f}_\text{shocked}$ and dynamical encounter parameters. Our findings are generally consistent with shocks being produced either by the direct collision of the ISM or by the chain of events provoked by the tidal impulse during the first passage.

We explore the effect of using different kinematic tracers (HI and H$\alpha$) on reconstructing the encounter parameters of the Mice major galaxy merger (NGC 4676A/B). We observed the Mice using the SparsePak Integral Field Unit (IFU) on the WIYN telescope, and compared the H$\alpha$ velocity map with VLA HI observations. The relatively high spectral resolution of our data (R $\approx$ 5000) allows us to resolve more than one kinematic component in the emission lines of some fibers. We separate the H$\alpha$-[N II] emission of the star-forming regions from shocks using their [N II]/H$\alpha$ line ratio and velocity dispersion. We show that the velocity of star-forming regions agree with that of the cold gas (HI), particularly, in the tidal tails of the system. We reconstruct the morphology and kinematics of these tidal tails utilizing an automated modeling method based on the Identikit software package. We quantify the goodness of fit and the uncertainties of the derived encounter parameters. Most of the initial conditions reconstructed using H$\alpha$ and HI are consistent with each other, and qualitatively agree with the results of previous works. For example, we find 210$\pm^{50}_{40}$ Myrs, and 180$\pm^{50}_{40}$ Myrs for the time since pericenter, when modeling H$\alpha$ and HI kinematics, respectively. This confirms that in some cases, H$\alpha$ kinematics can be used instead of HI kinematics for reconstructing the initial conditions of galaxy mergers, and our automated modeling method is applicable to some merging systems.

Important but rare and subtle processes driving galaxy morphology and star-formation may be missed by traditional spiral, elliptical, irregular or S\'ersic bulge/disk classifications. To overcome this limitation, we use a principal component analysis of non-parametric morphological indicators (concentration, asymmetry, Gini coefficient, $M_{20}$, multi-mode, intensity and deviation) measured at rest-frame $B$-band (corresponding to HST/WFC3 F125W at 1.4 $< z <$ 2) to trace the natural distribution of massive ($>10^{10} M_{\odot}$) galaxy morphologies. Principal component analysis (PCA) quantifies the correlations between these morphological indicators and determines the relative importance of each. The first three principal components (PCs) capture $\sim$75 per cent of the variance inherent to our sample. We interpret the first principal component (PC) as bulge strength, the second PC as dominated by concentration and the third PC as dominated by asymmetry. Both PC1 and PC2 correlate with the visual appearance of a central bulge and predict galaxy quiescence. PC1 is a better predictor of quenching than stellar mass, as as good as other structural indicators (S\'ersic-n or compactness). We divide the PCA results into groups using an agglomerative hierarchical clustering method. Unlike S\'ersic, this classification scheme separates compact galaxies from larger, smooth proto-elliptical systems, and star-forming disk-dominated clumpy galaxies from star-forming bulge-dominated asymmetric galaxies. Distinguishing between these galaxy structural types in a quantitative manner is an important step towards understanding the connections between morphology, galaxy assembly and star-formation.

We develop and test an automated technique to model the dynamics of interacting galaxy pairs. We use Identikit (Barnes & Hibbard 2009, Barnes 2011) as a tool for modeling and matching the morphology and kinematics of the interacting pairs of equal-mass galaxies. In order to reduce the effect of subjective human judgement, we automate the selection of phase-space regions used to match simulations to data, and we explore how selection of these regions affects the random uncertainties of parameters in the best-fit model. In this work, we use an independent set of GADGET SPH simulations as input data to determine the systematic bias in the measured encounter parameters based on the known initial conditions of these simulations. We test both cold gas and young stellar components in the GADGET simulations to explore the effect of choosing HI vs. H$\alpha$ as the line of sight velocity tracer. We find that we can group the results into tests with good, fair, and poor convergence based on the distribution of parameters of models close to the best-fit model. For tests with good and fair convergence, we rule out large fractions of parameter space and recover merger stage, eccentricity, pericentric distance, viewing angle, and initial disc orientations within 3$\sigma$ of the correct value. All of tests on prograde-prograde systems have either good or fair convergence. The results of tests on edge-on discs are less biased than face-on tests. Retrograde and polar systems do not converge and may require constraints from regions other than the tidal tails and bridges.