Source author record

Sungsoo S. Kim

Sungsoo S. Kim appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2020arXiv

Effect of Bars on Evolution of SDSS Spiral Galaxies

We explore the significance of bars in triggering central star formation (SF) and AGN activity for spiral galaxy evolution using a volume-limited sample with $0.020<z<0.055$, $M_{\rm r}<-19.5$, and $σ>70\rm km s^{-1}$ selected from SDSS DR7. On a central SF rate-$σ$ plane, we measure the fraction of galaxies with strong bars in our sample and also the AGN fractions for barred and non-barred galaxies, respectively. The comparison between the bar and AGN fractions reveals a causal connection between the two phenomena of SF quenching and AGN activity. A massive BH and abundant gas fuels are sufficient conditions to trigger AGNs. We infer that the AGNs triggered by satisfying the two conditions drive the strong AGN feedback, suddenly suppressing the central SF and leaving the SF sequence. We find that in galaxies where either of the two conditions is not sufficient, bars are a great help for the AGN triggering, accelerating the entire process of evolution, which is particularly evident in pseudo-bulge galaxies. All of our findings are obtained only when plotted in terms of their central velocity dispersion and central SFR (not galactic scale SFR), indicating that the AGN-driven SF quenching is confined in the central kpc region.

preprint2020arXiv

Making top-heavy IMFs from canonical IMFs near the Galactic Centre

We show that dynamical evolution in a strong (Galactic Centre-like) tidal field can create clusters that would appear to have very top-heavy IMFs. The tidal disruption of single star forming events can leave several bound `clusters' spread along 20 pc of the orbit within 1-2 Myr. These surviving (sub)clusters tend to contain an over-abundance of massive stars, with low-mass stars tending to be spread along the whole `tidal arm'. Therefore observing a cluster in a strong tidal field with a top-heavy IMF might well not mean the stars formed with a top-heavy IMF.

preprint2019arXiv

Direct effects of the environment on AGN triggering in SDSS spiral galaxies: merger-AGN connection

We examine whether galaxy environments directly affect triggering nuclear activity in Sloan Digital Sky Survey (SDSS) local spiral galaxies using a volume-limited sample with the $r$-band absolute magnitude $M_{r} < -19.0$ and $0.02 < z < 0.055$ selected from the SDSS Data Release 7. To avoid incompleteness of the central velocity dispersion $σ$ of the volume-limited sample and to fix the black hole mass affecting AGN activity, we limit the sample to a narrow $σ$ range of $130$ km s$^{-1}<σ<200$ km s$^{-1}$. We define a variety of environments as a combination of neighbour interactions and local density on a galaxy. After the central star formation rate (which is closely related to AGN activity level) is additionally restricted, the direct impact of the environment is unveiled. In the outskirts of rich clusters, red spiral galaxies show a significant excess of the AGN fraction despite the lack of central gas. We argue that they have been pre-processed before entering the rich clusters, and due to mergers or strong encounters in the in-fall region, their remaining gases efficiently lose angular momentum. We investigate an environment in which many star-forming galaxies coexist with a few starburst-AGN composite hosts having the highest [OIII] luminosity. We claim that they are a gas-rich merger product in groups or are group galaxies in-falling into clusters, indicating that many AGN signatures may be obscured following the merger events.

preprint2019arXiv

The Gas Accretion History of Low Mass Halos within the Cosmic Web from Cosmological Simulations

Using high resolution hydrodynamical cosmological simulations, we study the gas accretion history of low mass halos located in a field-like, low density environment. We track their evolution individually from the early, pre-reionization era, through reionization, and beyond until $z=0$. Before reionization, low mass halos accrete cool cosmic web gas at a very rapid rate, often reaching the highest gas mass they will ever have. But when reionization occurs, we see that almost all halos lose significant quantities of their gas content, although some respond less quickly than others. We find that the response rate is influenced by halo mass first, and secondarily by their internal gas density at the epoch of reionization. Reionization also fully ionises the cosmic web gas by z$\sim$6. As a result, the lowest mass halos (M$\sim$10$^6~h^{-1}$M$_\odot$ at $z=6$) can never again re-accrete gas from the cosmic web, and by $z\sim5$ have lost all their internal gas to ionisation, resulting in a halt in star formation at this epoch. However, more massive halos can recover from their gas mass loss, and re-accrete ionised cosmic web gas. We find the efficiency of this re-accretion is a function of halo mass first, followed by local surrounding gas density. Halos that are closer to the cosmic web structure can accrete denser gas more rapidly. We find that our lower mass halos have a sweet spot for rapid, dense gas accretion at distances of roughly 1-5 virial radii from the most massive halos in our sample ($>$10$^8~h^{-1}$M$_\odot$), as these tend to be embedded deeply within the cosmic web.

preprint2016arXiv

Low-end mass function of the Quintuplet cluster

The Quintuplet and Arches clusters, which were formed in the harsh environment of the Galactic Center (GC) a few million years ago, have been excellent targets for studying the effects of a star-forming environment on the initial mass function (IMF). In order to estimate the shape of the low-end IMF of the Arches cluster, Shin & Kim devised a novel photometric method that utilizes pixel intensity histograms (PIHs) of the observed images. Here, we apply the PIH method to the Quintuplet cluster and estimate the shape of its low-end IMF below the magnitude of completeness limit as set by conventional photometry. We found that the low-end IMF of the Quintuplet is consistent with that found for the Arches cluster--Kroupa MF, with a significant number of low-mass stars below 1 $M_{\odot}$. We conclude that the most likely IMFs of the Quintuplet and the Arches clusters are not too different from the IMFs found in the Galactic disc. We also find that the observed PIHs and stellar number density profiles of both clusters are best reproduced when the clusters are assumed to be at three-dimensional distances of approximately 100 pc from the GC.

preprint2015arXiv

Stochastic Model of the Spin Distribution of Dark Matter Halos

We employ a stochastic approach to probing the origin of the log-normal distributions of halo spin in N-body simulations. After analyzing spin evolution in halo merging trees, it was found that a spin change can be characterized by a stochastic random walk of angular momentum. Also, spin distributions generated by random walks are fairly consistent with those directly obtained from N-body simulations. We derived a stochastic differential equation from a widely used spin definition and measured the probability distributions of the derived angular momentum change from a massive set of halo merging trees. The roles of major merging and accretion are also statistically analyzed in evolving spin distributions. Several factors (local environment, halo mass, merging mass ratio, and redshift) are found to influence the angular momentum change. The spin distributions generated in the mean-field or void regions tend to shift slightly to a higher spin value compared with simulated spin distributions, which seems to be caused by the correlated random walks. We verified the assumption of randomness in the angular momentum change observed in the N-body simulation and detected several degrees of correlation between walks, which may provide a clue for the discrepancies between the simulated and generated spin distributions in the voids. However, the generated spin distributions in the group and cluster regions successfully match the simulated spin distribution. We also demonstrated that the log-normality of the spin distribution is a natural consequence of the stochastic differential equation of the halo spin, which is well described by the Geometric Brownian Motion model.

preprint2014arXiv

How does a low-mass cut-off in the stellar IMF affect the evolution of young star clusters?

We investigate how different stellar initial mass functions (IMFs) can affect the mass loss and survival of star clusters. We find that IMFs with radically different low-mass cut-offs (between 0.1 and 2 Msun) do not change cluster destruction time-scales as much as might be expected. Unsurprisingly, we find that clusters with more high-mass stars lose relatively more mass through stellar evolution, but the response to this mass loss is to expand and hence significantly slow their dynamical evolution. We also argue that it is very difficult, if not impossible, to have clusters with different IMFs that are initially "the same", since the mass, radius and relaxation times depend on each other and on the IMF in a complex way. We conclude that changing the IMF to be biased towards more massive stars does speed up mass loss and dissolution, but that it is not as dramatic as might be thought.

preprint2014arXiv

Low-End Mass Function of the Arches Cluster

The initial mass function (IMF) of the Arches cluster, which was formed a few million years ago in the harsh environment of the Galactic center (GC), has long been a target of interest to those who study the GC and the theory of star formation. The distinct star-forming conditions in the GC might have caused the cluster to have a shallower slope or an elevated lower mass cutoff in its IMF. But its mass function has been revealed only down to 1-2 Msun (the lower limit of resolved stars), and the low- end mass function of the Arches is still unknown. To estimate the unresolved part of the Arches mass function, we have devised a novel photometric method that involves the histogram of pixel intensities in the observed image, which contains information on the unresolved, faint stars. By comparing the pixel intensity histograms (PIHs) of numerous artificial images constructed from model IMFs with the observed PIH, we find that the best-fit model IMF for the Arches cluster has a cutoff mass less than or similar to 0.1 Msun and a shape very close to that of the Kroupa MF. Our findings imply that the IMF of the Arches cluster is similar to those found in the Galactic disk.

preprint2014arXiv

Systematic Effects on the Genus Topology of Large Scale Structure of the Universe

Large-scale structure of the universe is a useful cosmological probe of the primordial non-Gaussianity and the expansion history of the universe because its topology does not change with time in the linear regime in the standard paradigm of structure formation. However, when the topology of iso-density contour surfaces is measured from an observational data, many systematic effects are introduced due to the finite size of pixels used to define the density field, non-linear gravitational evolution, redshift-space distortion, shot noise (discrete sampling), and bias in the distribution of the density field tracers. We study the various systematic effects on the genus curve to a great accuracy by using the Horizon Run 2 simulation of a ΛCDM cosmology. We numerically measure the genus curve from the gravitationally evolved matter and dark matter halo density fields. It is found that all the non-Gaussian deviations due to the systematic effects can be modeled by using a few low-order Hermite polynomials from H0 to H4. We compare our results with the analytic theories whenever possible, and find many new terms in the Hermite series that are making significant contributions to the non-Gaussian deviations. In particular, it is found that the amplitude drop of the genus curve due to the non-linear gravitational evolution can be accurately modeled by two terms H0 and H2 with coefficients both proportional to σ_0^2, the mean-square density fluctuation.

preprint2014arXiv

The Processing of the Clumpy Molecular Gas in the Galactic Center and the Star-Formation

We present the Green Bank 100m Telescope (GBT) mapping observations of CS 1-0, and the Submillimeter Array (SMA) 157-pointings mosaic of the 0.86 mm dust continuum emission as well as several warm and dense gas tracers, in the central ~20 pc area in Galactic Center. The unprecedentedly large field-of-view and the high angular resolution of our SMA dust image allow the identification of abundant 0.1-0.2 pc scale dense gas clumps. We found that in the Galactic Center, the Class I methanol masers are excellently correlated with the dense gas clumps. However, on the ~0.1 pc scale, these dense gas clumps still have a extremely large linewidth (FWHM~10-20 km/s). Simple calculations suggest that the identified clumps can be possibly the pressurized gas reservoir feeding the formation of 1-10 solar-mass stars. These gas clumps may be the most promising candidates for ALMA to resolve the high-mass star-formation in the Galactic center.

preprint2013arXiv

ISM Processing in the Inner 20 pc in Galactic center

We present the Submillimeter Array (SMA) 157-pointing mosaic in 0.86 mm dust continuum emission with 5$".1\times4"$.2 angular resolution, and the NRAO Green Bank 100m Telescope (GBT) observations of the CS/C$^{34}$S/$^{13}$CS 1-0 and SiO 1-0 emission with $\le20"$$\times18"$ angular resolution. The dust continuum image marginally resolves at least several tens of 10-10$^{2}$ $M_{\odot}$ dense clumps in the 5$'$ field including the circumnuclear disk (CND) and the exterior gas streamers. There is very good agreement between the high resolution dust continuum map of the CND and all previous molecular line observations. As the dust emission is the most reliable optically thin tracer of the mass, free from most chemical and excitation effects, we demonstrate the reality of the abundant localized structures within the CND, and their connection to external gas structures. From the spectral line data, the velocity dispersions of the dense clumps and their parent molecular clouds are $\sim$10-20 times higher than their virial velocity dispersions. This supports the idea that the CND and its immediate environment may not be stationary or stable structures. Some of the dense gas clumps are associated with 22 GHz water masers and 36.2 GHz and 44.1 GHz CH$_{3}$OH masers. However, we do not find clumps which are bound by the gravity of the enclosed molecular gas. Hence, the CH$_{3}$OH or H$_{2}$O maser emission may be due to strong (proto)stellar feedback, which may be dispersing some of the gas clumps.

preprint2013arXiv

Topology of Luminous Red Galaxies from the Sloan Digital Sky Survey

We present measurements of the genus topology of luminous red galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS) Data Release 7 catalog, with unprecedented statistical significance. To estimate the uncertainties in the measured genus, we construct 81 mock SDSS LRG surveys along the past light cone from the Horizon Run 3, one of the largest N-body simulations to date that evolved 7210^3 particles in a 10815 Mpc/h size box. After carefully modeling and removing all known systematic effects due to finite pixel size, survey boundary, radial and angular selection functions, shot noise and galaxy biasing, we find the observed genus amplitude to reach 272 at 22 Mpc/h smoothing scale with an uncertainty of 4.2%; the estimated error fully incorporates cosmic variance. This is the most accurate constraint of the genus amplitude to date, which significantly improves on our previous results. In particular, the shape of the genus curve agrees very well with the mean topology of the SDSS LRG mock surveys in the LCDM universe. However, comparison with simulations also shows small deviations of the observed genus curve from the theoretical expectation for Gaussian initial conditions. While these discrepancies are mainly driven by known systematic effects such as those of shot noise and redshift-space distortions, they do contain important cosmological information on the physical effects connected with galaxy formation, gravitational evolution and primordial non-Gaussianity. We address here the key role played by systematics on the genus curve, and show how to accurately correct for their effects to recover the topology of the underlying matter. In a forthcoming paper, we provide an interpretation of those deviations in the context of the local model of non-Gaussianity.

preprint2012arXiv

A Second-order bias model for the Logarithmic Halo Mass Density

We present an analytic model for the local bias of dark matter halos in a LCDM universe. The model uses the halo mass density instead of the halo number density and is searched for various halo mass cuts, smoothing lengths, and redshift epoches. We find that, when the logarithmic density is used, the second-order polynomial can fit the numerical relation between the halo mass distribution and the underlying matter distribution extremely well. In this model the logarithm of the dark matter density is expanded in terms of log halo mass density to the second order. The model remains excellent for all halo mass cuts (from M_{cut}=3\times10^{11}$ to $3\times10^{12}h^{-1}M_{\odot}$), smoothing scales (from $R=5h^{-1}$Mpc to $50h^{-1}$Mpc), and redshift ranges (from z=0 to 1.0) considered in this study. The stochastic term in the relation is found not entirely random, but a part of the term can be determined by the magnitude of the shear tensor.

preprint2012arXiv

The Challenge of the Largest Structures in the Universe to Cosmology

Large galaxy redshift surveys have long been used to constrain cosmological models and structure formation scenarios. In particular, the largest structures discovered observationally are thought to carry critical information on the amplitude of large-scale density fluctuations or homogeneity of the universe, and have often challenged the standard cosmological framework. The Sloan Great Wall (SGW) recently found in the Sloan Digital Sky Survey (SDSS) region casts doubt on the concordance cosmological model with a cosmological constant (i.e. the flat LCDM model). Here we show that the existence of the SGW is perfectly consistent with the LCDM model, a result that only our very large cosmological N-body simulation (the Horizon Run 2, HR2) could supply. In addition, we report on the discovery of a void complex in the SDSS much larger than the SGW, and show that such size of the largest void is also predicted in the LCDM paradigm. Our results demonstrate that an initially homogeneous isotropic universe with primordial Gaussian random phase density fluctuations growing in accordance with the General Relativity, can explain the richness and size of the observed large-scale structures in the SDSS. Using the HR2 simulation we predict that a future galaxy redshift survey about four times deeper or with 3 magnitude fainter limit than the SDSS should reveal a largest structure of bright galaxies about twice as big as the SGW.

preprint2011arXiv

Nuclear Star-Forming Ring of the Milky Way: Simulations

We present hydrodynamic simulations of gas clouds in the central kpc region of the Milky Way that is modeled with a three-dimensional bar potential. Our simulations consider realistic gas cooling and heating, star formation, and supernova feedback. A ring of dense gas clouds forms as a result of X1-X2 orbit transfer, and our potential model results in a ring radius of ~200 pc, which coincides with the extraordinary reservoir of dense molecular clouds in the inner bulge, the Central Molecular Zone (CMZ). The gas clouds accumulated in the CMZ can reach high enough densities to form stars, and with an appropriate choice of simulation parameters, we successfully reproduce the observed gas mass and the star formation rate (SFR) in the CMZ, ~2x10^7 Msun and ~0.1 Msun/yr. Star formation in our simulations takes place mostly in the outermost X2 orbits, and the SFR per unit surface area outside the CMZ is much lower. These facts suggest that the inner Galactic bulge may harbor a mild version of the nuclear star-forming rings seen in some external disk galaxies. Furthermore, from the relatively small size of the Milky Way's nuclear bulge, which is thought to be a result of sustained star formation in the CMZ, we infer that the Galactic inner bulge probably had a shallower density profile or stronger bar elongation in the past.

preprint2010arXiv

Galaxy Clustering Topology in the Sloan Digital Sky Survey Main Galaxy Sample: a Test for Galaxy Formation Models

We measure the topology of the main galaxy distribution using the Seventh Data Release of the Sloan Digital Sky Survey, examining the dependence of galaxy clustering topology on galaxy properties. The observational results are used to test galaxy formation models. A volume-limited sample defined by $M_r<-20.19$ enables us to measure the genus curve with amplitude of $G=378$ at $6h^{-1}$Mpc smoothing scale, with 4.8\% uncertainty including all systematics and cosmic variance. The clustering topology over the smoothing length interval from 6 to $10 h^{-1}$Mpc reveals a mild scale-dependence for the shift ($Δν$) and void abundance ($A_V$) parameters of the genus curve. We find substantial bias in the topology of galaxy clustering with respect to the predicted topology of the matter distribution, which varies with luminosity, morphology, color, and the smoothing scale of the density field. The distribution of relatively brighter galaxies shows a greater prevalence of isolated clusters and more percolated voids. Even though early (late)-type galaxies show topology similar to that of red (blue) galaxies, the morphology dependence of topology is not identical to the color dependence. In particular, the void abundance parameter $A_V$ depends on morphology more strongly than on color. We test five galaxy assignment schemes applied to cosmological N-body simulations of a $Λ$CDM universe to generate mock galaxies: the Halo-Galaxy one-to-one Correspondence model, the Halo Occupation Distribution model, and three implementations of Semi-Analytic Models (SAMs). None of the models reproduces all aspects of the observed clustering topology; the deviations vary from one model to another but include statistically significant discrepancies in the abundance of isolated voids or isolated clusters and the amplitude and overall shift of the genus curve. (Abridged)

preprint2009arXiv

Galactic Warps in Triaxial Halos

We study the behaviors of galactic disks in triaxial halos both numerically and analytically to see if warps can be excited and sustained in triaxial potentials. We consider the following two scenarios: 1) galactic disks that are initially tilted relative to the equatorial plane of the halo (for a pedagogical purpose), and 2) tilted infall of dark matter relative to the equatorial plane of the disk and the halo. With numerical simulations of 100,000 disk particles in a fixed halo potential, we find that in triaxial halos, warps can be excited and sustained just as in spherical or axisymmetric halos but they show some oscillatory behaviors and even can be transformed to a polar-ring system if the halo has a prolate-like triaxiality. The non-axisymmetric component of the halo causes the disk to nutate, and the differential nutation between the inner and outer parts of the disk generally makes the magnitude of the warp slightly diminish and fluctuate. We also find that warps are relatively weaker in oblate and oblate-like triaxial halos, and since these halos are the halo configurations of disk galaxies inferred by cosmological simulations, our results are consistent with the fact that most of the observed warps are quite weak. We derive approximate formulae for the torques exerted on the disk by the triaxial halo and the dark matter torus, and with these formulae we successfully describe the behaviors of the disks in our simulations. The techniques used in deriving these formulae could be applied for realistic halos with more complex structures.

preprint2008arXiv

Incidence of high-amplitude delta Scuti type variable stars

An order-of-magnitude estimate for the incidence of high-amplitude delta Scuti-type variable stars (HADS) in the delta Scuti area of the H-R diagram is calculated. Using a model for the stellar distribution in the Milky Way, we calculate the number of stars that are expected to fall in the delta Scuti area of the H-R diagram within the magnitude range and sky coverage of the ROTSE Survey for Variables I (RSV1). The incidence of the HADS phenomenon is then obtained by comparing the number of stars calculated by the model and the actual, observed number of HADS in the RSV1. We find that ~0.3 % of the stars that lie in the delta Scuti area of the H-R diagram within the RSV1 observational limits exhibit the HADS phenomenon. This number is much lower than the incidence of the low-amplitude delta Scuti stars (LADS), ~>1/3, implying that the HADS phenomenon takes place in a very small fraction of stars and/or its duration is very short, compared to the LADS.

preprint2006arXiv

Theoretical Isochrones with Extinction in the K Band. II. J - K versus K

We calculate theoretical isochrones in a consistent way for five filter pairs near the J and K band atmospheric windows (J-K, J-K', J-Ks, F110W-F205W, and F110W-F222M) using the Padova stellar evolutionary models of Girardi et al. We present magnitude transformations between various K-band filters as a function of color. Isochrones with extinction of up to 6 mag in the K band are also presented. As found for the filter pairs composed of H & K band filters, we find that the reddened isochrones of different filter pairs behave as if they follow different extinction laws, and that the extinction curves of Hubble Space Telescope NICMOS filter pairs in the color-magnitude diagram are considerably nonlinear. Because of these problems, extinction values estimated with NICMOS filters can be in error by up to 1.3 mag. Our calculation suggests that the extinction law implied by the observations of Rieke et al for wavelengths between the J and K bands is better described by a power-law function with an exponent of 1.66 instead of 1.59, which is commonly used with an assumption that the transmission functions of J and K filters are Dirac delta functions.

preprint2003arXiv

The NIRSPEC Brown Dwarf Spectroscopic Survey. I. Low-Resolution Near-Infrared Spectra

We present the first results of a near-infrared (0.96-2.31 micron) spectroscopic survey of M, L, and T dwarfs obtained with NIRSPEC on the Keck II telescope. Our new survey has a resolving power of R = 2000 and is comprised of two major data sets: 53 J-band (1.14-1.36 micron) spectra covering all spectral types from M6 to T8 with at least two members in each spectral subclass (wherever possible), and 25 flux-calibrated spectra from 1.14 to 2.31 microns for most spectral classes between M6 and T8. Sixteen of these 25 objects have additional spectral coverage from 0.96-1.14 microns to provide overlap with optical spectra. Spectral flux ratio indexes for prominent molecular bands are derived and equivalent widths (EWs) for several atomic lines are measured. We find that a combination of four H2O and two CH4 band strengths can be used for spectral classification. Weak (EW~1-2 angstrom) atomic lines of Al I and Ca I disappear at the boundary between M and L types.

preprint2002arXiv

Massive Stars in the Arches Cluster

We present and use new spectra and narrow-band images, along with previously published broad-band images, of stars in the Arches cluster to extract photometry, astrometry, equivalent width, and velocity information. The data are interpreted with a wind/atmosphere code to determine stellar temperatures, luminosities, mass-loss rates, and abundances. We have doubled the number of known emission-line stars, and we have also made the first spectroscopic identification of the main sequence for any population in the Galactic Center. We conclude that the most massive stars are bona-fide Wolf-Rayet (WR) stars and are some of the most massive stars known, having M_{initial} > 100 Msun, and prodigious winds, Mdot > 10^{-5} Msun yr^{-1}, that are enriched with helium and nitrogen; with these identifications, the Arches cluster contains about 5% of all known WR stars in the Galaxy. We find an upper limit to the velocity dispersion of 22 kms^{-1}, implying an upper limit to the cluster mass of 7(10^4) Msun within a radius of 0.23 pc; we also estimate the bulk heliocentric velocity of the cluster to be v_{cluster,odot} approximately +95 kms^{-1}.

preprint1997arXiv

Two-Component Fokker-Planck Models for the Evolution of Isolated Globular Clusters

Two-component (normal and degenerate stars) models are the simplest realization of clusters with a mass spectrum because high mass stars evolve quickly into degenerates, while low mass stars remain on the main-sequence for the age of the universe. Here we examine the evolution of isolated globular clusters using two-component Fokker-Planck (FP) models that include heating by binaries formed in tidal capture and in three-body encounters. Three-body binary heating dominates and the postcollapse expansion is self-similar, at least in models with total mass M <= 3 x 10^5 M_\odot, initial half-mass radius r_{h,i} >= 5 pc, component mass ratio m_2/m_1 <= 2, and number ratio N_1/N_2 <= 300 when m_2=1.4 M_\odot. We derive scaling laws for ρ_c, v_c, r_c, and r_h as functions of m_1/m_2, N, M, and time t from simple energy-balance arguments, and these agree well with the FP simulations. We have studied the conditions under which gravothermal oscillations (GTOs) occur. If E_{tot} and E_c are the energies of the cluster and of the core, respectively, and t_{rh} and t_c are their relaxation times, then ε\equiv (E_{tot}/t_{rh})/(E_c/t_{rc}) is a good predictor of GTOs: all models with ε>0.01 are stable, and all but one with ε< 0.01 oscillate. We derive a scaling law for εagainst N and m_1/m_2 and compared with our numerical results. Clusters with larger m_2/m_1 or smaller N are stabler.