The physical origins of the mechanical properties of Fe-rich Si alloys are investigated by combining electronic structure calculations with statistical mechanics means such as the cluster variation method, molecular dynamics simulation, etc, applied to homogeneous and heterogeneous systems. Firstly, we examined the elastic properties based on electronic structure calculations in a homogeneous system and attributed the physical origin of the loss of ductility with increasing Si content to the combined effects of magneto-volume and D0 3 ordering. As a typical example of a heterogeneity forming a microstructure, we focus on grain boundaries, and segregation behavior of Si atoms is studied through high-precision electronic structure calculations. Two kinds of segregation sites are identified: looser and tighter sites. Depending on the site, different segregation mechanisms are revealed. Finally, the dislocation behavior in the Fe-Si alloy is investigated mainly by molecular dynamics simulations combined with electronic structure calculations. The solid-solution hardening and softening are interpreted in terms of two kinds of energy barriers for kink nucleation and migration on a screw dislocation line. Furthermore, the clue to the peculiar work hardening behavior is discussed based on kinetic Monte Carlo simulations by focusing on the preferential selection of slip planes triggered by kink nucleation.npj Computational Materials (2017) 3:10 ; doi:10.1038/s41524-017-0012-4 INTRODUCTION Mechanical properties of alloys originate from the behavior of electrons and atoms on the microscopic scale; however, the strength of atomic bonding does not directly relate to macroscopic mechanical properties such as the yield strength, fracture toughness, and fatigue resistance. This is because the microstructure on the mesoscale, which is composed of various defects such as impurities, grain boundaries, and dislocations, mediates or enhances the response of alloys against external forces, leading to fairly nonlinear and multiscale phenomena. Facing such a nonlinear nature associated with the strength and plastic deformation of alloys, identifying the controlling factors for the mechanical properties is a significantly difficult task, and clarifying the underlying physics at different length and time scales still remains a major challenge in materials science.The authors of the present article took up the challenge of this complicated problem with their particular theoretical and computational tools unique to each characteristic scale range. These tools are electronic structure calculations, the cluster variation method (CVM) of statistical mechanics, and molecular dynamics (MD) simulations, which may cover the atomic to mesoscopic scale ranges. By virtue of high-performance computers, some of the results obtained by the present calculations achieve a higher precision and reveal a clearer picture than those obtained by ordinary experiments.Among the various mechanical properties, the main focus of the present study is the solid-solut...