Context. An automatic tool to derive structural parameters of semi-resolved star clusters located in crowded stellar fields in nearby galaxies is needed for homogeneous processing of archival frames. Aims. We have developed a program that automatically derives the structural parameters of star clusters and estimates errors by accounting for individual stars and variable sky background. Methods. Models of observed frames consist of the cluster's surface brightness distribution, convolved with a point spread function; the stars, represented by the same point spread function; and a smoothly variable sky background. The cluster's model is fitted within a large radius by using the Levenberg-Marquardt and Markov chain Monte Carlo algorithms to derive structural parameters, the flux of the cluster, and individual fluxes of all well-resolved stars. Results. FitClust, a program to derive structural parameters of semi-resolved clusters in crowded stellar fields, was developed and is available for free use. The program was tested on simulated cluster frames, and was used to measure clusters of the M31 galaxy in Subaru Suprime-Cam frames. Conclusions. Accounting for bright resolved stars and variable sky background significantly improves the accuracy of derived structural parameters of star clusters. However, their uncertainty remains dominated by the stochastic noise of unresolved stars.

We improved the stochastic 2-D galaxy disk models (Mineikis & Vansevi\v{c}ius 2014a) by introducing enriched gas outflows from galaxies and synthetic color-magnitude diagrams of stellar populations. To test the models, we use the HST/ACS stellar photometry data in four fields located along the major axis of the galaxy M33 (Williams et al. 2009) and demonstrate the potential of the models to derive 2-D star formation histories in the resolved disk galaxies.

Context. When trying to derive the star cluster physical parameters of the M33 galaxy using broad-band unresolved ground-based photometry, previous studies mainly made use of simple stellar population models, shown in the recent years to be oversimplified. Aims. In this study, we aim to derive the star cluster physical parameters (age, mass, and extinction; metallicity is assumed to be LMC-like for clusters with age below 1\,Gyr and left free for older clusters) of this galaxy using models that take stochastic dispersion of cluster integrated colors into account. Methods. We use three recently published M33 catalogs of cluster optical broad-band photometry in standard $UBVRI$ and in CFHT/MegaCam $u^{*}g'r'i'z'$ photometric systems. We also use near-infrared $JHK$ photometry that we derive from deep 2MASS images. We derive the cluster parameters using a method that takes into account the stochasticity problem, presented in previous papers of this series. Results. The derived differential age distribution of the M33 cluster population is composed of a two-slope profile indicating that the number of clusters decreases when age gets older. The first slope is interpreted as the evolutionary fading phase of the cluster magnitudes, and the second slope as the cluster disruption. The threshold between these two phases occurs at $\sim$300\,Myrs, comparable to what is observed in the M31 galaxy. We also model by use of artificial clusters the ability of the cluster physical parameter derivation method to correctly derive the two-slope profile for different photometric systems tested.

This study is the third of a series that investigates the degeneracy and stochasticity problems present in the determination of physical parameters such as age, mass, extinction, and metallicity of partially resolved or unresolved star cluster populations situated in external galaxies when using broad-band photometry. This work tests the derivation of parameters of artificial star clusters using models with fixed and free metallicity for the WFC3+ACS photometric system. Then the method is applied to derive parameters of a sample of 203 star clusters in the Andromeda galaxy observed with the HST. Following Papers I \& II, the star cluster parameters are derived using a large grid of stochastic models that are compared to the observed cluster broad-band integrated WFC3+ACS magnitudes. We derive the age, mass, and extinction of the sample of M31 star clusters with one fixed metallicity in agreement with previous studies. Using artificial tests we demonstrate the ability of the WFC3+ACS photometric system to derive the metallicity of star clusters. We show that the metallicity derived using photometry of 36 massive M31 star clusters is in a good agreement with the metallicity previously derived using spectroscopy taken from literature.

We investigate the effects of stochasticity on the observed galaxy parameters by comparing our stochastic star formation two-dimensional (2-D) galaxy evolution models with the commonly used one-dimensional (1-D) models with smooth star formation. The 2-D stochastic models predict high variability of the star formation rate and the surface photometric parameters across the galactic disks and in time.

We have developed a fast numerical 2-D model of galaxy disk evolution (resolved along the galaxy radius and azimuth) by adopting a scheme of parameterized stochastic self-propagating star formation. We explore the parameter space of the model and demonstrate its capability to reproduce 1-D radial profiles of the galaxy M33: gas surface density, surface brightness in the i and GALEX FUV passbands, and metallicity.

Stochasticity of bright stars introduces uncertainty and bias into derived structural parameters of star clusters. We have simulated a grid of cluster $V$-band images, observed with Subaru Suprime-Cam with age, mass, and size representing a cluster population in the M31 galaxy and derived their structural parameters by fitting King model to the surface brightness distribution. We have found that clusters less massive than $10^4 M_\odot$ show significant uncertainty in their core and tidal radii for all ages, while clusters younger than 10 Myr have their sizes systematically underestimated for all masses. This emphasizes the importance of stochastic simulations to asses the true uncertainty of structural parameters in studies of semi-resolved and unresolved clusters.

This paper is the second of a series that investigates the stochasticity and degeneracy problems that hinder the derivation of the age, mass, extinction, and metallicity of unresolved star clusters in external galaxies when broad-band photometry is used. While Paper I concentrated on deriving age, mass, and extinction of star clusters for one fixed metallicity, we here derive these parameters in case when metallicity is let free to vary. The results were obtained using several different filter systems ($UBVRI$, $UBVRIJHK$, GALEX+$UBVRI$), which allowed to optimally reduce the different degeneracies between the cluster physical parameters. The age, mass, and extinction of a sample of artificial star clusters were derived by comparing their broad-band integrated magnitudes with the magnitudes of a large grid of cluster models with various metallicities. A large collection of artificial clusters was studied to model the different degeneracies in the age, mass, extinction, and metallicity parameter space when stochasticity is taken into account in the cluster models. We show that, without prior knowledge on the metallicity, the optical bands ($UBVRI$) fail to allow a correct derivation of the age, mass, and extinction because of the strong degeneracies between models of different metallicities. Adding near-infrared information ($UBVRI$+$JHK$) slightly helps in improving the parameter derivation, except for the metallicity. Adding ultraviolet data (GALEX+$UBVRI$) helps significantly in deriving these parameters and allows constraining the metallicity when the photometric errors have a Gaussian distribution with standard deviations 0.05 mag for $UBVRI$ and 0.15 mag for the GALEX bands.

We developed a method for a fast modeling of broad-band UBVRI integrated magnitudes of unresolved star clusters and used it to derive their physical parameters (age, mass, and extinction). The method was applied on M33 galaxy cluster sample and consistency of ages and masses derived from unresolved observations with the values derived from resolved stellar photometry was demonstrated. We found that interstellar extinction causes minor age-extinction degeneracy for the studied sample due to a narrow extinction range in M33.

Context. Stochasticity and physical parameter degeneracy problems complicate the derivation of the parameters (age, mass, and extinction) of unresolved star clusters when using broad-band photometry. Aims. We develop a method to simulate stochasticity and degeneracies, and to investigate their influence on the accuracy of derived physical parameters. Then we apply it to star cluster samples of M31 and M33 galaxies. Methods. Age, mass and extinction of observed star clusters are derived by comparing their broad-band UBVRI integrated magnitudes to the magnitudes of a large grid of star cluster models with fixed metallicity Z=0.008. Masses of stars for a cluster model are randomly sampled from the initial mass function. Models of star clusters from the model grid, which have all of their magnitudes located within 3 observational errors from the magnitudes of the observed cluster, are selected for the computation of its age, mass, and extinction. Results. In the case of the M31 galaxy, the extinction range is wide and the age-extinction degeneracy is strong for a fraction of its clusters. Because of a narrower extinction range, the age-extinction degeneracy is weaker for the M33 clusters. By using artificial cluster sample, we show that age-extinction degeneracy can be reduced significantly if the range of intrinsic extinction within the host galaxy is narrow.

We present a model of chemical and spectrophotometric evolution of disk galaxies based on a stochastic self-propagating star formation scenario. The model incorporates galaxy formation through the process of accretion, chemical and photometric evolution treatment, based on simple stellar populations (SSP), and parameterized gas dynamics inside the model. The model reproduces observational data of a late-type spiral galaxy M33 reasonably well. Promising test results prove the applicability of the model and the adequate accuracy for the interpretation of disk galaxy properties.