In this article, we aim to provide a general and complete understanding of semi-supervised (SS) causal inference for treatment effects, using two such estimands as prototype cases. Specifically, we consider estimation of: (a) the average treatment effect and (b) the quantile treatment effect, in an SS setting, which is characterized by two available data sets: (i) a labeled data set of size n, providing observations for a response and a set of potentially high dimensional covariates, as well as a binary treatment indicator; and (ii) an unlabeled data set of size N , much larger than n, but without the response observed. Using these two data sets, we develop a family of SS estimators which are guaranteed to be: (1) more robust and (2) more efficient, than their supervised counterparts based on the the labeled data set only. Moreover, beyond the "standard" double robustness results (in terms of consistency) that can be achieved by supervised methods as well, we further establish root-n consistency and asymptotic normality of our SS estimators whenever the propensity score in the model is correctly specified, without requiring specific forms of the nuisance functions involved. Such an improvement in robustness arises from the use of the massive unlabeled data, so it is generally not attainable in a purely supervised setting. In addition, our estimators are shown to be semiparametrically efficient also as long as all the nuisance functions are correctly specified. Moreover, as an illustration of the nuisance function estimation, we consider inverse-probability-weighting type kernel smoothing estimators involving possibly unknown covariate transformation mechanisms, and establish in high dimensional scenarios novel results on their uniform convergence rates. These results should be of independent interest. Numerical results on both simulated and real data validate the advantage of our methods over their supervised counterparts with respect to both robustness and efficiency.