We used the Green Bank Telescope to detect molecular lines observed toward Mid-Infrared (MIR) bubbles N62, N65, N90, and N117. The bubbles were selected from Watson et al. (2016) who detected non-Gaussian CS (1-0) emission lines toward the bubbles. Two of the bubbles are adjacent to infrared dark clouds (IRDCs); we examined these sources for evidence of interaction between the bubble rim and IRDC. The other two bubbles contain YSOs interior to the bubble rim; in these sources we observed the gas near the YSOs. We detect CS (1-0) emission toward all of the sources, and in several pointings the CS emission shows non-Gaussian line shapes. HC$_3$N (5-4), C$^{34}$S (1-0), CH$_3$OH (1-0), and SiO (v=0) (1-0) were also detected in some pointings. We calculate column densities and abundances for the detected molecules. We compare the velocity of optically-thick CS emission with the velocity of the other, optically thin lines to look for evidence of infall. We find that even in pointings with non-Gaussian CS emission, our detections do not support an infall model. We interpret the kinematics of the gas in N62, N65, and N117 as likely evidence of multiple clouds along the line of sight moving at slightly offset velocities. We do not detect evidence of bubble rims interacting with IRDCs in N62 or N90. The gas interior to bubbles appears more disrupted than the gas in the IRDCs. N65 shows significantly stronger emission lines than the other sources, as well as the most complicated non-Gaussian line shapes.

We re-examine the Lem\^aitre-Tolman-Bondi (LTB) solutions with a dust source admitting symmetry centers. We consider as free parameters of the solutions the initial value functions: $Y_i$, $\rho_i$ and $\Ri$, obtained by restricting the curvature radius, $Y\equiv \sqrt{g_{\theta\theta}}$, the rest mass density, $\rho$, and the 3-dimensional Ricci scalar of the rest frames, $\R$, to an arbitrary regular Cauchy hypersurface, $\Ti$, marked by constant cosmic time ($t=t_i$). Using $Y_i$ to fix the radial coordinate and the topology (homeomorphic class) of $\Ti$, and scaling the time evolution in terms of an adimensional scale factor $y=Y/Y_i$, we show that the dynamics, regularity conditions and geometric features of the models are determined by $\rho_i$, $\Ri$ and by suitably constructed volume averages and contrast functions expressible in terms of invariant scalars defined in $\Ti$. These quantities lead to a straightforward characterization of initial conditions in terms of the nature of the inhomogeneity of $\Ti$, as density and/or curvature overdensities (``lumps'') and underdensities (''voids'') around a symmetry center. In general, only models with initial density and curvature lumps evolve without shell crossing singularities, though special classes of initial conditions, associated with a simmultaneous big bang, allow for a regular evolution for initial density and curvature voids. Specific restrictions are found so that a regular evolution for $t>t_i$ is possible for initial voids. A step-by-step guideline is provided for using the new variables in the construction of LTB models and for plotting all relevant quantities.