The Model of Aerosol Dynamics, Reaction, Ionization and Dissolution (MADRID) with three improved gas/particle mass transfer approaches (i.e., bulk equilibrium (EQUI), hybrid (HYBR), and kinetic (KINE)) has been incorporated into the Weather Research and Forecast/Chemistry Model (WRF/Chem) (referred to as WRF/Chem‐MADRID) and evaluated with a 5‐day episode from the 2000 Texas Air Quality Study (TexAQS2000). WRF/Chem‐MADRID demonstrates an overall good skill in simulating surface/aloft meteorological parameters and chemical concentrations of O3 and PM2.5, tropospheric O3 residuals, and aerosol optical depths. The discrepancies can be attributed to inaccuracies in meteorological predictions (e.g., overprediction in mid‐day boundary layer height), sensitivity to meteorological schemes used (e.g., boundary layer and land‐surface schemes), inaccurate total emissions or their hourly variations (e.g., HCHO, olefins, other inorganic aerosols) or uncounted wildfire emissions, uncertainties in initial and boundary conditions for some species (e.g., other inorganic aerosols, CO, and O3) at surface and aloft, and some missing/inactivated model treatments for this application (e.g., chlorine chemistry and secondary organic aerosol formation). Major differences in the results among the three gas/particle mass transfer approaches occur over coastal areas, where EQUI predicts higher PM2.5 than HYBR and KINE due to improperly redistributing condensed nitrate from chloride depletion process to fine PM mode. The net direct, semi‐direct, and indirect effects of PM2.5 decrease domainwide shortwave radiation by 11.2–14.4 W m−2 (or 4.1–5.6%) and near‐surface temperature by 0.06–0.14°C (or 0.2–0.4%), lead to 125 to 796 cm−3 cloud condensation nuclei at a supersaturation of 0.1%, produce cloud droplet numbers as high as 2064 cm−3, and reduce domainwide mean precipitation by 0.22–0.59 mm day−1.