Planar defects such as free surfaces and grain boundaries in metals and alloys play important roles affecting many material properties such as fracture toughness, corrosion resistance, wetting, and catalysis. Their interactions with point defects and solute elements also play critical roles on governing the microstructural evolution and associated property changes in materials. This work seeks to use atomistic modeling to obtain a fundamental understanding of many surface and interface related properties and phenomena, namely: orientation-dependent surface energy of elemental metals and alloys, segregation of solute elements at grain boundaries and their impact on grain boundary cohesive strength, and the controversial sluggish diffusion in both the bulk and grain boundaries of high entropy alloys. First, an analytical formula is derived, which can predict the surface energy of any arbitrary (h k l) crystallographic orientation in both body-centered-cubic (BCC) and face-centered-cubic (FCC) pure metals, using only two or three low-index (e.g., (100), (110), (111)) surface energies as input. This analytical formula is validated against 4357 independent single element surface energies reported in literature or calculated by the present author, and it proves to be highly accurate but easy to use. This formula is then expanded to include the simple-cubic (SC) structure and tested against 4542 surface energies of metallic alloys of different cubic structures, and good agreement is achieved for most cases. Second, the effect of segregation of substitutional solute elements on grain boundary cohesive strength in BCC Fe is studied. It is found that the bulk substitution energy can be used as an effective indicator to predict the embrittlement or strengthening potency induced by the solute segregation at grain boundaries. Third, the controversial vacancy-mediated sluggish diffusion in an equiatomic FeNiCrCoCu FCC high entropy alloy is studied. Many literature studies have postulated that the compositional complexity in high entropy alloys could lead to sluggish diffusion. To test this hypothesis, this work compares the vacancy-mediated self-diffusion in this model high entropy alloy with a hypothetical single-element material (called average-atom material) that has similar average properties as the high entropy alloy but without the compositional complexity. The results show that the self-diffusivities in the two bulk systems are very similar, suggesting that the compositional complexity in the high entropy alloy may not be sufficient to induce sluggish diffusion in bulk high entropy alloys. Based on the knowledge learned from the bulk alloy, the exploration of the possible sluggish diffusion has been extended to grain boundaries, using a similar approach as in the study of self-diffusion in bulk. Interestingly, the results show that sluggish diffusion is evident at a Σ5(210) grain boundary in the high entropy alloy due to the compositional complexity, especially in the low temperature regime, which is different from the bulk diffusion. The underlying mechanisms for the sluggish diffusion at this grain boundary is discussed. / Doctor of Philosophy / Human beings have utilized metals and alloys for over ten millennia and learned much from them. Based on the accumulated knowledge, they have countless applications in our current daily life. However, there is still much to learn for improving our current technology and even opening new opportunities. Throughout most of history, our understanding of these materials was largely obtained through empirical experimentation and refining them into theories and scientific laws. Nowadays, due to the advancements in computer simulations, we can learn more by modeling the behaviors of metals and alloys at the length and time scales that are either be too arduous, costly, or currently impossible experimentally. This work aims at using computer modeling to study some important surface/interface related physical behaviors and properties in metals and alloys at the atomistic scale. First, this work intends to develop a robust surface energy model in an analytical form for any crystallographic orientation. Surface energy is an important material property for many surface-related processes such as fracturing, wetting, sintering, catalysis, and crystalline particle shape. Surface energy is different at different surface orientations, and predicting this difference is important for understanding these surface phenomena. Second, the effect of solute segregation on grain boundary cohesive strength is studied. Most commonly used metallic materials consist of many small crystalline grains and the borders between them are called grain boundaries, which are weak spots for fracture. The minimum energy required to split a boundary is called the grain boundary cohesive strength. The presence of solutes or impurities at grain boundaries can further alter the cohesive strength. A better understanding of this phenomena will eventually help us develop more fracture-resistant materials. The third project deals with the possible sluggish/retarded diffusion in high entropy alloys, which contain five or more principal alloying elements and have many unique mechanical, radiation-resistant, and corrosion-resistant properties. Many researchers attribute these unique properties to the slow species diffusion in these alloys, but its existence is still controversial. This work studies the atomic-level diffusion mechanisms in an FeNiCrCoCu high entropy alloy both in bulk (grain interior) and at grain boundaries in order to determine if sluggish diffusion is present and its causes.